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j/j ' Abstract. In recent decades, novel magnetism of d- and /-electron compounds has 

been discussed very intensively both in experimental and theoretical research fields of 
condensed matter physics. It has been recognized that those material groups are in 
the same category of strongly correlated electron systems, while the low-energy physics 
'"^ , of d- and /-electron compounds has been separately investigated rather in different 

C ' manners. One of common features of both d- and /-electron systems is certainly 

the existence of active orbital degree of freedom, but in /-electron materials, due to 
the strong spin-orbit interaction in rare-earth and actinide ions, the physics seems 
to be quite different from that of d-electron systems. In general, when the number 

^ , of internal degrees of freedom and relevant interactions is increased, it is possible to 

C^^ ' obtain rich phase diagram including large varieties of magnetic phases by using several 

^y ■ kinds of theoretical techniques. However, we should not be simply satisfied with the 

^b , reproduction of rich phase diagram. It is believed that more essential point is to seek 

^^ ' for a simple principle penetrating complicated phenomena in common with d- and /- 

electron materials, which opens the door to a new stage in orbital physics. In this sense, 
it is considered to be an important task of this article to explain common features of 

j^ ' magnetism in d- and /-electron systems from a microscopic viewpoint, using a key 

concept of orbital ordering, in addition to the review of the complex phase diagram 
of each material group. As a typical d-electron complex material exhibiting orbital 
order, first we focus on perovskite manganites, in which remarkable colossal magneto- 

Q [ resistance effect has been intensively studied. The manganites provide us a good stage 

O ' to understand that a simple mechanism works for the formation of complex spin, 

charge, and orbital ordering. We also explain intriguing striped charge ordering on the 
orbital-ordered background in nickelates and the effect of orbital ordering to resolve 

; ^ ' spin frustration in geometrically frustrated eg electron systems. Note that orbital 

5^ , ordering phenomena are also found in i2g electron systems. Here we review recent 

advances in the imderstanding of orbital ordering phenomenon in Ca2Ru04. Next 
we discuss another spin-charge-orbital complex system such as /-electron compound. 
After the detailed explanation of the construction of microscopic models on the basis 
of a j-j coupling scheme, we introduce a d-electron-like scenario to understand novel 
magnetism in some actinide compounds with the HoCoGa5-type tetragonal crystal 
structure. Finally, we show that complicated multipole order can be understood from 
the spin-orbital model on the basis of the j-j coupling scheme. As a typical material 
with multipole order, we pick up Np02 which has been believed to exhibit peculiar 
octupole order. Throughout this review, it is emphasized that the same orbital physics 
works both in d- and /-electron complex materials in spite of the difference between d 
and / orbitals. 
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1. Introduction 

It has been widely recognized that orbital degree of freedom plays a key role in 
understanding of novel magnetism observed in transition metal oxides PQ |21 El El- A 
typical material of such spin-charge-orbital complex is the manganese oxide, exhibiting 
remarkable colossal magneto-resistance (CMR) phenomena 0. In the recent decade, the 
study of manganites has been one of the most important areas of research in condensed 
matter physics. In one word, the CMR effect is considered to occur when the manganite 
ground-state changes from insulating to ferromagnetic (FM) metallic, after a small 
magnetic field is applied. Based on the concept of two-phase competition, the CMR 
behavior has been successfully qualitatively reproduced in computational simulations, 
for instance, employing resistor- network models [Hj. In the two phases, the appearance 
of the FM metallic phase in manganites has been usually rationalized by the so-called 
double-exchange (DE) mechanism [Zj, based on a strong Hund's rule coupling between 
mobile Cg electrons and localized t2g spins. On the other hand, the insulating phase in 
manganites is basically understood by the coupling between degenerate Cg electrons and 
Jahn- Teller (JT) distortions of the MnOe octahedra 011], leading to the various types 
of charge and/or orbital orders observed experimentally. 

The rich phase diagram of manganites has been revealed due to competition and 
interplay among spin, charge, and orbital degrees of freedom, but a recent trend is 
to unveil further new phases both from experimental and theoretical investigations. 
A typical example can be found in the undoped perovskite manganite, RMnOa with 
rare earth ion R, which is the mother compound of CMR manganites. For R=La, 
it has been understood clearly that the A-type antiferromagnetic (AF) phase appears 
JHIE] with the C-type ordering of (3x^— r^)- and (3y^— r^)-orbitals jTUI. Here "A-type" 
denotes a layered antiferro structure with ferro-order in the ab plane and antiferro- 
order along the c axis, while "C-type" indicates a chain-type antiferro structure with 
antiferro-order in the ab plane and ferro-order along the c axis jH] . See Fig. IHl^a) for 
each structure. Theoretically, the A-type ordering has been explained by using several 
kinds of techniques [ni[iaiISllIllliaiIHllIIllIHllinil2Dlinil2ai2nil21l2ai2Ell2II- By 
substituting La by alkaline earth ions such as Sr and Ca, holes are effectively doped into 
Cg-electron band and due to the DE mechanism, the FM metallic phase appears with 
its concomitant CMR effect. Most of the discussion in manganites has centered on the 
many phases induced by doping with holes the A-type AF state, at different values of 
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their bandwidths. In this framework, it is imphcitly assumed that the undoped material 
is always in the A-type AF state. 

However, recently, a new AF phase has been reported as the ground state in the 
undoped limit for R=Ho [2HI 121] • This phase is called the "E-type" spin structure 
following the standard notation in this context |H]. See Fig. Mjo) for the E-type spin 
structure. It is surprising that a new phase can be still found even in the undoped 
material, previously considered to be well understood. In addition, the nature of the 
states obtained by lightly doping this E-phase is totally unknown, and new phenomena 
may be unveiled experimentally in the near future. This is believed to open an exciting 
new branch of investigations in manganites jHOl i^, since novel phases appear to be 
hidden in the vast parameter space of these compounds. A clear example has been 
recently provided by the prediction of a FM charge-ordered (CO) phase at x=l/2 |32tl33]. 
which may have been found experimentally jSHinS]- These facts indicate the importance 
of both experimental and theoretical efforts to unveil new phases in manganites, in 
addition to the explanation of the complex phases already observed. Such efforts have 
also been made to find new phases in other transition metal oxides, for instance, 
ruthenates and nickelates, as we will see later in this article. Concerning RMnOs 
with hexagonal structure, quite recently, "multiferroics" has been another keyword to 
understand exotic magnetic phenomena emerging from the multi-phase competition |36j . 
We believe that it is useful to review the nature of spin, charge, and orbital ordered 
phases of manganites and other transition metal oxides from a unified viewpoint, even 
though it is true that more work remains to be done to fully understand transition 
metal oxides, in particular, unusual magneto-transport properties of manganese oxides 
and appearance of unconventional superconductivity. 

A trend to seek for new magnetic as well as superconducting phases has been also 
found in the /-electron system, which is another type of spin-charge-orbital complex 
|37| EH]- Among so many kinds of /-electron materials, in recent years, /-electron 
compounds with IIoCoGa5-type tetragonal crystal structure, frequently referred to as 
"115", have been intensively investigated both in experimental and theoretical research 
fields of condensed matter physics. Such vigorous activities are certainly motivated 
by "high" temperature superconductivity observed in some 115 compounds. First, 
unconventional superconductivity has been found in Ce-based 115 compounds, CeTIns 
(T=Rh, Ir, and Co). A surprising point is that CeCoIus exhibits the superconducting 
transition temperature Tc=2.3K [nHj, which was the highest among yet observed for 
heavy fermion materials at ambient pressure when it was discovered. On the other 
hand, Celrlns shows rc=0.4K jlH! which is much less than that of CeCoIus. Note that 
CeRhlus is antiferromagnet with a Neel temperature Tn=3.8K at ambient pressure, 
while under high pressure, it becomes superconducting with Tc=2.1K ^^. 

After the discovery of superconductivity in Ce-115, the rapid expansion of the 
research frontier to transuranium systems has been accelerated by the discovery of 
superconductivity of Pu-based 115 compounds, PuTGas (T=Co and Rh). It has been 
reported that Tc of PuCoGas is 18. 5K |121 EH]; which is amazingly high value even 
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compared with other well-known intermetallic compounds. The coefficient of electronic 
specific heat 7 is estimated as 7=77mJ/mol-K^, moderately enhanced relative to that 
for normal metals, suggesting that PuCoGas should be heavy-fermion superconductor. 
In PuRhGas, superconductivity has been also found [Sj. Although the value of 
Tc=8.7K is lower than that of PuCoGa5, it is still high enough compared with other 
heavy-fermion superconductors. Quite recently, high quality single crystal PuRhGas 
has been synthesized jl3] and the Ga-NQR measurement has revealed that d-wave 
superconductivity is realized in PuRhGas jlH]. The Ga-NMR measurement of PuCoGas 
is consistent with this conclusion ^W,. PuIrGas has been also synthesized, but it is 
considered to be paramagnetic at ambient pressure. At least up to now, there is no 
indication of superconductivity even under the pressure of 9.5GPa down to 1.4K [15] . 

Besides such high temperature superconductivity of Ce-115 and Pu-115 compounds, 
interesting magnetic properties have been reported for UTGas, where T is a transition 
metal ion. For several transition metal ions T, UTGas are AF metals or Pauli 
paramagnets |iniEniEIlE2lEniEllESlEniEZlEHlEnilSn!. Among them, neutron 
scattering experiments have revealed that UNiGas exhibits the G-type AF phase, while 
UPdGas and UPtGas have the A-type AF state jHElEn]- See Fig. ESI for the magnetic 
structure. Note that G-type indicates a three-dimensional Neel state 0. On the other 
hand, for UTGas with T=Co, Rh, Ir, Fe, Ru, and Os, magnetic susceptibility is almost 
independent of temperature, since these are Pauli paramagnets. It is quite interesting 
that the magnetic structure is different for U-115 compounds which differ only by the 
substitution of transition metal ions. 

Recently, Np-115 compounds NpTGas (T=Fe, Co, Rh, and Ni) have been also 
synthesized and several kinds of physical quantities have been successfully measured 
inHinainilMlinSlinniEIlinHlinnilZni- in particular, the de Haas- van Alphen (dHvA) 
effect has been observed in NpNiGas ^^, which is the first observation of dHvA signal in 
transuranium compounds. Quite recently, dHvA measurement has been also successfully 
performed in plutonium compound Pulua over beyond several kinds of difficulties |71j . 
For NpCoGas and NpRhGas, the dHvA oscillations have been also detected and plural 
number of cylindrical Fermi surfaces are found [H^inH]- For NpFeGas, the magnetic 
moment at Fe site has been suggested in neutron scattering experiments J66.] and it 
has been also detected by ^^Fe Mossbauer spectroscopy [HHl- The magnetic structure of 
Np-115 compounds also depends sensitively on transition metal ion JUUIIH^ 70j: G-AF 
for NpNiGas, A-AF for NpCoGas and NpRhGas, and C-AF for NpFeGas. See Fig. EBl 
for the magnetic structure. Note also that in the neutron scattering experiment for 
NpNiGas, the signal suggesting the G-AF grows, after the FM transition occurs |66j . 
This G-AF structure is due to canted magnetic moments of Np ions. It is characteristic 
of U-115 and Np-115 compounds that the magnetic properties are sensitive to the choice 
of transition metal ions. 

The appearance of several kinds of AF states in U-115 and Np-115 compounds 
reminds us of the magnetic phase diagram of CMR manganites. Thus, we envisage a 
scenario to understand the complex magnetic structure of actinide compounds based on 
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an orbital degenerate model similar to that of manganites. However, one must pay due 
attention to the meanings of "spin" and "orbital" in /-electron systems. Since they are 
tightly coupled with each other through a strong spin-orbit interaction, distinguishing 
them is not straightforward in comparison with rf-electron systems. This point can create 
serious problems when we attempt to understand microscopic aspects of magnetism and 
superconductivity in /-electron compounds. Thus, it is necessary to carefully define the 
terms "orbital" and "spin" for / electrons in a microscopic discussion of magnetism and 
superconductivity in /-electron compounds. 

From a conceptual viewpoint, in general, / electrons are more localized in 
comparison with d electrons, but when we turn our attention from 4/ to 5/, electronic 
properties are changed gradually from localized to itinerant nature, leading to rich 
phenomena which have been recently investigated intensively. In this sense, it has been 
also highly requested to push forward the microscopic research on / electron systems. 
However, in sharp contrast to d electron systems, the existence of strong spin-orbit 
coupling has been a problem to develop a theoretical study in the same level as d- 
electron research. 

In order to overcome such problems, it has been recently proposed to employ a j-j 
coupling scheme to discuss microscopic aspects of magnetism and superconductivity of 
/-electron systems j721 ES], on the basis of the relativistic band-structure calculation 
results [ZmZSllZnilZZllZHllZOlEIIllHIl. There are a couple of advantages of the j-j 
coupling scheme. First, it is quite convenient for the inclusion of many-body effects 
using standard quantum-field theoretical techniques, since individual /-electron states 
are clearly defined. In contrast, in the LS coupling scheme we cannot use such standard 
techniques, since Wick's theorem does not hold. Second we can, in principle, include 
the effects of valence fluctuations. In some uranium compounds, the valence of the 
uranium ion is neither definitely U^"^ nor U^"'", indicating that the /-electron number 
takes a value between 2 and 3. In the j-j coupling scheme this is simply regarded as 
the average number of / electron per uranium ion. 

As we will discuss later in detail, with the use of the j-j coupling scheme, it is 
possible to establish the microscopic Hamiltonian for /-electron systems. In particular, 
under crystalline electric field of cubic symmetry, the /-electron model on the basis 
of the j-j coupling scheme can be reduced to the orbital degenerate model, which is 
equivalent to the Cg electron Hubbard model [72]. The common microscopic model for 
d- and /-electron compounds is a simple explanation for the appearance of magnetic 
structure in common with manganites and actinide 115 compounds, as actually analyzed 
by using numerical calculations jH2l IHHj- In order to understand unconventional 
superconductivity of Ce-115 and Pu-115 materials, the orbital degenerate model has 
been also analyzed with the use of a fluctuation-exchange approximation to include 
effectively spin and orbital fluctuations [HH ESI IHE] ■ A possible scenario for odd-parity 
triplet superconductivity induced by the Hund's rule interaction has been discussed 
based on the orbital degenerate model on the non-Bravais lattice such as honeycomb 
lattice jHZl IHHl- Novel magnetism and exotic superconductivity of filled skutterudite 
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materials (HHl ED] have been discussed by using the microscopic model on the basis of 
the j-j coupling scheme pT| 1^ IHSl HI] • We would like to emphasize that the same 
orbital physics should work between d- and /-electron systems, after the application of 
the j-j coupling scheme to /-electron materials. This is an important clue to establish 
the unified picture, which penetrates complex phenomena both in d- and /-electron 
compounds. 

Another advantage of the model on the basis of the j-j coupling is that it is 
possible to develop a microscopic theory for multipole ordering of /-electron materials. 
Recently, ordering of high-order multipole such as octupole has been intensively 
discussed for Ce^Lai-^Bg [HSl EEl EH EHl EHl IMIl EIIl IIII21 CUSl EH and for NpOa 
[inSllinnilinillinHllinillinillllllliailll to reconcile experimental observations which 
seem to contradict one another at first glance. Very recently, a possibility of octupole 
ordering has been also proposed for SmRu4Pi2 [ IT^ 1115] . Here we note that spin and 
orbital degrees of freedom correspond to dipole and quadrupole moment, respectively. 
The microscopic aspects of octupole ordering has not been discussed satisfactorily in the 
context of electronic models. Rather, phenomenological theories have been developed 
under the assumption that octupole ordering occurs. Note that direct detection of 
octupole ordering is very difficult, since the octupole moment directly couples to neither 
a magnetic field nor lattice distortions. However, those phenomenological theories have 
successfully explained several experimental facts consistently, e.g., induced quadrupole 
moments in octupole ordered states in Cea;Lai_a;B6 pUlllUUtllUlj and Np02 |107tlllUJ . 
On the other hand, it has been highly required to proceed to microscopic theory, in 
order to understand the origin of multipole ordering in /-electron systems over beyond 
the phenomenological level. 

Concerning this issue, recently, it has been clearly shown that the model based on 
the j-j coupling scheme also works for the explanation of octupole ordering |113t I116| 
I117[ I118[ I119J . It is possible to obtain the multipole interaction from the /-electron 
model, by applying the same procedure to derive the orbital-dependent superexchange 
interaction in d-electron systems. Namely, we can provide the microscopic basis for 
multipole interaction. It is also possible to show the stability of octupole ordered phase 
depending on the lattice structure. In this scenario, octupole is the combined degree of 
freedom of spin and orbital. In another word, octupole is considered to be characterized 
by the anisotropic spin densities. Thus, the difference in anisotropic up- and down-spin 
densities naturally provide the magnetic moment of octupole. It is one of progresses 
in orbital physics that higher-order multipole can be also understood in the context of 
spin and orbital ordering phenomena. 

In this article, in Sec. 2, we review the construction of the microscopic models for d- 
electron systems in detail. Then, we arrive at three kinds of Hamiltonians, depending on 
the active orbitals, which are (i) Cg-electron doubly degenerate model, (ii) t2g-electron 
triply degenerate model, and (ni) Cg-orbital double-exchange model coupled with t2g 
localized spins. In Sec. 3, in order to explain the orbital ordering in manganites, we 
focus on the theoretical results on the model (iii). In particular, the topological aspects 
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of orbital ordering will be discussed in detail. Then, we explain characteristic features 
of orbital ordering depending on hole doping x in subsections for x=0.0, x=0.5, x>0.5, 
and x<0.5. In Sec. 4, we discuss the results of models (i) and (ii), by introducing 
nickelates and ruthenates as typical materials, respectively. In Sec. 5, we will move 
on to the model construction for /-electron materials with the use of the j-j coupling 
scheme. We will explain the spirit of the j-j coupling scheme in detail, in comparison 
with the LS coupling scheme. Then, we establish the orbital degenerate model for /- 
electron systems. In Sec. 6, we will show theoretical results on the microscopic model 
for /-electron systems. We will explain possible orbital ordering of 115 compounds and 
octupole ordering of Np02 from a microscopic viewpoint. Finally, in Sec. 7, we will 
summarize this article. 

2. Model Hamiltonian for rf-electron systems 

Before proceeding to the description of theoretical results on orbital ordering phenomena 
in d-electron compounds, first it is necessary to define the model Hamiltonian for d- 
electron systems. The model should be composed of three parts as 

-n = -nkin + -nioc + -nintcr-sitc; (1) 

where ii^kin expresses the kinetic term of d electrons, iJioc denotes the local term for d 
electrons, and ifinter-site indicates the inter-site interaction among d electrons, which is 
not fully included by the combination of i^km and i/ioc- Among them, the local term 
H\oc includes three important ingredients, written as 

H\oc = -f^CEF + Hc\-c\ + -f^cl-ph) (2) 

where Hqef is the crystalline electric field (CEF) term, -f/^ci-ci denotes the Coulomb 
interactions among d electrons, and -f^ci-ph indicates the coupling term between d 
electrons and lattice distortions. The full Hamiltonian includes several competing 
tendencies and couplings, but as shown below, the essential physics can be obtained 
using relatively simple models, deduced from the complicated full Hamiltonian. 

2.1. Crystalline electric field effect 

In order to construct the model Hamiltonian for rf-electron systems, let us start our 
discussion at the level of the atomic problem, in which just one electron occupies a 
certain orbital in the "id shell of a transition metal ion. In the next subsection, we will 
include the effect of Coulomb interactions among d electrons. For an isolated ion, a 
five-fold degeneracy exists for the occupation of the 3(i orbitals, but this degeneracy 
is partially lifted by the CEF potential from anions surrounding the transition metal 
ion. Since it is the electrostatic potential field acting on one electron state even in the 
complicated crystal structure, the CEF potential should be, in any case, given in the 
second-quantized form as 

-f^CEF = 2_^ ^m,m"^imo-C^im'o-5 (3) 

i,(T,m,,m' 
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where dima is the annihilation operator for a d-electron with spin a in the m-orbital 
at site i, m(=— 2, ■ ■ ■ , 2) is the z-component of angular momentum i{=2), and Ajn,m' is 
the coefficient of the CEF potential, depending on the crystal structure and the angular 
momentum. 

The explicit form of Am,m' has been analyzed in detail by the ligand field theory. 
Here we briefiy explain the derivation of Am,m' by following the procedure of Hutchings 
|12U1 I121J . The CEF potential is given by the sum of electro-static potential from anions 
surrounding the transition metal ion, written by 

^CEF(r) = $:^— -, (4) 

where Z is the valence of anion, e is elementary electric charge, Rj denotes the position 
of i-th anion, and r indicates the position of 3d electron around the nucleus of transition 
metal ion. Then, the CEF coefficient is evaluated by 



Am,m' = / rfr^;^„(r)\/cEF(r)^„f^'(r), (5) 

where \&„fm(i") is the wavefunction of d electron, expressed as 

^nimir) = Rneir)Ye^{e, cp), (6) 

with r={r,6,ip) in the polar coordinate. Here Rne{r) is the radial wavefunction and 
Yim is the spherical harmonics, which are obtained by solving the Schrodinger equation 
for the hydrogen-like potential problem. Note that n denotes the principal quantum 
number. In the actual situation, |r| is in the order of Bohr radius, while |Ri| is related 
to the lattice constant which is in the order of several angstroms. Thus, it is convenient 
to expand Vcef in terms of |r|/|Rj| as 

vcEF(r) = E:^E E whi^) yueMyue..^.). (7) 

i k=0 m=—k 

where Hi={Ri,9i,ipi) in the polar coordinate. 

For further calculations, it is necessary to set the actual crystal structure. As a 
typical example, in this review article, we pick up the perovskite structure composed of 
MOg octahedron, in which transition metal ion M is surrounded by six oxygen ions, as 
shown in Fig. ^ To make the situation general, we set the different values, a and b, for 
the size of octahedron in the xy plane and along the z axis, respectively. The positions 
of oxygen ions are, then, given by Ri=(a/2, 0, 0), R2=(0, — a/2, 0), R3=(— a/2, 0, 0), 
R4=(0, a/2,0), R5=(0, 0,6/2), and R6=(0, 0, -6/2). For the case of a = 6, the cubic 
symmetry is maintained, while for the case of a 7^ 6, the system is in the tetragonal 
symmetry. After some algebraic calculations using the explicit form of the wavefunction 
for d electrons, a general form for Am-m' is given by 

a^ c 
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Figure 1. MOg octahedron at site i. The position of the i-th oxygen ion labeled 
by Oi is R,;. The size of the octahedron is specified by a and b, which are given by 
a=|Ri — R3| = |R2 — R4I and 6=|R5 — Re], respectively. See the maintext for the 
definitions of R, . 



(r^) Ze^ /35 

H : \l —c''^\2m,2m')6m-m',±4, 



a^ a 



where 6m,m' is the Kronecker's deha, (r^) is given by 

{r')= / r'=|i?3d(r)|Vrfr, 
Jo 

and c*^'^-' {im, i'm') is the so-called Gaunt coefficient 



(8) 



(9) 



Att 



defined by 

sin ededifY;^{e, <^)n„_„.(^, <^)f,w(^, <^). (lo) 



2k + 1, 

The non-zero values for the Gaunt coefficients have been tabulated in the standard 
textbook |124j . By consulting with the table for the Gaunt coefficient, we explicitly 
obtain Am m' as 



A 



A. 



2,2 — ^-2,-2 — 2A20 + A40, 
A,i = ^-1,-1 = -^20 -4^40, 
Ao,o = -2A20 + 6A40, 

A2-2 = A_20 = A44, 



:ii) 



where 



A20 
A40 
A44 



2 Ze^ (r2) 

7 a a^ 
lZe^{r^)l3 4/aN5 



b 
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Figure 2. Views for d-electron orbitals. 



Note that ^20=0 and ^44=5^40 for the cubic case {a=b). 

For instructive purpose, here we have shown a part of tedious calculations to 
determine CEF potential, but for actual purpose, it is more convenient to consult 
with the table of Hutchings, in which all possible CEF parameters have been listed 
for arbitrary angular momentum J (both integer and half-integer). For the c?-electron 
case, we can simply refer the results for J=2 in the Hutchings table |12Uj . For instance, 
in the tetragonal CEF as discussed above, we easily obtain 



Ai,i = A_i,_i = -3fi° - 48S? 



12,2 - 
,1 = 
^0,0 = 
^2,-2 



4; 



-65° - 
= ^-2,2 



72Bl 
= 12BI 



(13) 



where BJ^ with integers n and m is the so-called CEF parameter, traditionally used 
in the ligand field theory. It is easy to express the CEF parameters 5^ by using 
ours as B^ = A20/3, B^ = A^q/U, and Bf = Am/12. Note again that B^=0 and 
Bf=5B^ for the cubic case. In actuality, we do not estimate the CEF parameters purely 
theoretically, but they are determined from the fitting of experimental results such as 
magnetic susceptibility and specific heat at high temperatures. 

For a cubic case (a=6), the five-fold degeneracy in d orbital is lifted into doubly- 
degenerate Cg-orbitals {dx2_y2 and d^z^-r^) and triply-degenerate t2g-orbitals {dxy, dyz, 
and dzx)- The shape of each orbital is illustrated in Fig. |21 The eigenenergy for Cg- 
orbitals is given by 



£:322_j,2 = ex'2-y2 = 6A^Q = 
while for t2g orbitals, we obtain 

^xy ^yz ^zx ^^40 

Then, the CEF term can be written as 



a a^ 



-48B° = ^. 



3 a a^ 



"I'yai 



(14) 



(15) 



(16) 



i,o-,7 



Orbital ordering phenomena in d- and f -electron systems 



12 



where di^o- is the annihilation operator for a rf-electron with spin a in the 7-orbital at 
site i. The energy difference between those two levels is usually expressed as lODq in 
the traditional notation in the ligand field theory. It is explicitly written as 

5 Ze^ (r^) 

3 a a'^ 



WDq = lOA 



40 



120s? 



:i7) 



Let us try to estimate theoretically the value of WDq. For the purpose, it is convenient 
to transform the above equation into 



lODq 



lOZ 



Ob 
a 



{r') 



(18) 



where ob denotes the Bohr radius (aB=0.529A). Note that Z=2 for 0^~ ion and 
e^/(2aB)=l Ryd.=13.6 eV. By using the solution of the hydrogen-like potential problem, 
we obtain (r^)/a^=25515/Z^, where Zm denotes the atomic number of transition metal 
atom. For instance, if we simply set Zjv/=25 for manganese and a=2A as a typical value 
for perovskite structure, we obtain 10Dg=7.7meV, which is very small compared with 
the observed value, since lODq is found to be in the order of eV in actual transition metal 
oxides. In fact, the estimation by Yoshida suggests that lODq is about 10000-15000cm~^ 
(remember that 1 eV = 8063 cm^^) [125J. The electrostatic contribution obtained in 
the above naive estimation is considered to be much smaller than the observed CEF 
splitting. 

Here note that the energy level for the t2g-orbitals is lower than that for Cg-orbitals 
for perovskite structure. Qualitatively this can be intuitively understood as follows: 
The energy difference originates in the Coulomb interaction between the 3d electrons 
and the oxygen ions surrounding transition metal ion. As shown in Fig. El while the 
wave-functions of the Cg-orbitals are extended along the direction of the bond between 
transition metal ion and oxygen ions, those in the t2g-orbitals avoid this direction. Thus, 
an electron in t2g-orbitals is not heavily influenced by the Coulomb repulsion due to the 
negatively charged oxygen ions, and the energy level for t2g-orbitals is lower than that 
for Cg-orbitals. 

For a tetragonal case with a < 6, as observed in cuprates, we find the splitting of 



each degenerate orbital. The Cg orbital is spht into aig {3z^ 
orbitals, of which eigen energies are, respectively, given by 



and 6ig 



[X 



£3z2-r^ = QDq - 



4 Ze^ (r2) 



7 



a a 



1 



5 Ze^ (r^) 



21 



a a 



1- 



and 



6Z,,+ i^»^<'-^> 



7 



a a^ 



+ 



5 Ze^ (r^) 



21 



a a^ 



where Dq is defined as 

1 Ze^ (r^) 
6 a a"' 



Dq 



2 
3^3 



3U/ 



On the other hand, t2g orbital is split into Cg degenerate [yz and zx) and b2g 
orbitals, of which eigen energies are, respectively, given by 



-yz 



-ADq- 



2 Ze^ (r^) 
7 a d^ 



1- 



+ 



10 Ze^ (r^) 
63 a a^ 



-yl 

(19) 
(20) 

(21) 
{xy) 

(22) 
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and 



e^y = -ADq + 



"^ ^ ' 7 a a2 



'-'b 



2 /»,4\ r 



63 a a^ 



' b 



a\5l 



(23) 



Note that also for the tetragonal case, the CEF term is written as Eq. (fTT)|) . 

The splitting in the tetragonal case with a < 6 is intuitively understood from 
the shape of orbitals, on the basis of the point that the energy difference originates 
in the Coulomb interaction between the 3d electrons and the oxygen ions surrounding 
transition metal ion. While the wave- function of the Sz'^—r'^ orbital is extended along 
the z-direction, x^—y"^ orbital shape is extended in the x-y plane. Since apical oxygens 
move to the z-direction for the case of a < 6, the energy loss for 3^^— r^ orbital becomes 
small, indicating that the 32;^— r^ orbital is lower. Concerning the splitting of t2g orbitals, 
we note that xy orbital is extended in the xy plane, while yz and zx orbitals have some 
extension along z axis. Then, for a < b, the penalty from the electrostatic potential 
would be smaller for yz and zx, compared with that for xy. Then, we intuitively consider 
that the energy level for xy orbital is higher than those for yz and zx orbitals. 

2.2. Coulomb interactions 

Now we include the effect of Coulomb interactions among d electrons in the level of 
atomic problem. In the localized ion system, the Coulomb interaction term among 
d-electrons is generally given by 

-"el-el — T; / , / , / _, ^mim2.mam4^irma,"'im9.(79 i™'30'2^im4,cTiy l24j 

where the Coulomb matrix element among d electrons /"' is expressed as 

j J dVidV2%2m, (ri)*327„, (r2)^12^327„3 (r2)^327774 (n)- (25) 

Here gi2=g{\Yi — r2|) is the screened Coulomb interaction, which can be expanded in 
spherical harmonics 

912 = g{\ri - r2\) =^gt{,ri,r2)Yi,^{ei,(t)i)Yi^{e2,(t)2), (26) 

with ri=(ri, 6*1, 0i) and r2=(r2, 6*2, ^2) in the polar coordinate. The complicated 
integrals can be partly performed and the result is given by using the Gaunt coefficients 

as 

^mi7772,7n3m4 = '^'ni+'n2,m3+m4 ^l C^^'H^^l, "^3)c^''^("^4, ^2)^^=. (27) 

A:=0,2,4 

Note that the sum is limited by the Wigner-Eckart theorem to A;=0, 2, and 4. Here we 
define F^, which is the radial integral for the fc-th partial wave, called Slater integral or 
Slater-Condon parameter |126tll27j . given by 

F^d= rldr, rldr2RUn)RUr2)gkin,r2). (28) 

Jo Jo 



77717772,77137774 
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Figure 3. Configurations of d" electrons for n=1^9. For n—4^7, we show both the 
high- and low-spin states. 



By using Slater-Condon parameters, it is more convenient to define Racah parameters, 
given by 

B=i9Fi-5F^)/AAl, (29) 

C = 5F^763, 

where A, B, and C are Racah parameters |12Hj . All possible Coulomb integrals 
are expressed by these Racah parameters |128j . In usual, the Racah parameters are 
determined from the experimental results in the high-energy region. 

By diagonalizing HcEF+He\-e\, let us discuss rf^-electron configuration under the 
cubic CEF, where the superscript n denotes the local d-electron number. It is possible to 
obtain analytic result partly, but most of the solutions can be obtained only numerically. 
As summarized in Fig. 01 the results are understood from the electron configuration 
obtained by accommodating n electrons in Cg and t2g orbitals due to the competition 
between the CEF splitting lODq and the Hund's rule interaction. For n=l, we find 
6-fold degenerate ground state, originating from the three t2g orbitals with spin degree 
of freedom. For n=2, 9- fold degeneracy is found in the ground state. When we 
accommodate a couple of electrons among three t2g orbitals, there are three possibilities 
for the formation of total spin 5=1 between two of three orbitals. Since each spin S=l 
has three degenerate states with 5*2=— 1, 0, and -|-1, in total we obtain 9-fold degeneracy. 
For n=3, 4- fold degeneracy is found in the ground state. Now three t2g electrons form 
S=3/2 spin due to Hund's rule, which have four degenerate states with ^^=±3/2 and 
±1/2. 

For the case of n=4, we find two types of ground states, depending on the balance 
between lODq and Coulomb interaction. For small lODq, 10-fold degeneracy is observed 
in the ground state, while the degeneracy is changed to 9-fold for large lODq. These 
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are high- and low-spin states, labeled by "H" and "L", respectively. In the high-spin 
configuration, due to the strong Hund's rule interaction, total spin S=2 is formed, but we 
consider that one electron is added to the d^ configuration. The fourth electron should 
be accommodated in Cg sector, which provide an extra double degeneracy. Since S=2 
state include five degenerate states with Sz=±2, ±1, and 0, in total we obtain 10-fold 
degeneracy. The eigen energy for (i^(H) is exactly given by £^[(i^(H)]=674 — 215 — 12B^. 
On the other hand, for the low-spin state, we consider that all four electrons are included 
among t2g orbitals to gain the CEF potential energy. Namely, S=l spin is formed in the 
t2g sector, leading to 9-fold degeneracy as in the case of d'^. The energy for d'^ih) can be 
obtained only numerically, but in the limit of infinite lODg, the energy is asymptotically 
given by E[d'^{L)]=QA — 15B + 5C — 192i?4. Due to the comparison of two energies, we 
obtain a rough condition for the change of the spin state as 6B + 5C > lODq, which 
indicates clearly that the high-spin state is chosen when the Hund's rule interaction is 
stronger than the effect of CEF potential. 

In the high-spin state for n=4, as observed in Mn^+ ion, three electrons are 
accommodated in the t2g orbitals, while one electron exists in the Cg orbital. In such a 
case, it is necessary to treat both t2g and Cg electrons in the same model Hamiltonian. On 
the other hand, if the low-spin state for n=4 is favored due to the large CEF potential, 
it is enough to consider only the t2g orbital, as in the cases of n=l~3. In fact, for Ad 
electron systems, the effect of CEF potential becomes strong due to the large ion radius. 
For instance, in ruthenates, Ru^^ ion includes four Ad electrons, which form the low-spin 
state with S=l, as deduced from the AF phase of Ca2Ru04. 

For n=5, we find 6-fold degeneracy in the region of small lODq. This is originating 
from the high-spin S=5/2 state with Sz=^5/2, ±3/2 and ±1/2, in which each orbital 
is occupied by one electron. For very large lODq, the ground state degeneracy is still 
six, but the wavefunction is drastically changed in this low-spin S'=l/2 state. It is 
understood that five electrons (or one hole) occupies the t2g orbitals. Namely, there 
are three kinds of S=l/2 state, leading to six- fold degeneracy in total. In principle, the 
intermediate spin state with 5=3/2 is possible, but the region for the intermediate state 
seems to be very limited in the parameter space. 

For n=6, the ground state for small lODq has 15-fold degeneracy. This is the high- 
spin S=2 state, in which two electrons in the Cg sector, while four electrons in the t2g 
section. As shown in the configuration of d^(ii) in Fig. El there are three possibility 
for down spin electron in the tg orbitals. Namely, there are three kinds of 5=2 states, 
leading to 15-fold degeneracy in total. For large lODq, we obtain the singlet ground 
state, which is the low-spin state with S=0, in which t2g orbitals are fully occupied by 
six electrons. 

For n=7, the high-spin state in the region of small lODq has 12-fold degeneracy, by 
accommodating two electrons in the Cg and five electrons in the t2g orbitals. As shown 
in the configuration of (i^(H) in Fig. El there are three possibility for one hole in the tg 
orbitals. Thus, there are three kinds of 5=3/2 states, leading to 12-fold degeneracy in 
total. In the low-spin state, on the other hand, the ground state degeneracy is four. It 
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is easily understood, when we put one more electron in the Cg orbital in addition to the 
d^{L) configuration. 

Both for n=8 and 9, since t2g orbitals are fully occupied and thus, active orbital is 
Cg. For n=8, the ground state with S=l composed of a couple of Cg electrons is triply 
degenerate. For n=9, we find one hole in the Cg orbitals, there are four fold degeneracy in 
the ground state. Note that for n=9, when the tetragonal CEF effect is strong enough, 
as actually found in perovskite copper oxides, only x^ — y"^ orbital becomes active. 

In general, it is necessary to consider all five orbitals in the Hamiltonian for 
the complete analysis of electronic properties of transition metal oxides, even if the 
full Hamiltonian includes several competing tendencies and couplings. However, the 
essential physics can be frequently obtained even using relatively simple models, by 
focusing on the active orbital, emerging out of the competition between the CEF 
potential and the Coulomb interaction. In this sense, we immediately arrive at two 
possibilities: One is t2g model for the cases of d^ ~ d^ and d'^i^L) ~ d^{L). Another is Cg 
model for the cases of (i^(L), d^, and d^. For other cases, it is necessary to consider both 
Cg and t2g orbitals in the same Hamiltonian simultaneously. However, for the cases of 
(i'^(H), it is possible to define a simplified model, which will be explained later for Mn^+ 
ion with the high-spin state. 

For the purpose to express the above models, it is convenient to simplify the 
Coulomb interaction term in a qualitatively correct form, since the representation using 
Racah parameters are exact, but too complicated for further analysis. Under the cubic 
CEF potential, as described in the previous subsection, two Cg and three t2g orbitals 
provide appropriate basis. Then, the Coulomb interaction term should be expressed as 

-Hel-el n / J / ^ ^^^ -'7i'^li72a-2;730-2,740"i"i'yio-i'^i720-2 i730-2'^i74o-n I'^UJ 

i 71727374 o"lo'2 

where the modified Coulomb matrix element I is given by the proper combination of 
the original one f^. Note that the orbital index 7 denotes Cg and t2g orbitals, which are 
defined by linear combination of the original d-electron orbitals labeled by m. 

By using the modified Coulomb matrix element, it is useful to define the so-called 
"Kanamori parameters", U, U', J, and J' |129j . Among them, U is the intra-orbital 
Coulomb interaction, given by 

'-^7 -'7cr,7cr';7cr',7(T; I'^-'-J 

with 0^0'. f/' is the inter-orbital Coulomb interaction, expressed by 

1^1' -'70",7'o'';7'o'',7(T) \y^) 

with 7 7^ 7'. J is the inter-orbital exchange interaction, leading to the Hund's rule 
coupling, written as 

7)7' -'7fT,7'cr';7f7',7'cr7 l*^^/ 

with 7 7^ 7'. Finally, J' is the pair-hopping amplitude between different orbitals, given 
by 

77' ^^(T,'ycr' -,^'(7' ,■^'0-1 \'^^) 



Orbital ordering phenomena in d- and f -electron systems 



7 


i 


f^;,y 


Jry^^l 


xy, yz, zx 


xy, yz, zx 


A-25+C 


3B+C 


x''-y\ 'iz^-r'' 


a;2-2/^ 3^2 _^2 


A-45+C 


AB+C 


xy 


x'-y' 


A+4S+C 


C 


xy 


o 9 9 


A-AB+C 


4B+C 


yz^ zx 


9 9 

X —y 


A-2B+C 


3B+C 


yz, zx 


o 9 9 


A+2B+C 


B+C 



17 



Table 1. Expressions for U' and J by using Racah parameters A, B, and C. Note 
that U^A+4:B+3C for each orbital. For more information, see Refs. |128| and jl29| . 



with 7 7^ 7' and a 7^ a'. 

In Table. 1, we list the values for possible Kanamori parameters. Here we have four 
comments on the relation among Kanamori parameters, (i) A relation J=J' holds, 
which is simply due to the fact that each of the parameters above is given by an 
integral of the Coulomb interaction sandwiched with appropriate orbital wave functions. 
Analyzing the form of those integrals the equality between J and J' can be deduced, (ii) 
For completeness, we explicitly show the orbital indices for the definition of Kanamori 
parameters. In fact, as found in Table. 1, the inter-orbital Coulomb interactions between 
Cg and t2g orbitals depend on the kind of orbitals, while there is no orbital dependence 
in Coulomb interactions among Cg or t2g orbitals. (iii) We point out that among Cg or 
t2g orbitals, another relation U=U'+J+J' holds in any combination of orbitals. This 
relation is needed to recover the rotational invariance in orbital space, (iv) Since the 
largest energy scale among the Kanamori parameters is U , the orbitals are not doubly 
occupied by both up- and down-spin electrons. Thus, only one electron can exist in 
each orbital of the degenerate Cg or t2g sector. Furthermore, in order to take advantage 
of J, the spins of electrons point along the same direction. This is the so-called "Hund's 
rule". 

Then, for the case with active Cg orbital degree of freedom, we can define the 
Coulomb interaction term as 



^e?-el = ^ XI '^i7T'^i7i + ^' Xl "'^"'b 



1.7 



+ JY^ d\^ad\ha'di^<r'diha + J' J^l^aT^ai^^ibiC^ibT + h.c). 



(35) 



where the subscripts, a and b, denote x^ — y'^ and ?>z'^ — r^ orbitals, respectively, rai^o-= 
d\di^a, and ni^=Yl,(j "-170-- ^s easily understood from Table. 1, we find U'=A — AB + C 
and J=AB + C for Cg orbital. 

On the other hand, for the case with active t2g orbital, the Coulomb interaction 
term is expressed by 

y 



i^el-el = UY,^h^^hi + —2^ ^i7^ir 



1.7 



i.77^7' 
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J J' 

+ 2" 5-/ X] (^i-ya^^iYa'di-ya'diY^ + ^ zJ ^i7T'^i7i^i7'i^i7'T' (36) 

i,cr,cr' 7^7' i,7^7' 

where 7 and 7' run in the t2g orbitals, U'=A — 2B + C, and J=3B + C. 

Next we consider the more comphcated case, in which the Coulomb interaction 
term is not expressed only among Cg or t2g orbitals. A typical situation can be found in 
the high-spin state of d'^ configuration, which is relevant to manganites in the cubic 
perovskite structure. Due to the Hund's rule, tetravalent manganese ion includes 
three d electrons to form 5=3/2 in the t2g orbitals. By adding one more electron 
to Mn^+ with three up-spin t2g-electrons, let us consider the configuration for the Mn^"*" 
ion. As mentioned above, there are two possibilities due to the balance between the 
crystalline-field splitting and the Hund's rule coupling, but in 3d electron system such 
as manganites, the high-spin state is realized, since the Hund's rule coupling dominates 
over lODq. 

In order to simplify the model without loss of essential physics, it is reasonable to 
treat the three spin-polarized t2g-electrons as a locahzed "core-spin" expressed by Si 
at site i, since the overlap integral between t2g and oxygen p orbital is small compared 
with that between Cg and p orbitals. Moreover, due to the large value of the total 
spin S=3/2, it is usually approximated by a classical spin. This approximation has 
been tested by using computational techniques |130| 1131] . Thus, the effect of Coulomb 
interaction between Cg and t2g electrons can be effectively included by the strong Hund's 
rule coupling between the Cg-electron spin and localized t2g-spins, given by 

-f^Hund = —Ju 2_^ ^i ■ Sj; (37) 

i 

where Si= '^^a/B 'H-ya'^apdi^fs, Jh(>0) is the Hund's rule coupling between localized t2g- 
spin and mobile Cg-electron, and cr={o'x,o'y,az) are the Pauli matrices. The magnitude 
of Jh is considered to be of the order of J. Here note that Si is normalized as |Si|=l. 
Thus, the direction of the classical t2g-spin at site i is defined as 

Si = (sin6'i cos0i, sin6'isin0i, cos6'i), (38) 

by using the polar angle 6i and the azimuthal angle 0i. 

Note that the effect of the Coulomb interaction is not fully taken into account only 
by -^Hund, since there remains the direct electrostatic repulsion between Cg-electrons. 
Then, by further adding the Cg electron interaction term, we obtain the Coulomb 
interaction term for (i^(H) configuration as 

-f^cl-cl = Huund + -f^el-el- (39) 

Here we use "DE" in the superscript for the Hamiltonian, which is the abbreviation of 
"double exchange". We will explain later the reason why we use "DE", in the section 
explaining the result for manganites. 
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2.3. Electron-phonon interactions 

Another important ingredient in the model Hamiltonian for transition metal oxides is 
the lattice distortion coupled to Cg or t2g electrons. In particular, the orbital degeneracy 
is lifted by the Jahn- Teller distortion. It is difficult to explain all possible types of 
electron-phonon coupling term, but here we explain the coupling between degenerate Cg 
orbital and the distortion of the octahedron, since later we will focus on the magnetic 
properties of manganites in the perovskite structure. The coupling of the same type of 
distortions of octahedron to t2g orbitals is also discussed in this subsection. 

The basic formalism for the study of electrons coupled to Jahn- Teller modes has 
been set up by Kanamori |132j . He focused on cases where the electronic orbitals are 
degenerate in the undistorted crystal structure, as in the case of manganese ion in an 
octahedron of oxygens. As explained by Kanamori, the Jahn- Teller effect in this context 
can be simply stated as follows. When a given electronic level of a cluster is degenerate 
in a structure of high symmetry, this structure is generally unstable, and the cluster will 
present a distortion toward a lower symmetry ionic arrangement \\\V^\ . In the case of 
Mn^+, which is orbitally doubly degenerate when the crystal is undistorted, a splitting 
will occur when the crystal is distorted. The distortion of the MnOe octahedron is 
"cooperative" since once it occurs in a particular octahedron, it will affect the neighbors. 
The basic Hamiltonian to describe the interaction between electrons and Jahn- Teller 
modes was written by Kanamori. In the adiabatic approximation, it is given in the 
form of 

i^jT = g $^(Q2ir.i + Q3iT-.i) + fcjT ^(Q^i + Q^i)/2, (40) 

i i 

where g is the coupling constant between the Cg-electrons and distortions of the MnOg 
octahedron, Q2i and Qsi are normal modes of vibration of the oxygen octahedron that 
remove the degeneracy between the electronic levels, and fcjT is the spring constant for 
the Jahn- Teller mode distortions. We note that the sign of g depends on the kind of ion 
and electron configuration. The pseudospin operators are defined as 

Txi = ^{dladiha + dlji^), r,i = ^{dlji^ - rfL^iba)- (41) 

a (T 

In the expression of ifjx, a Tyi-term does not appear for symmetry reasons, since 
it belongs to the A2u representation. The non-zero terms should correspond to the 
irreducible symmetric representations of EgXEg, namely. Eg and Aig. The former 
representation is expressed by using the pseudo spin operators r^i and r^i as discussed 
here, while the latter, corresponding to the breathing mode, is discussed later. 
Following Kanamori, (^21 and Q^i are explicitly given by 

Q2i = {Xn - Xsi - Y2i + Y^i)/V2, (42) 

and 

Q3i = (2Z5i - 2^6; - Xn + Xsi - F2i + ni)/v^, (43) 
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where X^i, Y^i, and Z^i are the displacement of oxygen ions from the equihbrium 
positions along the x-, y-, and ;z-direction, respectively. Explicitly, they are included in 
the positions of oxygen ions as Rii=(a/2+Xii, 0, 0), R2i=(0, —0/2+^21, 0), R3i=(— a/2+ 
Xgi, 0, 0), R4i=(0, a/2 + F4i, 0), R5i=(0, 0, a/2 + Zg;), and R6i=(0, 0, -a/2 + Zgi). Here 
we consider the deformation of octahedron from the cubic case. The convention for the 
labeling /i of coordinates is shown in Fig. ^ 

In order to solve the local Jahn- Teller problem, it is convenient to scale the phononic 
degrees of freedom as 

Q2i = i9/kjT)q2h Qsi = (5'/^JT)g3i, (44) 

where g/kj^ is the typical length scale for the Jahn- Teller distortion, which is of the 
order of 0.1 A, namely, a few percent of the lattice constant. When the Jahn- Teller 
distortion is expressed in the polar coordinate as 

gsi = Qi sin ^i, ggj = g; cos ^i, (45) 

the ground state is easily obtained as (— sin[^i/2]cil^^ -|- cos[^i/2]c/i^^)|0) with the use of 
the phase ^i. The corresponding eigenenergy is given by — -Ejt, where Ejt is the static 
Jahn- Teller energy, defined by 

EjT = //(2fcjT). (46) 

Note here that the ground state energy is independent of the phase ^i. Namely, the shape 
of the deformed isolated octahedron is not uniquely determined in this discussion. In the 
Jahn- Teller crystal, the kinetic motion of Cg electrons, as well as the cooperative effect 
between adjacent distortions, play a crucial role in lifting the degeneracy and fixing the 
shape of the local distortion. 

To complete the electron-phonon coupling term, it is necessary to consider the 
breathing mode distortion, coupled to the local electron density as 

i i 

where ni=^ ^ dj dj^a and the breathing-mode distortion Qn is given by 

Qii = (Xn - Xsi + r2i - Kii + ^5i - ^6i)/v^, (48) 

and fcbr is the associated spring constant. Note that, in principle, the coupling constants 
of the Cg electrons with the Qi, Q2, and Q3 modes could be different from one another. 
For simplicity, here it is assumed that those coupling constants take the same value. 
On the other hand, for the spring constants, a different notation for the breathing 
mode is introduced, since the frequency for the breathing mode distortion has been 
found experimentally to be different from that for the Jahn- Teller mode. This point 
will be briefly discussed later. Note also that the Jahn- Teller and breathing modes are 
competing with each other. As it was shown above, the energy gain due to the Jahn- 
Teller distortion is maximized when one electron exists per site. On the other hand, the 
breathing mode distortion energy is proportional to the total number of Cg electrons per 
site, since this distortion gives rise to an effective on-site attraction between electrons. 
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By combining the Jahn- Teller mode and breathing mode distortions, the electron- 
phonon term for the Cg model is summarized as 

i^e'lV = ^JT + ^br. (49) 

This expression depends on the parameter P=khT/kjT, which regulates which distortion, 
the Jahn- Teller or breathing mode, plays a more important role. This point will be 
discussed in a separate subsection. 

When we simply ignore buckling and rotational modes and consider only the Jahn- 
Teller-type distortions of octahedron, the electron-lattice coupling for t2g orbitals is 
given by |134j 

^cl-ph = 9 $^(Q.i^ixy + Qxi^iy. + QyiUi,.) + (fcjT/2) J2{Qli + Qsi), (50) 
i i 

where Q^i={-l/2)Q3i+{V3/2)Q2u Qj;i=(-l/2)g3i-(V3/2)Q2i, and Q,i=Q3i- Later 
we will consider the effect of this type of electron-phonon coupling on the magnetic 
structure of t2g electron systems. 

Let us now consider the cooperative effect. Although we have not explicitly 
mentioned thus far, the distortions at each site are not independent, since all oxygens are 
shared by neighboring MnOg octahedra, as easily understood by the explicit expressions 
of Qii, Q2i, and Qsi presented before. A direct and simple way to consider such 
cooperative effect is to determine the oxygen positions Xii, X^i, 1^21, ^5i, -^si, and Zgi, 
by using computational method such as the Monte Carlo simulations and numerical 
relaxation techniques. To reduce the burden on the numerical calculations, for instance, 
the displacements of oxygen ions are assumed to be along the bond direction between 
nearest neighboring manganese ions. In other words, the displacement of the oxygen ion 
perpendicular to the Mn-Mn bond, i.e., the buckling mode, is usually ignored. As shown 
later, even in this simplified treatment, several interesting results have been obtained 
for the spin, charge, and orbital ordering in manganites. 

Rewriting Eqs. (j^2|) . P3|l . and P8|) in terms of the displacement of oxygens from 
the equilibrium positions, we express the distortions as 

Qii = (A,i + Ayi + A,i)/v^, 

Q2i = (Axi - Ayi)/v^, (51) 

g3i = (2A,i-Axi-Ayi)/y6, 

where Aai is given by 

Aai = ut- <_a, (52) 

with uf being the displacement of oxygen ion at site i from the equilibrium position 
along the a-axis. In the cooperative treatment, the {m}'s are directly optimized in the 
numerical calculations P3ll26j . On the other hand, in the non-cooperative calculations, 
{Q}'s are treated instead of the {u}'s. In the simulations, variables are taken as {Q}'s 
or {"uj's, depending on the treatments of lattice distortion. 

Finally, we briefly comment on the effect of macroscopic distortion. In the above 
treatment, we assume the cubic symmetry for the equilibrium positions of oxygens. 
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Table 2. Expressions for hopping amplitude between d orbitals. As for details, see 
Ref. |135| . We use a direction cosine as {l,m,n) for the direction of hopping from 7 
to 7' orbitals. We show the contributions from the hoppings through a- and 7r-bonds. 



but in actuality, the crystal structure is frequently deviated from the cubic symmetry. 
Although we cannot determine the stable crystal structure in the present treatment, the 
effect of macroscopic distortions is included as offset values for the distortions, which 
are given by 



Ql 



(0) 



(5Lx + 5Ly + 5L,)/v^, 
(5Lx - 5Ly)/v^, 
Qf = {25L, - 6L^ - SLy)/V6, 



Qf 



(53) 



where 6La=La — L, the non-distorted lattice constants are La, and L=(Lx + Ly + Lz)/3. 
Note that L^ is determined from the experimental results. By adding Qjj, in the right- 
hand side of Q^i in Eq. (jHH), it is possible to consider the effect of the deviation from 
the cubic symmetry 



2.4- Electron hopping 

In previous subsections, we have discussed the local electron state due to CEF potential 
and Coulomb interaction. We have also considered an additional electron-phonon 
coupling term. Since the possible local terms have been completed, let us consider the 
intersite effect in the next step. In actual materials, there are several kinds of intersite 
effects. Among them, we consider the electron hopping between adjacent d electron 
orbitals. For transition metal oxides, such a hopping process occurs through the oxygen 
ion, but in the formalism, it is enough to consider the direct hopping (i-electron orbitals. 
Effect of oxygen will be discussed later. 

Fortunately, the hopping amplitudes have been already evaluated and tabulated for 
any combination of electron orbitals. We can simply consult with the table of Slater- 
Koster integral |135j . which is the general scheme for the overlap integral between 
adjacent electron orbitals. The kinetic term for Cg or t2g electrons is expressed in a 
common form as 



^Mn = EE^77'<.^: 



i+a7'(T) 



(54) 



i,a,(T 7,7' 
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where F is the irreducible representation, a is the vector connecting nearest-neighbor 
sites, and t^ , is the nearest-neighbor hopping amphtude between 7- and 7'-orbitals 
along the a-direction. 

For the cubic lattice composed of transition metal ion, we consider the three axis 
directions, x=(l,0,0), y=(0,l,0), and z=(0,0,l). Then, t^y for Cg orbital is given in a 
2x2 matrix form as 



t^ = ti 



for x-direction, 



t^ = ti 



for y-direction, and 



t^ = t, 



3/4 -73/4 
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for z-direction, where ti = {dda) and (dda) is one of Slater integrals |135j . It should be 
noted that the signs in the hopping amplitudes between different orbitals are different 
between x- and y-directions. 

On the other hand, for t2g orbitals, we obtain the hopping amplitudes in a 3x3 
matrix form as 
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(59) 



(60) 



for 2;-direction, where t2={dd7i). 

Now let us consider explicitly the effect of oxygen orbitals. Since oxygen ion 
is placed in the middle of transition metal ions in the cubic perovskite, the main 
hopping process should occur via oxygen 2p orbitals. Thus, the d-electron hopping 
can be expressed by (pda) or {pdn), which is the overlap integral between d and p 
orbitals, divided by the energy difference between d and p orbitals. It is possible 
to calculate the (i-electron hopping via oxygen 2p orbitals, by consulting again the 
Slater-Koster table for the overlap integral between d and p orbitals. However, due to 
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the symmetry argument, we easily understand that the form of hopping amphtude is 
invariant, after the redefinition of ti=—{dpa)'^/{ep — Ed) and t2=—{dpT^Y/{ep — Ed) for 
Cg and t2g electron cases, respectively, where Ed and Ep are the energy levels for d and p 
electrons, respectively. 

We have considered the nearest neighbor hopping, but in actuality, it is necessary 
to consider higher neighbors in order to reproduce the Fermi surface observed in the 
experiments such as de Haas-van Alphen measurements. However, there is no essential 
difficulty for the consideration of higher neighbors by consulting the Slater-Koster table. 

2.5. Intersite interaction term 

In the previous subsections, we have considered the local term and kinetic motion of 
d electrons. Of course, due to the combination of these terms, intersite interaction 
terms effectively appear. In particular, in the strong-coupling limit, orbital-dependent 
superexchange terms can be obtained, leading to the complex magnetic structure with 
orbital ordering. Such effects are considered automatically, as long as we consider the 
problem within the electronic model. 

However, in the model for the high-spin state of (i^ electron configuration, it is 
necessary to explicitly add an extra term between localized t2g spin. As explained above, 
due to the Hund's rule coupling between Cg and t2g interaction, we can easily understand 
that Cg electrons can move smoothly without any energy loss, if spins of Cg and t2g are 
in the same direction. Namely, in order to gain the kinetic energy of Cg, t2g spins array 
ferromagnetically. This is a simple explanation for the appearance of ferromagnetism 
in the orbital degenerate system, in particular, in manganites. However, there should 
exist AF intersite coupling between neighboring t2g spins due to the superexchange 
interaction, given by 

-f^inter-sitc = -^AF 2_^ ^i ■ Sj, (61) 

(iJ> 

where (i, j) denotes the pair of nearest neighbor sites and Jaf is the AF coupling between 
neighboring t2g spins. This term should be added to the model for manganites. As we 
will explain later, the competition between FM tendency to gain kinetic energy and AF 
energy to gain magnetic energy is a possible source of complex magnetic structure in 
manganites. 

In addition to the effective coupling among localized spins, sometimes we consider 
another intersite effect originating from long-range Coulomb interaction, even if it is 
screened in actual materials. In order to include such effect in an effective manner, we 
also add 

^Stcr-sitc = V^E^'^J' (62) 

(iJ> 

where V denotes the nearest-neighbor Coulomb interaction. Since V has a tendency 
to stabilize the charge ordering, there occurs competition between striped spin order 
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and bipartite charge ordering. This is another source of complex spin-charge-orbital 
ordering. 

2.6. Summary 

We have completed the preparation of all components for the model Hamiltonian. As a 
short summary of this section, we show three types of model Hamiltonians due to the 
appropriate combination of several terms. 

For the system with active Cg or t2g orbital degree of freedom, we can consider the 
orbital degenerate model, expressed in a common form as 

Hr = -n]^in + -ncl_cl + -"cl-ph + -"intcr-site) (63) 

where V denotes Cg or t2g- Note that the inter-site Coulomb interaction term is explicitly 
added here, but depending on the nature of the problem, this term may be ignored. 

For the case of (i^(H) in which both Cg and t2g orbitals are included, we can define 
the following model: 

-f^DE = -f^kfn + -f^cl-cl + -^cl-ph + -^inter-sitc + -^intcr-sitc- (64) 

This expression is believed to define an appropriate starting model for manganites, but 
unfortunately, it is quite difficult to solve such a Hamiltonian. In order to investigate 
further the properties of manganites, further simplifications are needed. This point will 
be discussed in detail in the next section. 

3. Orbital physics in manganites 

In the complicated phase diagram for manganites, there appear so many magnetic 
phases. A key issue to understand such richness is the competition between itinerant 
and localized tendencies contained in manganites. As mentioned in the model 
construction, Cg electrons can gain the kinetic energy when the background t2g spins 
array ferromagnetically, leading to a metallic FM phase in manganites. On the other 
hand, in order to gain magnetic energy between localized t2g spins, there occurs AF 
phase with insulating tendency. In one word, the competition between FM metallic and 
AF insulating phases is the origin of complex phase diagram of manganites. 

As we will review very briefiy in the next subsection, the metallic tendency has been 
discussed in the concept of double exchange for a long time, and the essential point has 
been considered to be well understood. However, the tendency toward insulating phase 
has not been satisfactorily understood, mainly due to the complexity of multi-degrees 
of freedom such as spin, charge, and orbital. In particular, "orbital ordering" is the 
remarkable feature, characteristic to manganites with active Cg orbital. In this section, 
spin, charge, and orbital structure for the typical hole doping in the phase diagram of 
manganites is focused by stressing the importance of orbital ordering. 
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Figure 4. (a) Sketch of the double exchange mechanism which involves two Mn ions 
and one O ion. (b) The intermediate state of the process (a), (c) The mobility of 
Bg-electrons improves if the localized spins are polarized. 



3.1. Concept of double-exchange 

Since the historical review of theoretical and experimental works on manganites 
have been found in some articles and textbooks, we simply refer literatures such as 
Refs. jSl in Ej • However, it is instructive to mention here the meaning of "double 
exchange" (DE), which is important basic concept for manganites, by following the 
previous review article [H]. 

In the earth stage of the research on manganites, it was the task to clarify the 
qualitative aspects of the experimentally discovered relation between transport and 
magnetic properties, namely the increase in conductivity upon the polarization of the 
spins. The concept of "double exchange" was proposed by Zener as a way to allow 
for charge to move in manganites by the generation of a spin polarized state. The 
DE process has been historically explained in two somewhat different ways. Originally, 
Zener considered the explicit movement of electrons, as shown in Fig.|^a). This process 
is schematically written as Mn'^J'02|,3|Mn''"'" -^ Mn'^+Oi|^3j^Mn2t |136j . where 1, 2, and 
3 label electrons that belong either to the oxygen between manganese, or to the Cg-level 
of the Mn-ions. In this process, there are two simultaneous motions involving electron 
2 moving from the oxygen to the right Mn-ion, and electron 1 from the left Mn-ion to 
the oxygen. This is the origin of the name of "double exchange" . 

The second way to visualize DE processes was presented in detail by Anderson and 
Hasegawa |137j . and it involves a second-order process in which the two states described 
above go from one to the other using an intermediate state Mn^tOai Mng^", as shown 
in Fig. mb). In this context, the effective hopping for the electron to move from one 
Mn-site to the next is proportional to the square of the hopping involving the p-oxygen 
and d-manganese orbitals (tpd). Following Anderson and Hasegawa, let us consider a 
two-site problem, in which one itinerant electron with hopping amplitude t between 
sites 1 and 2 is coupled with localized spin S at each site. The coupling is assumed 
to be ferromagnetic and the magnitude is defined as J. For t=0, the local ground 
state is labelled by S'-|-l/2 with the energy of —JS per site, while the excited state is 
specified as 5*— 1/2 with the energy of J{S + 1) per site. For ty^O with t<^J, as observed 
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in manganese ions, the effective hopping aniphtude t^s is in proportion to the overlap 
integral between 15+1/2, S, Sq) and \S, 5 + 1/2, Sq), where l^i, 5*2, Sq) indicates the state 
with spin Si at site i and total spin Sq. Note that 5*0 is composed of two localized spins 
and one itinerant spin. The overlap integral is evaluated as (5+1/2, S, So\S, 5 + 1/2, Sq) 
=2(5 + 1)W{S, 1/2, So, S;S + 1/2, 5 + 1/2), where W is the so-called Racah coefficient 
for the combination of three spins. By using some relations for the Racah coefficient 
[I3H1, we obtain \t^s/t\ = {So + l/2)/(25 + 1). 

For 5=1/2, we can intuitively understand the meaning of the reduction factor. 
When the two localized spins are parallel, we obtain 5o=3/2 due to the large J and the 
reduction factor is unity. This is understood by the fact that the local triplet formed 
by the large J is not destroyed in the course of electron hopping motion. However, 
when the two localized spins are anti-parallel, 5o=l/2 and the reduction factor is 1/2. 
In this case, it is necessary to reconstruct the local triplet state after the hopping of 
electron, leading to the effective reduction of the electron hopping. For the case of 
large 5, the localized spins are considered classical. When we define an angle 6 between 
nearest-neighbor ones, the above overlap integral is easily evaluated by rotating the 
itinerant electron basis so as to be parallel to the localized spin direction. Then, we 
obtain tefj=tcos(^/2), as shown by Anderson and Hasegawa. If 6=0 the hopping is the 
largest, while if 6=7t, corresponding to an AF background, then the hopping cancels. 

Note that the oxygen linking the Mn-ions is crucial to understand the origin of 
the word "double" in this process. Nevertheless, the majority of the theoretical work 
carried out in the context of manganites simply forgets the presence of the oxygen and 
uses a manganese-only Hamiltonian. It is interesting to observe that FM states appear 
in this context even without the oxygen. It is clear that the electrons simply need a 
polarized background to improve their kinetic energy, as shown in Fig. ^c), in similar 
ways as the Nagaoka phase is generated in the one-band Hubbard model at large U/t. 
This tendency to optimize the kinetic energy is at work in a variety of models and 
the term double-exchange appears unnecessary. However, in spite of this fact, it has 
become customary to refer to virtually any FM phase found in manganese models as 
"DE induced" or "DE generated" , forgetting the historical origin of the term. In this 
review, a similar convention will be followed, namely the credit for the appearance of 
FM phases will be given to the DE mechanism, although a more general and simple 
kinetic-energy optimization is certainly at work. This is also the reason why we have 
used the abbreviation "DE" in the model for manganites. 

Early theoretical work on manganites carried out by Goodenough |139j explained 
many of the features observed in the neutron scattering experiments on Lai^-xCaxMnOs 
by Wollan and Koehler [8], notably the appearance of the A- type AF phase at x=0 and 
the CE-type phase at x=0.5. The approach of Goodenough was based on the notion of 
"semicovalent exchange". Analyzing the various possibilities for the orbital directions 
and generalizing to the case where Mn^"*" ions are also present, Goodenough arrived to 
the A- and C15-type phases of manganites very early in the theoretical study of these 
compounds. In this line of reasoning, note that the Coulomb interactions are important 
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to generate Hund-like rules and the oxygen is also important to produce the covalent 
bonds. The lattice distortions are also quite relevant in deciding which of the many 
possible states minimizes the energy. However, it is interesting to observe that in more 
recent theoretical work described below in this review, both the A- and CE-type phases 
can be generated without the explicit appearance of oxygens in the models and also 
without including long-range Coulomb terms. 

3.2. Topological aspects of orbital ordering 

In Sec. 2, we have already set the model Hamiltonian for manganites -^de- Before 
proceeding to the exhibition of the theoretical results on this model, we explain the 
essential point of manganites from a purely theoretical viewpoint. We believe that 
it is an important issue to establish a simple principle penetrating the complicated 
phenomena. 

3.2.1. A simplified model In order to extract the essential feature of manganites, let 
us define a minimal model, since i^oE is still a complex model. We note that there 
exist two important ingredients which should be kept even in the minimal model. One 
is the existence of orbital degree of freedom and another is a competition between FM 
metallic and AF insulating tendencies. In order to minimize the model by keeping these 
two issues, first we simply ignore the interaction terms, -ffei-ei, -f^ei-ph' ^^"^ -^^Ster-site- 
Second, we take an infinite limit of the Hund's rule coupling, Jh, between Cg electron 
and t2g spins. Then, the direction of Cg-electron spin perfectly follows that of t2g spin. 
We can suppress the spin index, if we define the spinless operator at each site in which 
the spin direction is fixed as that of t2g spin at each site. Namely, the model is virtually 
expressed by using spinless operators with orbital degree of freedom. Then, we obtain 
a simplified double-exchange model as 

H = -Y. A,i+at^,yc|^Ci+ay + Jay J2 ^i^i^ (65) 

where Cj^ is the annihilation operator for spinless d electron in the 7 orbital at site 
i, the hopping amplitudes are given in Eqs. (j3^ . (j^ . and (jHZI), where ti=— to with 
to={pdaY / {Sp — Ed), and Dij is the so-called double-exchange factor, given by 

A,j = cos(^i/2) cos(6'j/2) + sin(^i/2) sin(6'j/2)e^^('^--'^j). (66) 

Here 9i and 0i denote the polar and azimuthal angles of t2g spin at site i, respectively. 
This factor expresses the change of hopping amplitude due to the difference in angles 
between t2g-spins at sites i and j. Note that the effective hopping in this case is a 
complex number (Berry phase), contrary to the real number widely used in a large 
number of previous investigations. As for the effect of the Berry phase in the case of 
the one-orbital DE model, readers should refer Ref. |14nj . 

Furthermore, when we assume the Ising t2g spins, the double-exchange factor L>i j 
denotes or 1 depending on the spin configuration. Namely, Dij=l for FM t2g spin 
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Figure 5. Ground-state energy vs. Jaf for (a) x=0.0 and (b) x=0.5 of a simplified 
model. The unit cell is taken as a 4 x 4 lattice, (c) Schematic views for spin structure 
of FM, E-type AF, CE-type AF, AF2, and G-AF phases. 



pair at sites i and j, while Di^=0 for AF pair. One may feel that the model seems to be 
oversimplified, but as will see later, we can grasp an essential point of complex phases 
of manganites, since the competition between FM metallic and AF insulating natures 
is correctly included in this model. 



3.2.2. Band-insulating state First let us consider two limiting situations, Jaf=0 and 
Jaf ^ if)- For simplicity, a two-dimensional (2D) square lattice is taken here. For 
Jaf=0, it is easy to understand the appearance of FM metallic phase, since the 
Hamiltonian includes only the kinetic term of electrons in this limit. On the other 
hand, for Jaf ^ ^O) AF insulating phase should appear due to the magnetic energy gain. 
Then, what happens for intermediate value of Jaf? Naively thinking, it is possible to 
consider the mixture of FM and AF regions, in order to gain both kinetic and magnetic 
energies. For instance, we can consider the C-type AF phase, in which one-dimensional 
(ID) FM chains are antiferromagnetically coupled with each other. However, there is 
no special reason to fix the shape of the FM region as straight ID chain. It may be 
possible to have zigzag shape for the FM region. 

In order to determine the optimal shape of the FM region, we perform simple 
simulations for magnetic structure in the 2D lattice |141j . Since the model does not 
include explicit many-body interaction among Cg electrons, we can solve the problem on 
the periodic lattice composed of an appropriate unit cell such as 4 x 4 cluster. Namely, 
we prepare all possible patterns for t2g spin configuration and evaluate the ground 
state energy on each magnetic structure by changing the value of Jaf- Note that the 
calculations have been done on the momentum space by introducing the Bloch phase 
factor at the boundary. 
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The results for x=0.0 and 0.5 are shown in Figs.Efa) andO^b), where x denotes hole 
doping from the case of quarter filling, i.e., one electron per site with two orbitals. Note 
that x=0.0 corresponds to the case in which all sites are occupied by trivalent manganese 
ions, while x=0.5 indicates the situation in which half of the sites are occupied by Mn^"*". 
For x=0.0, there are four regions in Fig. Efa). As mentioned above, we obtain 2D FM 
phase at Jaf=0, while for large Jaf, the G-AF phase appears, as expected. In the 
intermediate region, we observe the striped spin phase, characterized by zigzag FM 
path, as shown in Fig. Efc). This structure is just the E-type AF phase. Note that 
the magnetic energy is cancelled in the E-type phase, since the numbers of FM and 
AF bonds are exactly equal. We also find a narrow window for another AF phase, 
called "AF2", between the E-type and G-type AF phases. In the present context, the 
appearance of this phase is not important. 

For x=0.5, as shown in Fig.El^b), the situation looks similar with the case of x=0.0, 
but among four magnetic phases, the spin structure in the region labeled by "CE" differs 
from that of the E-type phase at x=0.0. Namely, the period of the zigzag is different 
from that for x=0.0, as shown in Fig. El^c). We emphasize here that the spin structure 
in the intermediate coupling is composed of a bundle of spin FM chains, each with the 
zigzag geometry, and with AF interchain coupling. This is just the CI5-type AF phase. 
Note that the magnetic energy is also cancelled, since the numbers of FM and AF bonds 
are exactly equal in the CE-type phase. 

Now we consider the reason why such complicated structure appears. For the time 
being, let us discuss what happens if the zigzag geometry of CEl- or E-type is assumed, 
and how it compares with a straight line. A straightforward way is to calculate the 
energy band for the Cg electron system on the zigzag ID path, since Cg electrons can 
move only in the FM region in the simplified model due to the double exchange factor. 
First we consider the C-type AF phase characterized by straight ID path, even though 
it is not the ground state. As shown in Sec. 2, the hopping amphtudes of Cg electrons 
depend on the direction, but as easily checked by a diagonalization of 2 x 2 matrix, 
due to the cubic symmetry, the energy band does not depend on the chain direction. 
Then, by choosing the hopping direction along the ^-axis, we easily obtain Ek=0 and 
— 2to cos k, since there is non-zero hopping amplitude only between 3z'^ — r^ orbitals 
along this direction. 

To solve the present one-body problem in the zigzag ID case, unit cells are defined 
as shown in Figs.in^a) andlHfb), in which the hopping amplitudes change with a period of 
two or four lattice spacings for E- and CE-type structure, respectively, since the hopping 
direction changes as {■ ■ ■ ,x,y,x,y, ■ ■ ■} and {■ ■ ■ ,x,x,y,y,- ■ ■} along the zigzag chain, 
with t^^=—tj^j^ ioi ji ^ u according to the values of the hopping amplitudes discussed 
before. This difference in sign, i.e., the phase change, is essential for this problem. To 
make this point clear, it is useful to transform the spinless Cg-electron operators by using 
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Figure 6. Definition of the unit cell for (a) E-type and (b) CE-type ID zigzag chain. 
Momentum is defined along the zigzag chain. Energy band structure for (c) E-type 
and (d) CE-type structure. Solid and broken curves denote occupied and unoccupied 
bands, respectively. Note that the flat band in (d) has four-fold degeneracy. 



a unitary matrix as 





After simple algebra, iifkin is rewritten as 



'kin 



to/2 5^(alai+a + A'A+a + e^'^^aJA+a + e-'f-f3Jai+^j, 



3«</'a^,t/ 



-«0a,Qt, 



(67) 



(6^ 



where the phase (pa depends only on the hopping direction, and it is given by 0x=— 0, 
4>y=(p, and (pz=7!', with 0=7r/3. Note that the Cg-electron picks up a phase change when 
it moves between different neighboring orbitals. In this expression, the effect of the 
change of the local phase is correctly included in the Hamiltonian. 

Here we solve the problem in the E-type structure. Details of the solutions of the 
CE-type structure have found in the previous review [^. To introduce the momentum 
k along the zigzag chain, the Bloch's phase e^"^^ is added to the hopping term between 
adjacent sites. Then, the problem is reduced to finding the eigenvalues of a 4 x 4 matrix, 
given by 



h 



(^aA 




/^.l\ 




= ek 




\ ^/32 / 




V ^/32 / 



(69) 



where ipaj and ipjjj are the basis function for a- and /^-electrons at the j-site of the unit 
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cell, respectively, and the Hainiltonian matrix h is given by 



h = -to 



/ 





cosk 


cos(A; — 








cos{k + 0) 


cos A; 


cosfc 


cos(fc + (p) 








I cosfA; — 6' 


cosk 









(70) 



/ 



To make the Hamiltonian in a block diagonalized form, we introduce two kinds 
of bonding and antibonding states as $f = (^qi ± ^/3i ± ^a2 + ^/32)/2 and $^= 
ii^ai T V'/3i i '^«2 — "0/32 )/2. For the states 1 and 2, we obtain two eigen equations 
as 



and 



ito/2) 



ito/2) 



3 cos k ^/S sin k 
v3sinA; — cosfc 





cosfc 
-yosin A; 



-v3sin A; 
-3 cos k 



$0 



E 



(2) 



*2 



respectively. We can easily diagonalize each 2x2 matrix and obtain 



El;' = (to/2) (- COS k ± Vcos2 A; + 3) 



and 



E 



(2) 



(71) 



(72) 



(73) 



(74) 



,^ _ (to/2) (cos k ± Vcos2 A; + 3) . 

For the case of CE-type AF phase, it is necessary to solve eigen value problem in 
the 8x8 matrix. Readers interested in the detail of the calculations can refer the review 
article by Dagotto, Hotta, and Moreo [HJ. Here we show only the results. Eight eigen 
energies have been obtained as 

Ek = 0, ±to^/2 + cos{2k), ±to^2-cos{2k), (75) 

where the flat band 6^=0 has four-fold degeneracy. 

The band structures for E-type and CE-type zigzag path are shown in Figs. IHl^c) 
andin^d). Note that the solid curves denote the occupied bands. The most remarkable 
feature is that the system is band-insulating, with a bandgap of value to for x=0.0 PUI I^ 
and x=0.5 p^41j . Remark that the band-insulator state at x=0.5 was independently 
obtained in Refs. J1421 11431 I144J . This band insulating state, without any explicit 
potential among the electrons moving along the zigzag chains, is caused by the phase 
difference between t^^, and t^^. Since to is at least of the order of lOOOK, this band- 
insulating state is considered to be very robust. 

Intuitively, the band-insulating state of the zigzag AF structure originates in the 
presence of a standing-wave state due to the interference between two traveling waves 
running along the x- and y-directions. In this interference picture, the nodes of the 
wavefunction can exist on the "corner" of the zigzag structure, and the probability 
amplitude becomes larger in the "straight" segment of the path. Thus, even a weak 
potential can produce the charge and orbital ordering based on this band-insulating 
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phase. In fact, if some potential is included into such an insulating phase, the system 
maintains its "insulating" properties, and a spatial modulation in the charge density 
appears. For x=0.5, since the charge density should be increased in the sites 2 and 4 
in Fig. ini^b), it is easy to understand that checker-board type charge ordering occurs, 
when some potential is further included in the model. 

Concerning the orbital shape, we should point out that the Cg-electron orbital is 
maximally polarized along the transfer direction in order to gain the kinetic energy. This 
effect may be called an "orbital double-exchange" in the sense that the orbitals align 
along the axis direction to make the transfer of the electron smooth, similarly as the FM 
alignment of t2g spins in the usual DE mechanism. Namely, on the straight-line part in 
the X- and y-direction, the orbital is polarized as 3x^ — r^ and 3y^ — r^, respectively. 

On the other hand, for the case of x=0, there is no straight-line part in the E-type 
zigzag structure. In this case, rather the cooperative Jahn- Teller effect is essential to 
determined the orbital shape, since each site is Jahn- Teller active. We cannot determine 
the orbital ordering pattern within the present simple discussion. It is necessary to 
consider a more realistic Hamiltonian. The actual orbital ordering will be discussed 
later. 

3.2.3. Topological number In the previous subsection, we have emphasized that the 
shape of zigzag path plays an important role for the determination of the CE- and 
E-type AF phases. Now let us consider the quantity to specify the zigzag shape. For 
the purpose, we include the coupling between Cg electrons and JT phonons. Effect of 
Coulomb interaction will be discussed later. The model is given by 

H = — 2_^ -Di,i+ai^^'Ci^Ci+a7' + J A¥ ^ Si ■ Sj 

+ EjT Xl[2(?2irrf + gsiT.i) + ql + g|i] , (76) 

i 

where Ejt is the static Jahn-Teller energy, r^i= c\^Cib+c\^^Cia, and T^i= c[^Cia-c|^Cib. By 
using the phase ^i defined in Eq. (I45j) . which is the angle to specify the Jahn-Teller 
distortion, it is convenient to transform Cia and Qb into the "phase-dressed" operators, 
Cia and Cib, as 

^ ~^^A= m) ('"']■ (77) 

where the unitary matrix -R(^i) is given by 

^«''-^ l,-smK,/2] cosB,/2lj- <"**' 

Note that if ,^i is increased by 2tt, the SU(2) matrix itself changes its sign. This is the 
same phenomenon found in spin wavefunction, since in general, spinor is isomorphic to 
the wavefunction of a two-level system. In 1950's, Longuet-Higgins et al. have pointed 
out that the electron wavefunction of Jahn-Teller molecule changes its sign for the 
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27r-rotation in the parameter space in the adiabatic approximation |145j . Note that the 
total wavefunction, given by the product of electron and phonon wavefunctions, is always 
single- valued, since the phonon part also changes its sign for the 27r-rotation. The spinor- 
like wavefunction for the electron part appears due to the adiabatic approximation for 
the JT system. It has been also mentioned that the change of sign is regarded as the 
effect of the Berry phase |146j . 

In order to keep the transformation unchanged upon a 27r-rotation in ^i, a phase 
factor e*^'/^ is needed. This is also regarded as the effect from the phonon wavefunction. 
In the expression for the ground state of the single JT molecule, namely, the single-site 
problem discussed before, this phase factor has not been considered explicitly, since the 
electron does not hop around from site to site and the phases do not correlate with each 
other. It was enough to pay attention to the fact that the electron wavefunction at a 
single site is double-valued. However, in the JT crystal in which Cg electrons move in the 
periodic array of the JT centers, the addition of this phase factor is useful to take into 
account the effect of the Berry phase arising from the circular motion of Cg-electrons 
around the JT center J147[ll48j . It could be possible to carry out the calculation without 
including explicitly this phase factor, but in that case, it is necessary to pay due attention 
to the inclusion of the effect of the Berry phase. The qualitative importance of this effect 
will be explained later. 

Note also that the phase ^i determines the electron orbital set at each site. In the 
previous section, the single-site problem was discussed and the ground-state at site i 
was found to be 

|"b") = [-sin(6/2)rfL + cos(6/2)4j|0), (79) 

which is referred to as the "b" -orbital, namely the combination with the lowest-energy 
at a given site. The excited-state or "a"-orbital is simply obtained by requesting it to 
be orthogonal to "b" as 

I "a") = [cos(6/2)rfL + sin(6/2)rfL]|0)- (80) 

For instance, at ^i=27r/3, "a" and "b" denote the dy2_^2- and d3^2_r2-orbitals, 
respectively. It should be noted here that (i3^2_j,2 and d^y2_^2 never appear as the local 
orbital set. Sometimes those were treated as an orthogonal orbital set to reproduce 
the experimental results, but such a treatment is an approximation, since the orbital 
ordering is not due to the simple alternation of two arbitrary kinds of orbitals. 

Using the above described transformations, the model Eq. (fTUI) is rewritten after 
some algebra as 

H = -^ ^i,i+a^^y (i, i + a)cf^Ci+ay + Jaf ^ Si ■ Sj 

ia77' (i,j) 

+ EjT Xl[2?i(^ia - ^ib) + qli (81) 

i 

where ni^= Cj Ci^ and the hopping amphtude is changed as 

^77'(iJ) = -^(6)7*?^^'-^" (^j)r?'7'- (82) 
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In order to characterize the shape of the zigzag path, it is useful to formulate 
the change of the phase. For the purpose, we use the concept of "the Berry-phase 
connection" and define "the winding number" by following Ref . |149j . The phase-dressed 
operator, q^, naturally introduces a connection form, the Berry-phase connection, as 

uno] 

' z(0|craVctjO) ■ dr 2(0|craV4|0) " dr 



z(0|crbVctjO) • dr 2(0|crbV4|0) ■ dr 



Adr = 

1 / V^-dr -zV^ • dr 



where |0) denotes the vacuum state. 

Now we consider a 2D sheet of the JT crystal. In two space dimensions, there is a 
topologically conserved current for an arbitrary vector field v as 

f = ^J2'"''^^''>^^ (84) 

where j is the current and e'^'^^ is the antisymmetric Levi-Civita tensor. This current 
obeys the continuity equation, expressed as 

J2d,j'' = 0. (85) 

Therefore, if j^ and j^ are zero at the boundary of a closed surface S, the quantity Q 
(Chern number), defined by 

Q^ [dh f (86) 

Js 

is conserved, where indices 1 and 2 indicate the space-components and indicates the 

time-component of a vector, respectively. 

For time-independent solutions with which we are now concerned, we have doj'^ = 0. 

Thus, Q is conserved for an arbitrary surface S. Substituting a = Tt{A) = V^ for v 

into Eq. (|84p . we obtain the topologically conserved quantity, or "the winding number" 

as 



/ d^r f = — d^r{dia2 - ^201) 



w 

Is ^" Js 

= — (b drsL= — (t) dr-V^ = m, (87) 

where m is an integer representing the number for the twisting-around of the JT 
distortions along a path C enclosing S. Because of the conserved nature, we will used 
the winding number w to label a state hereafter. 

In the system with zigzag AF structure, C is considered to be a closed loop for the 
ID path in the periodic boundary condition. In this case, the winding number W may be 
decomposed into two terms as w=Wg+Wt- The former, Wg, is the geometric term, which 
becomes (1) corresponding to the periodic (anti-periodic) boundary condition in the 



Orbital ordering phenomena in d- and f -electron systems 36 




Figure 7. A typical building block for a zigzag ID FM path for an eg electron with 
orbital ordering. 



Cg-electron wavefunction. The discussion on the kinetic energy leads us to conclude that 
the state with 'Wg=0 has lower energy than that with 'Wg=l for x > 1/2, in agreement 
with the two-site analysis p^Slj. Thus, Wg is taken as zero hereafter. 

In order to understand that only the number of corners in the zigzag path, N^, 
determines the topological term Wt, let us consider the transfer of a single Cg electron 
along the path shown in Fig. [7| As mentioned above, on the straight-line part in 
the a;-(|/-)direction, the phase is fixed at 6'x=27r/3 (6'y=47r/3), because the Cg-electron 
orbital is polarized along the transfer direction. Thus, Wt does not change on the 
straight-line part of the path. However, when the electron passes the vertex a (/?), the 
phase changes from 6'x to 9y {9y to ^x), indicating that the electron picks up a phase 
change of 27r/3 (47r/3). Since these two vertices appear in pairs, Wt{=w) is evaluated 
as Wt=(iVv/2)(27r/3+47r/3)/(27r)=A^v/2. The phases at the vertices are assigned as an 
average of the phases sandwiching those vertices, ^a=vr and 6'^=0, to keep Wg invariant. 
Then, the phases are determined at all the sites, once ^xj 9y, 6a, and 6j3 are known. 

Finally, we note that the problem in the zigzag one-dimensional chain provided us 
with a typical example to better understand the importance of the additional factor 
g*€i/2 jj^ front of the 2x2 SU(2) unitary matrix to generate the phase dressed operator 
at each site. As clearly shown above, the "a" and "b" orbitals should be chosen as 
"a" =1/2 _ ^2 ^^^ "b"=3x2 - r2 at site 2, and ''a''=z^ - x^ and "b"=3y2 - r^ at site 
4, respectively. Namely, ^2=27r/3 and ^4=47r/3. The reason for these choices of ^i is 
easily understood due to the fact that the orbital tends to polarize along the hopping 
direction to maximize the overlap. Thus, to make the Hamiltonian simple, it is useful 
to fix the orbitals at sites 2 and 4 as ^2=2vr/3 and ^4=47r/3. Here, the phase factor 6*^'/^ 
in the basis function is essential to reproduce exactly the same solution as obtained in 
the discussion above. As already mentioned, in a single-site problem, this phase factor 
can be neglected, since it provides only an additional phase to the whole wave function. 
However, if the Cg-electron starts moving from site to site, the accumulation of the phase 
difference between adjacent sites does not lead just to an additional phase factor to the 
whole wave function. In fact, if this additional phase is accidentally neglected, the band 
structure will shift in momentum space as k—>-k + vr, indicating that the minimum of 
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the lowest-energy band is not located at ^=0, but at k=7i, as already pointed out by 
Koizumi et al. [148]. Of course, this can be removed by the redefinition of k by including 
"the crystal momentum" , but it is not necessary to redefine k, if the local phase factors 
are correctly included in the problem. 

3.2.4- Stability of zigzag structure Concerning the stability of the zigzag AF phase, 
from the results in Fig. IHl we can understand the following two points: (i) The zigzag 
structure of E- and CE-type shows the lowest energy compared with other zigzag paths 
with the same periodicity and compared with the straight ID path, (ii) The energy 
of the zigzag AF phase becomes lower than that of the FM or other AF phases in the 
parameter region of Jaf around JAF~0.1to- However, in the calculation, we have just 
assumed the periodicity of four lattice spacing, but it is unclear whether such period 
actually produces the global ground state or not. As emphasized above, it is true that 
the zigzag ID FM chain has a large band-gap, but this fact does not guarantee that this 
band-insulating phase is the lowest-energy state. In other words, we cannot exclude the 
possibihty that the zigzag structure with another periodicity becomes the global ground 
state. 

Unfortunately, it is quite difficult to carry out the direct comparison among the 
energies for all possible states, since there are infinite possibilities for the combinations 
of hopping directions. Instead, to mimic the periodic change of the phase 0a in the 
hopping process, let us imagine a virtual situation in which a JT distortion occurs in 
the ID Cg-electron system, by following Koizumi at al. jl47j . To focus on the effect of the 
local phase, it is assumed that the amplitude of the JT distortion g^ is independent of the 
site index, i.e., qi = q, and only the phase ^i is changed periodically. For simplicity, the 
phase is uniformly twisted with the period of M lattice spacings, namely, C,j=Jx{2tt)/M 
for l<j<M. Since the periodic change of the hopping direction is mimicked by the phase 
change of the JT distortion, tfj^ is simply taken as the unit matrix to^/ii/ to avoid the 
double-counting of the effect of the phase change. If the potential amphtude is written 
as f =2g£'jT, the Hamiltonian for the present situation is given by 

H= -to ^(49a + clbCjb + h.c.) 

+ vY^ [sin ^i (4^Qb + cJt^Qa) + COS ^i (cJ^Qa - cJ^Qb) ] , (88) 

i 

where the spinless Cg-electron operator is used, since the ID FM chain is considered 
here. The potential term for the JT distortion is ignored, since it provides only a 
constant energy shift in this case. By using the transformation Eqs. (f77j) and (f75j) . the 
Hamiltonian is rewritten as 

H=-to $:[e^«-^^)/^{cos ^(49a + 4£.b) 

+ sin ' ^ (49b - 4cia)} + h.c] + V 5]](4cia - 4cib). (89) 
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The Hainiltonian in inomentum space is obtained by the Fourier transform as 

^ = XI t^'^' cos(vr/M)(4a4a + 4bCfcb) + isk sin(7r/M)(4aCfcb - 4bCfca)] 

k 

+ ^Xl(^fcaCfca-4bCfcb), (90) 

k 

where 5^=— 2to cos /c, Sk=—2tQsmk, and the periodic boundary condition is imposed. 
Note that in this expression, k is the generahzed quasimomentum, redefined as 
k — n/M — > k, to incorporate the additional phase ir/M which appears to arise from 
a fictitious magnetic field (see Ref. |148j ). The eigenenergies are easily obtained by 
diagonalization as 



E^ = Sk cos(7r/M) ± ^v^ + sl^ixi^i'n/M) 



?± 



= (l/2)[£fc+7r/M + ^k-n/M ± J ^v'^ + {Sk+w/M — ^k-K/uY ]• (91) 

Since this is just the coupling of two bands, ek+n/M and ek-ixjM-, it is easily understood 
that the energy gain due to the opening of the bandgap is the best for the filling of 
n=2/M , where ra=l— x with doping x. In other words, when the periodicity M is equal 
to 2/ra, the energy becomes the lowest among the states considered here with several 
possible periods. Although this is just a proof in an idealized special situation, it is 
believed that it captures the essence of the problem. 

Here the effect of the local phase factor e*^'/^ should be again noted. If this factor 
is dropped, the phase n/M due to the fictitious magnetic field disappears and the 
eigenenergies are given by the coupling of ^fc+Tr+Tr/Af and Ek+TT-Tr/M, which has been 
also checked by the computational calculation. This "vr" shift in momentum space 
appears at the boundary, modifying the periodic boundary condition to anti-periodic, 
even if there is no intention to use anti-periodic boundary condition. Of course, this is 
avoidable when the momentum k is redefined as fc + vr ^ A;, as pointed out in Ref. |148j . 
However, it is natural that the results for periodic boundary condition are obtained 
in the calculation using periodic boundary condition. Thus, also from this technical 
viewpoint, it is recommended that the phase factor e*^'/^ is added for the local rotation 
in the orbital space. 

3.2.5. Summary In summary, at x=0.5, the CE-type AF phase can be stabilized even 
without the Coulomb and/or the JT phononic interactions, only with large Hund and 
finite Jaf couplings. We have also pointed out the appearance of E-type phase due 
to the same mechanism. Of course, as we will see in the next subsection, Coulombic 
and JT phononic interactions are needed to reproduce the charge and orbital ordering. 
However, as already mentioned in the above discussion, because of the special geometry 
of the one-dimensional zigzag FM chain, for instance, at x=0.5, it is easy to imagine that 
the checkerboard type charge-ordering and (Sx^ — r^/3|/^ — r^) orbital-ordering pattern 
will be stabilized. Furthermore, the charge confinement in the straight segment, i.e., 
sites 2 and 4 in Fig. IHl^b), will naturally lead to charge stacking along the z-axis, with 
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stability caused by the special geometry of the zigzag structure. Thus, the complex 
spin-charge-orbital structure for half-doped manganites can be understood intuitively 
simply from the viewpoint of its band-insulating nature. 

3.3. Spin, charge, and orbital ordering 

Now we review the theoretical results on spin, charge, and orbital ordering in undoped 
and doped manganites on the basis of realistic models. The Hamiltonian mainly used 
here is two-orbital double exchange model strongly coupled with Jahn- Teller phonons, 
explicitly given by 

H = -f^kin + -f^Hund + -f^inter-sitc + -f^cl-plr (92) 

In the infinite limit for Jh, we can further simplify the model into the following form. 

+ EjT Y, [2(9ii^i + q2ir,i + qsiT.i) + [3ql, + ql + ql] , (93) 

i 

where ni= 'Yli^y^i-y^ii ^^^ /3=khr/k,]T, the ratio of spring constants of breathing and 
Jahn- Teller phonons. Concerning the value of /3, from experimental results and band- 
calculation data for the energy of breathing and Jahn- Teller modes |152j , it is estimated 
as /? ~ 2 for manganites. It is convenient to introduce non-dimensional electron-phonon 
coupling constant A as 



A = V2EjT/to = g/Vk^o- (94) 

Here we simply drop the Coulomb interaction terms, but the reason will be discussed 
in the next subsection. 

This model is analyzed by using both the numerical techniques (Monte Carlo 
simulation and relaxation method) and mean-field approximation. Note that in the 
numerical simulations, depending on the non-cooperative and cooperative treatments, 
the variables are angles 6i and 0i which specifies t2g spin directions and coordinates 
{q} and oxygen positions {u}, respectively. In any case, all variables are classical and 
thus, there is no essential problems to perform almost exactly the simulation, within 
the limit of the power of computers. Recently, there has been an important progress in 
the simulation for the electron systems coupled with classical variables |153| 11541 I155J . 
In particular, Motome and Furukawa have developed the efficient simulation technique 
for the acceleration of the calculation and the increase of the precision J1541I15,'^ . 

3.3.1. Effect of Coulomb interaction Let us discuss briefly the effect of the Coulomb 
interaction. Since we consider the strong Hund's rule interaction between Cg electron 
and t2g localized spins, Cg electron spins tend to array in a site and thus, the effect of 
intra-orbital Coulomb interaction is automatically suppressed. However, inter-orbital 
Coulomb interaction still works between electrons with the same spin. Here we explain 
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the reason why we ignore the on-site Coulomb interaction. The effect of inter-site 
Coulomb interaction is discussed in the stabilization mechanism of charge stacking in 
the x=0.5 CE-type structure. 

In the spinless model, the inter-orbital Coulomb interaction term is written by 

-f^cl-cl = U'^ ^ia^ib, (95) 

i 

where ni^= c[ Cj-y and the present U' means f/' — J in the standard notation for the 
on-site Coulomb interaction. We also consider the inter-site Coulomb interaction term, 
given by 

^inter-sitc = ^ J^ '^i^j' (96) 

where ni=Y^ c[ Ci^. In order to understand the ignorance of on-site Coulomb interaction 
term, it is quite instructive to consider the mean-field approximation. As for the detail 
of the formulation, readers can refer the original paper |156j and the previous review 
j3]. Here we show the result of the mean- field Hamiltonian. 

HmF = — 2_^ -Di.i+ai^^'C;^Ci+a7' + <^AF 2_^ ^i ' ^j 

ia77' (i,j) 

+ EjT Xl[-2((^xi)^xi + {rzi)r,i) + (r,i)2 + (r,;)^] 

i 

+ J2m'/2){n,) + F^(r^i+a)](r^i - {n,)/2), (97) 

i a 

where (■ ■ ■) denotes the average value. The renormalized JT energy is given by 

EjT = ^jT + f/74, (98) 

and the renormalized inter-orbital Coulomb interaction is expressed as 

U' = U' - 4Eto, (99) 

where £'br=5'^/(2/i;br)- Physically, the former relation indicates that the JT energy is 
effectively enhanced by U'. Namely, the strong on-site Coulombic correlation plays the 
same role as that of the JT phonon, at least at the mean-field level, indicating that it 
is not necessary to include U' explicitly in the models, as emphasized in Ref. |156j . See 
also Ref. [21]. The latter equation for U' means that the one-site inter-orbital Coulomb 
interaction is effectively reduced by the breathing-mode phonon, since the optical-mode 
phonon provides an effective attraction between electrons. The expected positive value 
of U' indicates that Cg electrons dislike double occupancy at the site, since the energy 
loss is proportional to the average local electron number in the mean-field argument. 
Thus, to exploit the gain due to the static JT energy and avoid the loss due to the 
on-site repulsion, an Cg electron will singly occupy a given site. 
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Figure 8. (a) The four spin arrangements for Ferro, A-AF, C-AF, and G-AF. (b) 
Total energy vs Jaf on a 1? cluster at low temperature with Jh^S^o and A=1.5. The 
results were obtained using Monte Carlo and relaxational techniques, with excellent 
agreement among them, (c) Orbital order corresponding to the A-type AF state. For 
more details, the reader should consult Ref. pS] . 



3.3.2. x=0 First let us consider the mother material LaMnOs with one Cg electron per 
site. This material has the insulating AF phase with A-type AF spin order, in which 
t2g spins are ferromagnetic in the a-h plane and AF along the c-axis. For the purpose 
to understand the appearance of A-AF phase, it is enough to consider a 2 x 2 x 2 cube 
as a minimal cluster for undoped manganites. Results in a larger size cluster will be 
discussed later. 

Recent investigations by Hotta et al. jS^j have shown that, in the context of the 
model with Jahn- Teller phonons, the important ingredient to understand the A-type AF 
phase is Jaf, namely by increasing this coupling from 0.05 to larger values, a transition 
from a FM to an A-type AF exists The relevance of Jahn- Teller couplings at x=0.0 
has also been remarked in Ref. [27j . This can be visualized easily in Fig. |H1 where the 
energy vs. Jaf at fixed intermediate A and Jh is shown. Four regimes are identified: 
FM, A-AF, C-AF, and G-AF states that are sketched also in that figure. The reason 
is simple: As Jaf grows, the tendency toward spin AF must grow, since this coupling 
favors such an order. If Jaf is very large, then it is clear that a G-AF state must be 
the one that lowers the energy, in agreement with the Monte Carlo simulations. If Jaf 
is small or zero, there is no reason why spin AF will be favorable at intermediate A 
and the density under consideration, and then the state is ferromagnetic to improve the 
electronic mobility. It should be no surprise that at intermediate Jaf, the dominant state 
is intermediate between the two extremes, with A-type and C-type antiferromagnetism 
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becoming stable in intermediate regions of parameter space. 

It is interesting to note that similar results regarding the relevance of Jaf to stabilize 
the A-type order have been found by Koshibae et al. fH] in a model with Coulomb 
interactions. An analogous conclusion was found by Solovyev, Hamada, and Terakura 
[ini lEj and Ishihara et al. [T2j. Betouras and Fujimoto |157j . using bosonization 
techniques for the one- dimensional one-orbital model, also emphasized the importance of 
Jaf, similarly as did by Yi, Yu, and Lee based on Monte Carlo studies in two dimensions 
of the same model |158j . The overall conclusion is that there are clear analogies between 
the strong Coulomb and strong Jahn- Teller coupling approaches. Actually in the mean- 
field approximation, it was shown by Hotta, Malvezzi, and Dagotto |15(ij that the 
influence of the Coulombic terms can be hidden in simple redefinitions of the electron- 
phonon couplings (see also Ref. J2l]). In our opinion, both approaches (Jahn- Teller and 
Coulomb) have strong similarities and it is not surprising that basically the same physics 
is obtained in both cases. Actually, Fig. 2 of Maezono, Ishihara, and Nagaosa jl59j 
showing the energy vs. Jaf in mean-field calculations of the Coulombic Hamiltonian 
without phonons is very similar to our Fig. |HI^b), aside from overall scales. On the other 
hand, Mizokawa and Fujimori fUEj states that the A-type AF is stabilized only when 
the Jahn- Teller distortion is included, namely, the FM phase is stabilized in the purely 
Coulomb model, based on the unrestricted Hartree-Fock calculation for the d-p model. 

The issue of what kind of orbital order is concomitant with A-type AF order is an 
important matter. This has been discussed at length by Hotta et al. |23 , and the final 
conclusion, after the introduction of perturbations caused by the experimentally known 
difference in lattice spacings between the three axes, is that the order shown in Fig. |Hl^c) 
minimizes the energy. This state has indeed been identified in resonant X-ray scattering 
experiments JU], and it is quite remarkable that such a complex pattern of spin and 
orbital degrees of freedom indeed emerges from mean-field and computational studies. 
Studies by van den Brink et al. |2ni using purely Coulombic models arrived at similar 
conclusions. The orbital ordering has been also captured from the viewpoint of orbital 
density wave state by Koizumi et al. |147| I148J . The similar discussion has been done 
recently by Efremov and Khomskii |16Uj . 

Why does the orbital order occur here? In order to respond to this question, it 
is quite instructive to consider the situation perturbatively in the electron hopping. A 
hopping matrix only connecting the same orbitals, with hopping parameter t, is assumed 
for simplicity. The energy difference between Cg orbitals at a given site is -Ejt, which is a 
monotonous function of A. For simplicity, in the notation let us refer to orbital uniform 
(staggered) as orbital "FM" ("AF"). Case (a) corresponds to spin FM and orbital AF: 
In this case when an electron moves from orbital a on the left to the same orbital on the 
right, which is the only possible hopping by assumption, an energy of order Ej^ is lost, 
but kinetic energy is gained. As in any second order perturbative calculation the energy 
gain is then proportional to t'^/Ej'j<. In case (b), both spin and orbital FM, the electrons 
do not move and the energy gain is zero (again, the nondiagonal hoppings are assumed 
negligible just for simplicity). In case (c), the spin are AF but the orbitals are FM. This 
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Figure 9. (a) Ground-state phase diagram for undoped manganites by using the 
4x4x4 lattice. Sohd curves denote the mean-field results, while solid circles indicate 
the result for optimization, (b) Spin structure for E-AF phase in three-dimensional 
environment, (c) Fcrmi-surface lines of 2D Cg electron system at x=0.0. 



is like a one orbital model and the gain in energy is proportional to t^/(2JH). Finally, in 
case (d) with AF in spin and orbital, both Hund and orbital splitting energies are lost 
in the intermediate state, and the overall gain becomes proportional to t^/(2JH + -Ejt)- 
As a consequence, if the Hund coupling is larger than Ej^, then case (a) is the best, 
as it occurs at intermediate Ejic values. Then, the presence of orbital order can be 
easily understood from a perturbative estimation, quite similarly as done by Kugel and 
Khomskii in their pioneering work on orbital order |161j . 

Here readers may have a question in their mind. Where is the E-type phase 
emphasized in the previous subsection? In order to respond to this question, it is 
necessary to treat a larger-size cluster. In Fig. El^a), we show the phase diagram of 
undoped manganites with the direct comparison between the mean-field and numerical 
results in the 4x4x4 lattice |SI]. Solid curves are depicted from the cooperative mean- 
field approximation, while solid circles denote the result for optimization both for t2g 
spin directions and oxygen positions. The good agreements clearly indicate that the 
present mean-field procedure works quite well for undoped manganites. 

Now the phase diagram includes six phases, but there is a clear separation around 
at A ~ 1.5. For A > 1.5, there occurs a chain of transitions in the order of FM, A-AF, 
C-AF, and G-AF phases with increasing Jaf, already obtained in 2x2x2 calculations. 
Note that the boundary curve always indicates the first order transition. The present 
result shows that size effects are small in undoped strongly-coupled manganites, which is 
intuitively reasonable. The spin arrangement for each phase is shown in Fig. |Hl^a). Note 
that the FM phase is concomitant with orbital ordering (00), which will be discussed 
later, and this FM/00 phase is considered to be insulating. 

On the other hand, in the weak or intermediate coupling region for A < 1.5, there 
is a transition between the FM orbital-disordered (OD) phase and E-type AF phase. 
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The spin arrangement for E-AF phase is shown in Fig. Mjo). Along the zigzag chains, 
t2g spins order ferromagnetically, but they are antiparallel perpendicular to the zigzag 
direction. This is just the new AF phase, suggested by recent experiments on HoMnOa. 
As suggested in Fig.E^b), the spin directions are reversed from plane to plane. Note that 
the orbital structure in the E-AF phase is the same as that of the A-AF phase, namely, 
the staggered pattern of {3x^—r^)- and (3y^— r^)-like orbitals. It is easily understood 
that the orbital ordering is closely related to the cooperative Jahn- Teller distortion in 
undoped manganites and such a cooperative effect should be very strong irrespective of 
the t2g spin configuration. 

Note that near A~1.5, which is a realistic value for manganites, the A-AF phase 
is adjacent to the E-type state. This region could correspond to the actual situation 
observed in experiments for RMnOa: When the ionic radius of the R-site decreases, 
a Neel temperature Tn of the A-AF phase decreases as well, and eventually the E- 
AF phase is stabilized for R=Ho. Another interesting point of the phase diagram is 
that the E-type spin arrangement is the ground-state for a wide range of Jaf, even at 
A=0, indicating that the coupling with JT phonons is not a necessary condition for its 
stabilization. As pointed out in the previous subsection, the E-type phase is stable due 
to the zigzag geometry of the FM chains that induce a band-insulator. Namely, E-AF 
phase is always insulating irrespective of the value of A. 

Concerning the appearance of the E-AF phase, Kimura et al. have explained it on 
the basis of a frustrated spin system with FM nearest-neighbor and AF next-nearest- 
neighbor interactions within the Mn02 plane jSH]- They have found that the staggered 
orbital order associated with the GdFeOs-type distortion induced the anisotropic next- 
nearest-neighbor interaction, leading to unique sinusoidal and up-up-down-down AF 
order, i.e., E-type phase, in undoped manganites. In a conceptual level, the spin model 
is considered to be obtained in the strong coupling limit of the Cg-orbital degenerate 
double-exchange model. Thus, the band-insulating picture for the appearance of E-type 
phase in the present scenario is complementary to the result of Kimura et al., in the sense 
that the weak-coupling state is continuously connected to that in the strong-coupling 
limit. 

In addition to the explanation of the A-AF of LaMnOs and E-AF of HoMnOs, 
Kimura et al. have also examined systematically the magnetic and orbital structures in 
a series of RMnOs as a function of tr, the radius of rare-earth ion R. They have pointed 
out that the effect on the crystal structure by decreasing tr appears as the enhancement 
of the GdFeOs-type distortion, indicating the shortening of oxygen-oxygen distance. 
Then, the superexchange interaction between next-nearest-neighbor sites is enhanced 
due to the shortened oxygen-oxygen path, leading to the frustrated spin model with the 
competition between FM nearest-neighbor and AF next-nearest-neighbor interactions 
By analyzing the frustrated spin model on the staggered orbital-ordered background, 
Kimura et al. have explained the phase diagram of RMnOs. It is considered that 
the phase diagram can be also understood from the band-insulating picture, but for 
comparison with actual materials, it is necessary to include the effect of the GdFeOa- 
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type distortions, which has not been considered in the present modeL 

Near the transition region between A- and E-type AF phases, Salafranca and 
Brey have mentioned the importance of the competition between the nearest neighbor 
AF superexchange interaction and the double exchange induced long-range FM 
interaction |l(i2j . They concluded that such competition results in the appearance of 
incommensurate phases. These phases consist of a periodic array of domain walls. 

As discussed above, in the strong-coupling region, FM/00 insulating state appears, 
but when A decreases, 00 disappears and instead, an OD phase is observed. This is 
considered as a metallic phase, as deduced from the result of the density of states. Note 
that this metallic OD/FM phase is next to the insulating E-AF phase for A < 1.5, which 
is a new and important result in the study of undoped manganites jl63j . Namely, the 
competition between FM metallic and insulating phases is at the heart of the CMR 
phenomena, and then, by tuning experimentally the lattice parameters in RMnOa it 
may be possible to observe the magnetic-field induced metal-insulator transition even 
in undoped manganites. 

Let us consider the reason why the metallic phase can exist even at half-filling. 
To clarify this point, it is quite useful to depict the Fermi-surface lines. As shown in 
Fig. Ofc), the nesting vector is (vr, 0) or (0, vr), not (vr, tt). These nesting vectors are not 
compatible with the staggered orbital ordering pattern that is stabilized increasing A. 
This is one of the remarkable features of the multiorbital Cg-electron system, which is not 
specific to two dimensionality. In fact, in the results for the three-dimensional (3D) case, 
we also observe the signal of the metal- insulator transition at a finite value of A. In this 
case, the orbital ordering pattern becomes very complicated, but the pattern repeats 
periodically on lattice larger than 2x2x2. In the 3D case, an intrinsic incompatibility 
between the Fermi surface and the orbital ordering pattern is also found. Even without 
invoking the numerical results discussed before, the qualitative arguments related with 
the nesting effects in i^kin incompatible with staggered orbital ordering strongly suggests 
the presence of a metallic phase in two and three dimensions at small A. 

3.3.3. x=0.5 Now let us move to another important doping x=0.5. For half-doped 
perovskite manganites, the so-called CE-type AF phase has been established as the 
ground state in the 1950 's. This phase is composed of zigzag FM arrays of t2g- 
spins, which are coupled antiferromagnetically perpendicular to the zigzag direction. 
Furthermore, the checkerboard-type charge ordering in the x-y plane, the charge stacking 
along the z-axis, and (3x^ — r'^ / "iy"^ — r'^) orbital ordering are associated with this phase. 
A schematic view of CE-type structure with charge and orbital ordering is shown in 
Fig. [TUT a) for 2D case. In 3D, this patters repeats along the z-axis by keeping charge 
and orbital structure, but changing spin directions. 

Although there is little doubt that the famous CE-state of Goodenough is indeed the 
ground state of x=0.5 intermediate and low bandwidth manganites, only very recently 
such a state has received theoretical confirmation using unbiased techniques, at least 
within some models. In the early approach of Goodenough it was assumed that the 
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Figure 10. (a) Schematic view of CE-type structure at x=0.5 for 2D case, (b) Monte 
Carlo energy per site vs Jaf at density x=0.5, A=1.5, low temperature r=l/100, and 
Jh=oo, using the two-orbital model in two dimensions with non-cooperative Jahn- 
Teller phonons. As for AF2, see Fig. |H{c). (c) Phase diagram in the plane A- Jaf at 
x=0.5, obtained numerically using up to 8x8 clusters. All transitions arc of first-order. 
The notation is the standard one (CD — charge disorder, CO = charge order, 00 = 
orbital order, OD = orbital disorder). Results are reproduced from Ref. [32] where 
more details can be found. 



charge was distributed in a checkerboard pattern, upon which spin and orbital order 
was found. But it would be desirable to obtain the CE-state based entirely upon a more 
fundamental theoretical analysis, as the true state of minimum energy of a well-defined 
and realistic Hamiltonian. If such a calculation can be done, as a bonus one would find 
out which states compete with the CE-state in parameter space, an issue very important 
in view of the mixed-phase tendencies of Mn-oxides, which cannot be handled within 
the approach of Goodenough. 

One may naively believe that it is as easy as introducing a huge nearest-neighbor 
Coulomb repulsion V to stabilize a charge-ordered state at x=0.5, upon which the 
reasoning of Goodenough can be applied. However, there are at least two problems with 
this approach |164j . First, such a large V quite likely will destabilize the FM charge- 
disordered state and others supposed to be competing with the CE-state. It may be 
possible to explain the CE-state with this approach, but not others also observed at 
x=0.5 in large bandwidth Mn-oxides. Second, a large V would produce a checkerboard 
pattern in the three directions. However, experimentally it has been known for a long 
time [Sj that the charge stacks along the 2;-axis, namely the same checkerboard pattern 
is repeated along 2;-axis, rather than being shifted by one lattice spacing from plane 
to plane. A dominant Coulomb interaction V cannot be the whole story for x=0.5 
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low-bandwidth manganese oxides. 

The nontrivial task of finding a CE-state without the use of a huge nearest-neighbors 
repulsion has been recently performed by Yunoki, Hotta, and Dagotto |S2], using the 
two-orbital model with strong electron Jahn- Teller phonon coupling. The calculation 
proceeded using an unbiased Monte Carlo simulation, and as an output of the study, 
the CE-state indeed emerged as the ground-state in some region of coupling space. 
Typical results are shown in Figs. [TUT b) and lTUT c). In part (b) the energy at very low 
temperature is shown as a function of Jaf at fixed density x=0.5, Jh=oo for simplicity, 
and with a robust electron-phonon coupling A=1.5 using the two orbital model ifjx 
At small Jaf, a FM phase was found to be stabilized, according to the Monte Carlo 
simulation. Actually, at Jaf=0.0 it has not been possible to stabilize a partially AF- 
state at x=0.5, namely the states are always ferromagnetic at least within the wide range 
of A's investigated (but they can have charge and orbital order). On the other hand, as 
Jaf grows, a tendency to form AF links develops, as it happens at x=0.0. At large Jaf 
eventually the system transitions to states that are mostly antiferromagnetic, such as the 
so-called "AF(2)" state of Fig. fTUTb) (with an up-up-down-down spin pattern repeated 
along one axis, and AF coupling along the other axis), or directly a fully AF-state in 
both directions. 

However, the intermediate values of Jaf are the most interesting ones. In this case 
the energy of the two-dimensional clusters become fiat as a function of Jaf suggesting 
that the state has the same number of FM and AF links, a property that the CE-state 
indeed has. By measuring charge-correlations it was found that a checkerboard pattern 
is formed particularly at intermediate and large A's, as in the CE-state. Finally, after 
measuring the spin and orbital correlations, it was confirmed that indeed the complex 
pattern of the CE-state was fully stabilized in the simulation. This occurs in a robust 
portion of the A- Jaf plane, as shown in Fig. [TUT c). The use of Jaf as the natural 
parameter to vary in order to understand the CE-state is justified based on Fig. irUT c). 
since the region of stability of the CE-phase is elongated along the A-axis, meaning that 
its existence is not so much dependent on that coupling but much more on Jaf itself. 
It appears that some explicit tendency in the Hamiltonian toward the formation of AF 
links is necessary to form the CE-state. If this tendency is absent, a FM state if formed, 
while if it is too strong an AF-state appears. The x=0.5 CE-state, similar to the A-type 
AF at x=0.0, needs an intermediate value of Jaf for stabilization. The stability window 
is finite and in this respect there is no need to carry out a fine tuning of parameters to 
find the CE phase. However, it is clear that there is a balance of AF and FM tendencies 
in the CI5-phase that makes the state somewhat fragile. 

Note that the transitions among the many states obtained when varying Jaf are 
all of first order, namely they correspond to crossings of levels at zero temperature. 
The first-order character of these transitions is a crucial ingredient of the recent scenario 
proposed by Moreo et al. jl65j involving mixed-phase tendencies with coexisting clusters 
with equal density. Recently, first-order transitions have also been reported in the one- 
orbital model at x=0.5 by Alonso et al. |166| 1167] . as well as tendencies toward phase 
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Figure 11. (a) Energy per site as a function of Jaf for A=1.6 and Jh=oo for Hjt 
using a 4x4x4 lattice. The curves denote the mean- field results and the solid symbols 
indicate the energy obtained by the relaxation method. Thick solid, thick broken, thin 
broken, thick dashed, thin dashed, thin broken, and thin solid lines denotes FM, A- 
type, CE-type with WC structure, charge-stacked CE-type, C-type, and G-type states, 
respectively. Note that the charge-stacked CE-state is observed in experiments, (b) 
Phase diagram in the {Jaf,V) plane for 4x4x4 lattice. Note that the charge-stacked 
structure along the z-axis can be observed only in the CE-type AF phase. Results are 
reproduced from Ref. |15()| where more details can be found. 



separation. 

Let us address now the issue of charge-stacking (CS) along the 2;-axis. For this 
purpose simulations using 3D clusters were carried out |156j . The result for the energy 
vs. Jaf is shown in Fig. llir a). with Jh=oo and A=1.6 fixed. The CE-state with charge- 
stacking has been found to be the ground state on a wide Jaf window. The reason 
that this state has lower energy than the so-called "Wigner-crystal" (WC) version of 
the CE-state, namely with the charge spread as much as possible, is once again the 
influence of Jaf- With a charge stacked arrangement, the links along the 2-axis can all 
be simultaneously antiferromagnetic, thereby minimizing the energy. In the WC-state 
this is not possible. 

It should be noted that this charge stacked CE-state is not immediately destroyed 
when the weak nearest-neighbor repulsion V is introduced to the model, as shown in 
Fig. ITlT b). obtained in the mean- field calculations by Hotta, Malvezzi, and Dagotto 
[15()j . If V is further increased for a realistic value of Jaf, the ground state eventually 
changes from the charge stacked CE-phase to the WC version of the CE-state or the C- 
type AF phase with WC charge ordering. As explained above, the stability of the charge 
stacked phase to the WC version of the CE-state is due to the magnetic energy difference. 
However, the competition between the charge-stacked CE-state and the C-type AF phase 
with the WC structure is not simply understood by the effect of Jaf, since those two 
kinds of AF phases have the same magnetic energy. In this case, the stabilization of the 
charge stacking originates from the difference in the geometry of the one- dimensional 
FM path, namely a zigzag-path for the CE-phase and a straight-line path for the C-type 
AF state. As discussed above, the energy for Cg electrons in the zigzag path is lower than 
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that in the straight-hne path, and this energy difference causes the stabihzation of the 
charge stacking. In short, the stabihty of the charge-stacked structure at the expense 
of V is supported by "the geometric energy" as well as the magnetic energy. Note that 
each energy gain is just a fraction of to- Thus, in the absence of other mechanisms to 
understand the charge- stacking, another consequence of this analysis is that V actually 
must be substantially smaller than naively expected, otherwise such a charge pattern 
would not be stable. In fact, estimations given by Yunoki, Hotta, and Dagotto j22] 
suggest that the manganites must have a large dielectric function at short distances (see 
Ref. |168j ) to prevent the melting of the charge-stacked state. 

Note also that the mean-field approximations by Hotta, Malvezzi, and Dagotto 
|156j have shown that on-site Coulomb interactions U and U' can also generate a two- 
dimensional CE-state, in agreement with the calculations by van den Brink et al. |142j . 
Then, we believe that strong Jahn- Teller and Coulomb couplings tend to give similar 
results. This belief finds partial confirmation in the mean-field approximations of Hotta, 
Malvezzi, and Dagotto |15()j . where the similarities between a strong A and (f/, U') were 
investigated. Even doing the calculation with Coulombic interactions, the influence of 
Jaf is still crucial to inducing charge- stacking. The importance of this parameter has 
also been remarked by Mathieu, Svedlindh and Nordblad |169j based on experimental 
results. 

Many other authors carried out important work in the context of the CE-state at 
x=0.5. For example, with the help of Hartree-Fock calculations, Mizokawa and Fujimori 
|17()j reported the stabilization of the CE-state at x=0.5 only if Jahn- Teller distortions 
were incorporated into a model with Coulomb interactions. This state was found to 
be in competition with a uniform FM state, as well as with an A-type AF-state with 
uniform orbital order. In this respect the results are very similar to those found by 
Yunoki, Hotta and Dagotto jH2] using Monte Carlo simulations. In addition, using a 
large nearest-neighbor repulsion and the one-orbital model, charge ordering and a spin 
structure compatible with the zigzag chains of the CE state was found by Lee and Min at 
x=0.5 )171j . Jackeli, Perkins, and Plakida also obtained charge-ordering at x=0.5 using 
mean-field approximations and a large V [ITS] . Charge-stacking was not investigated 
by those authors. The CE-state in Pro.sCao.sMnOs was also obtained by Anisimov et 
al. using LSDA-I-U techniques |173j . 

3.3.4- x>0.5 In the previous subsection, the discussion focused on the CE-type AF 
phase at x=0.5. Naively, it may be expected that similar arguments can be extended to 
the regime x>l/2, since in the phase diagram for LaxCai-xMnOs, the AF phase has been 
found at low temperatures in the region 0.50<x<0.88. Then, let us try to consider the 
band-insulating phase for density x=2/3 based on the minimal model Eq. (jF)H|) without 
both the Jahn- Teller phononic and Coulombic interactions, since this doping is quite 
important for the appearance of the bi-stripe structure (see Refs. |174|I175] ). 

After several calculations for x=2/3, as reported by Hotta et al. |141j . the lowest- 
energy state was found to be characterized by the straight path, not the zigzag one. 
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Figure 12. (a) Path with w=l at x=l/2. Charge and orbital densities are calculated 
in the MFA for E}T—2t. At each site, the orbital shape is shown with its size in 
proportion to the orbital density, (b) The BS-structure path with w=2 at x=2/3. (c) 
The BS-structure path with w=3 at x=3/4. (d) The WC-structure path with w—l at 
x=2/3. (e) The WC-structure path with w^l at x=3/4. 



leading to the C-type AF phase which was also discussed in previous subsection. For 
a visual representation of the C-type state, see Fig. 4 of Ref. |176j . At first glance, the 
zigzag structure similar to that for x=0.5 could be the ground-state for the same reason 
as it occurs in the case of x=0.5. However, while it is true that the state with such a 
zigzag structure is a band-insulator, the energy gain due to the opening of the bandgap 
is not always the dominant effect. In fact, even in the case of x=0.5, the energy of the 
bottom of the band for the straight path is — 2to, while for the zigzag path, it is — \/3to. 
For x=l/2, the energy gain due to the gap opening overcomes the energy difference 
at the bottom of the band, leading to the band-insulating ground-state. However, for 
x=2/3 even if a band-gap opens the energy of the zigzag structure cannot be lower than 
that of the metallic straight-line phase. Intuitively, this point can be understood as 
follows: An electron can move smoothly along the one-dimensional path if it is straight. 
However, if the path is zigzag, "reflection" of the wavefunction occurs at the corner, and 
then a smooth movement of one electron is no longer possible. Thus, for small numbers 
of carriers, it is natural that the ground-state is characterized by the straight path to 
optimize the kinetic energy of the Cg electrons. 

However, in neutron scattering experiments a spin pattern similar to the CE-type 
AF phase has been suggested by Radaelli et al. |177j . In order to stabilize the zigzag 
AF phase to reproduce those experiments it is necessary to include the Jahn- Teller 
distortion effectively. As discussed by Hotta et al. |141j . a variety of zigzag paths 
could be stabilized when the Jahn- Teller phonons are included. In such a case, the 
classification of zigzag paths is an important issue to understand the competing "bi- 
stripe" vs. "Wigner- crystal" structures. The former has been proposed by Mori et 
al. |174[ 1175] . while the latter was claimed to be stable by Radaelli et al. |177j . As 
shown in the previous subsection, the shape of the zigzag structure is characterized by 
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Figure 13. (a) C- and E-type unit cells j8j. (b) The spin structure in the a-b plane 
at x=l/2. Open and solid symbols denote the spin up and down, respectively. The 
thick line indicates the zigzag FM path. The open and shaded squares denote the C- 
and E-type unit cells. At x=l/2, C-type unit cell occupies half of the two-dimensional 
plane, clearly indicating the "CE" type phase, (c) The spin structure at x=2/3 for 
Wigner-crystal type phase. Note that 66% of the two-dimensional lattice is occupied by 
C-type unit cell. Thus, it is called "C2/3Ei/3"-type AF phase, (d) The spin structure 
at x=2/3 for bi-stripe type phase. Note that 33% of the two-dimensional lattice is 
occupied by C-type unit cell. Thus, it is called "Ci/3E2/3"-type AF phase. 



the winding number w associated with the Berry-phase connection of an Cg-electron 
parallel-transported through Jahn- Teller centers, along zigzag one- dimensional paths. 
As shown above, the total winding number is equal to half of the number of corners 
included in the zigzag unit path. Namely, the winding number w is a good label to 
specify the shape of the zigzag one-dimensional FM path. 

After several attempts to include effectively the Jahn- Teller phonons, it was found 
that the bi-stripe phase and the Wigner crystal phase universally appear for w=x/{l—x) 
and w=l, respectively. Note here that the winding number for the bi-stripe structure 
has a remarkable dependence on x, reflecting the fact that the distance between adjacent 
bi-stripes changes with x. This x-dependence of the modulation vector of the lattice 
distortion has been observed in electron microscopy experiments |174l I175J . The 
corresponding zigzag paths with the charge and orbital ordering are shown in Fig. El In 
the bi-stripe structure, the charge is confined in the short straight segment as in the case 
of the CE-type structure at x=0.5. On the other hand, in the Wigner-crystal structure, 
the straight segment includes two sites, indicating that the charge prefers to occupy 
either of these sites. Then, to minimize the Jahn- Teller energy and/or the Coulomb 
repulsion, the Cg electrons are distributed with equal spacing. The corresponding spin 
structure is shown in Fig. ^1 A difference in the zigzag geometry can produce a 
significant different in the spin structure. The definitions for the C- and E-type AF 
structures jH] are shown in Fig. fT^ a) for convenience. At x=l/2, as clearly shown 
in Fig. IT!¥ b). half of the plane is filled by the C-type, while another half is covered 
by the E-type, clearly illustrating the meaning of "CE" in the spin structure of half- 
doped manganites. On the other hand, as shown in Figs. I13r c) and lT^ d). the bi-stripe 
and Wigner crystal structure have Ci„xEx-type and CxEi_x-type AF spin arrangements, 
respectively. Such zigzag-based AF structure has been discussed experimentally in single 
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Figure 14. Schematic views for spin, charge, and orbital ordering for (a) Wigner- 
crystal structure and and (b) bi-stripe type structures at x=2/3 in 3D environment. 
The open and soUd symbols indicate the spin up and down, respectively. The FM 
one-dimensional path is denoted by the thick line. The empty sites denote Mn^"*" ions, 



while the robes indicate the Mn'^+ ions in which 3x^ 
occupied. 



or 3y^ — r^ orbitals are 



layered manganites Ndi_xSri+xMn04 by Kimura et al. |178j . 

The charge structure along the 2;-axis for x=2/3 has been discussed by Hotta et 
al. |179j . As schematically shown in Figs. IT^ a) and (b), a remarkable feature can 
be observed. Due to the confinement of charge in the short straight segment for the 
bi-stripe phase, the charge stacking is suggested from our topological argument. On 
the other hand, in the Wigner- crystal type structure, charge is not stacked, but it is 
shifted by one lattice constant to avoid the Coulomb repulsion. Thus, if the charge 
stacking is also observed in the experiment for x=2/3, our topological scenario suggests 
the bi-stripe phase as the ground-state in the low temperature region. To establish the 
final "winner" in the competition between the bi-stripe and Wigner-crystal structure at 
x=2/3, more precise experiments, as well as quantitative calculations, will be further 
needed. 



3.3.5. x<0.5 Regarding densities smaller than 0.5, the states at x=l/8, 1/4 and 3/8 
have received considerable attention. See Refs. |18m I18H I182J . These investigations 
are still in a "fluid" state, and the experiments are not quite decisive yet, and for 
this reason, this issue will not be discussed in much detail here. However, without a 
doubt, it is very important to clarify the structure of charge-ordered states that may 
be in competition with the FM states in the range in which the latter is stable in 
some compounds. "Stripes" may emerge from this picture, as recently remarked in 
experiments |183[ I184[ I185| I186J and calculations [22] , and surely the identification of 
charge/orbital arrangements at x<0.5 will be an important area of investigations in the 
very near future. 

Here a typical result for this stripe-like charge ordering is shown in Fig. ^1 in 
which the lower-energy orbital at each site is depicted, and its size is in proportion to 
the electron density occupying that orbital. This pattern is theoretically obtained by 
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Figure 15. Numerical results for orbital densities in the FM phase for (a)x=l/2, 
(b)l/3, and (c)l/4 33 . The charge density in the lower-energy orbital is shown, and 
the size of the orbital is exactly in proportion to this density. 



the relaxation technique for the optimization of oxygen positions, namely including the 
cooperative Jahn- Teller effect. At least in the strong electron-phonon coupling region, 
the stripe charge ordering along the diagonal direction in the xy plane becomes the 
global ground-state. Note, however, that many meta-stable states can appear very close 
to this ground state. Thus, the shape of the stripe is considered to fluctuate both 
in space and time, and in experiments it may occur that only some fragments of this 
stripe can be detected. It should also be emphasized that the orbital ordering occurs 
concomitant with this stripe charge ordering. In the electron-rich region, the same 
antiferro orbital-order exists as that corresponding to x=0.0. On the other hand, the 
pattern around the diagonal array of electron-poor sites is quite similar to the building 
block of the charge/orbital structure at x=0.5. 

In Fig. Uni it is found that the same charge and orbital structure stacks along the 
6-axis. Namely, it is possible to cover the whole two-dimensional plane by some periodic 
charge-orbital array along the a-axis. If this periodic array is taken as the closed loop C, 
the winding numbers are w=l, 2, and 3, for x=l/2, 1/3, and 1/4, respectively. Note that 
in this case w is independent of the path along the a-axis. The results imply a general 
relation w={l — x)/x for the charge-orbital stripe in the FM phase, reflecting the fact 
that the distance between the diagonal arrays of holes changes with x. Our topological 
argument predicts stable charge-orbital stripes at special doping such as x=l/(l + w), 
with w an integer. 

This orbital ordering can be also interpreted as providing a "tt" -shift in the orbital 
sector, by analogy with the dynamical stripes found in cuprates |187j . although in copper 
oxides the charge/spin stripes mainly appear along the x- or ^/-directions. The study of 
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Figure 16. Schematic figure of the possible spin-charge-orbital structure at x=l/4 in 
the zigzag AF phase at low temperature and large electron-phonon coupling |33| . This 
figure was obtained using numerical techniques, and cooperative phonons, for Jh=oo 
and JAF=0.1io- For the non-cooperative phonons, basically the same pattern can be 
obtained. 



the similarities and differences between stripes in manganites and cuprates is one of tfie 
most interesting open problems in tlie study of transition metal oxides, and considerable 
work is expected in the near future. 

Finally, a new zigzag AF spin configuration for x<0.5 is here briefly discussed jSSj- 
In Fig. ^[ a schematic view of this novel spin-charge-orbital structure on the 8x8 
lattice at x=l/4 is shown, deduced using the numerical relaxation technique applied to 
cooperative Jahn- Teller phonons in the strong-coupling region. This structure appears 
to be the global ground state, but many excited states with different spin and charge 
structures are also found with small excitation energy, suggesting that the AF spin 
structure for x<0.5 in the layered manganites is easily disordered due to this "quasi- 
degeneracy" in the ground state. This result may be related to the "spin-glass" nature 
of the single layer manganites reported in experiments |188j . 

It should be noted that the charge-orbital structure is essentially the same as that 
in the two-dimensional FM phase, as shown in Fig. ^1 This suggests the following 
scenario for the layered manganites: When the temperature is decreased from the 
higher temperature region, first charge ordering occurs due to the cooperative Jahn- 
Teller distortions in the FM (or paramagnetic) region. If the temperature is further 
decreased, the zigzag AF spin arrangement is stabilized, adjusting itself to the orbital 
structure. Thus, the separation between the charge ordering temperature Tco and the 
Neel temperature T^ occurs naturally in this context. This is not surprising, since 
Tco is due to the electron- lattice coupling, while Tn originates in the coupling Jaf- 
However, if the electron-phonon coupling is weak, then Tco becomes very low. In this 
case, the transition to the zigzag AF phase may occur prior to the charge ordering. As 
discussed above, the Cg electron hopping is confined to one dimensional structures in the 



Orbital ordering phenomena in d- and f -electron systems 55 

zigzag AF environment. Thus, in this situation, even a weak couphng electron-phonon 
couphng can produce the charge-orbital ordering, as easily understood from the Peierls 
instability argument. Namely, just at the transition to the zigzag AF phase, the charge- 
orbital ordering occurs simultaneously, indicating that Tco=rN- Note also that in the 
zigzag AF phase, there is no essential difference in the charge-orbital structures for the 
non-cooperative and cooperative phonons, due to the one-dimensionality of those zigzag 
chains. 

3.4- Summary 

In this section, we have reviewed the theoretical results on spin, charge, and orbital 
ordering in manganites. We believe that the complicated ordering in manganites is 
caused by (i) competition between FM metallic and AF insulating phases and (ii) active 
Cg orbital degree of freedom. The existence of the FM metallic phase in the issue (i) has 
been understood by the double-exchange concept, while the variety of AF insulating 
originating from the point (ii) has not been considered satisfactorily in the standard 
double-exchange mechanism. Here we stress that the existence of active orbital does 
not simply indicate the increase of internal degrees of freedom in addition to spin and 
charge. We should remark an important effect of the orbital shape, leading to the 
geometrical pattern in the spin configuration. This point has been emphasized in this 
section in the context of topological aspect of orbital ordering. 

We have not mentioned another important characteristic issue of manganites, i.e., 
phase separation tendency, which is a driving force of colossal magneto-resistance 
phenomenon in manganites. The strong tendency of the phase separation is easily 
understood in the complex phase diagram including several kinds of first order transition. 
Readers should refer the previous review and textbook |Sl El, in which the phase- 
separation tendency and related physics have been explained in detail. 

4. Orbital physics in other (i-electron materials 

In the previous section, we have concentrated on the orbital physics of manganites. 
However, we can also observe orbital ordering phenomena in other transition metal 
oxides. Here we introduce possible orbital ordering in nickelates and ruthenates as 
typical materials of Cg and t2g electron systems, respectively, in the sense that Hg^ and 
Ht2 can be applied. Finally, we also discuss a potential role of orbital ordering in 
geometrically frustrated electron systems with orbital degeneracy. 

4-1- Cg electron systems 

The existence and origin of "striped" structures continues attracting considerable 
attention in the research field of transition metal oxides |1891 I19UJ . In a system with 
dominant electron-electron repulsion, the Wigner-crystal state should be stabilized, 
but in real materials, more complicated non-uniform charge structures have been 
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found. In Nd-based lightly-doped cuprates, neutron scattering experiments revealed 
incommensurate spin structures |19H 11921 I193J . where AF spin stripes are periodically 
separated by domain walls of holes. In La2-xSrxCu04, dynamical stripes are believed 
to exist along vertical or horizontal directions (Cu-0 bond direction). Note that for 
x<0.055, the spin-glass phase exhibits a diagonal spin modulation [1941 I195J . In 
nickelates, the charge-ordered stripes are along the diagonal direction 1961 11971 1198| 
I199J . In manganites, as mentioned in the previous section, evidence for striped charge- 
ordering also along the diagonal direction has been reported in the AF phase for x>l/2 
|174[ I175J . while short-range diagonal stripe correlations have been found in the FM 
phase at x<l/2 fn^ITT^ . 

In general, stripes can be classified into metallic or insulating. In La2_:rSr3;Cu04, 
the dynamical stripes exhibit metallic properties, but they are easily pinned by lattice 
effects and impurities. In Lai.6-3;Ndo.4Sr3;Cu04, stripes along the bond-direction are 
pinned by lattice distortions |189j . but they are still metallic. Intuitively, vertical or 
horizontal stripes could be associated with the formation of "rivers of holes" , to prevent 
individual charges from fighting against the AF background 120011187112011 . Such stripes 
should be metallic, even if they are pinned, since they are induced by the optimization 
of hole motion between nearest-neighbor Cu-sites via oxygens. 

However, in the diagonal stripes observed in manganites and nickelates, charges 
are basically localized, indicating that such insulating stripes are not determined just 
by the optimization of the hole motion. In the FM state of manganites, the hole 
movement is already optimal and, naively, charges should not form stripes. Obviously, 
an additional effective local potential must be acting to confine electrons into stripes. If 
such a potential originates in lattice distortions, it is expected to occur along the bond 
direction to avoid energy loss due to the conflict between neighboring lattice distortions 
sharing the same oxygens. Then, static stripes stabilized by lattice distortions tends to 
occur along the diagonal direction, as shown Fig. ^l which are stabilized by Jahn- Teller 
distortions [HH] . 

In simple terms, vertical or horizontal stripes in cuprates can be understood by 
the competition between Coulomb interaction and hole motion, while diagonal stripes 
are better explained as a consequence of a robust electron-lattice coupling. However, a 
difficulty has been found for theoretical studies of stripe formation in doped nickelates, 
since both Coulomb interaction and electron-lattice coupling appear to be important. 
Since the Ni^+ ion has two electrons in the Cg orbitals, on-site Coulomb interactions 
certainly play a crucial role to form spins S=l. When holes are doped, one electron 
is removed and another remains in the Cg orbitals, indicating that the hole-doped 
site should become JT active. Then, in hole-doped nickelates both Coulombic and 
phononic interactions could be of relevance, a fact not considered in previous theoretical 
investigations. 

In the following, we will review the recent results by Hotta and Dagotto |202j . 
The model for nickelates is the Cg orbital Hubbard Hamiltonian Eq. (jHSI), but another 
important ingredient is added here. Namely, the electron-lattice term is divided into 
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Figure 17. (a) Spin correlation 5(q) vs. Aforx=0. (b) Two kinds of local eg-electron 
arrangements for x=0. (c) AF spin pattern theoretically determined for A<3. 



couplings for the apical and in-plane oxygen motions. In layered nickelates, all NiOg 
octahedra are significantly elongated along the c-axis, splitting the Cg orbitals. This 
splitting from apical oxygens should be included explicitly from the start as the level 
splitting between a- and b-orbitals. Then, the model is defined as 

H = H^l + Hll^, + A J^K - nib)/2, (100) 

i 

where A is the level splitting. Later, the in-plane motion of oxygens should be studied 
by adding ifgp_ h- Note that the energy unit is also to in this subsection. 

First, consider the undoped case. The calculation is done for an 8-site tilted cluster, 
equivalent in complexity to a 16-site lattice for the single-band Hubbard model. Since at 
all sites the two orbitals are occupied due to the Hund's rule coupling, the JT distortions 
are not active and it is possible to grasp the essential ground-state properties using 
H . In Fig. ITTT a). the Fourier transform of spin correlations is shown vs. A, where 
5(g) = (l/iV)Ei,je^'^-('-J)(^fS'/), with S?= E-,(<T^i7T -4i^i7i)/2- As expected, a 
robust (vr, vr) peak can be observed for A<3, suggesting that the AF phase is stabilized 
by super-exchange interactions. The rapid decrease of S'(7r, vr) for A>3 is understood by 
comparing the energies for local triplet and singlet states, as shown in Fig. I17r b). The 
ground-state properties change at U'—J=U—A, leading to A=3J for the transition. 
The spin structure at x=0 is schematically shown in Fig. ITTT c) . 

Let us turn our attention to the case x=l/2. The 8-site tilted lattice is again used 
for the analysis, and the phase diagram Fig. IT^ a) is obtained for A=0.5. Since A 
of nickelates is half of that of cuprates from the lattice constants for CuOg and NiOg 
octahedra, it is reasonable to select A=0.5 in the unit of to- Increasing J, an interesting 
transformation from AF to FM phases is found. This is natural, since at large J the 
system has a formal similarity with manganite models, where kinetic-energy gains lead 
to ferromagnetism, while at small J the magnetic energy dominates. However, between 
the G-type AF for J^O and FM phase for Jp^U', unexpected states appear which are 
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Figure 18. (a) Ground-state phase diagram at x=l/2 witliout electron-phonon 
coupling, (b) S{q) for the CE- and E-type phases, at the couphngs indicated, (c) 
Spin and charge patterns for the CE- and E-type phases. These are schematic views, 
since focal charge-densities in practice are not exactly 1 and 2. 



mixtures of FM and AF phases, due to the competition between kinetic and magnetic 
energies. Typical spin correlations S{q) are shown in Fig. IT^ b). Note that peaks at 
g=(7r,0) and (7r/2,7r/2) indicate "C" and "E" type spin-structures, respectively. Double 
peaks at g=(7r,0) and (7r/2,7r/2) denote the CE-type structure, frequently observed in 
half-doped manganites j2U3j . In half-doped nickelates, the CE-phase is expressed as a 
mixture of type (I) and (II) in Fig. ITHT c). depending on the positions of the 5"=! and 
5*= 1/2 sites, although the "zigzag" FM chain structure is common for both types. The 
E-type phase is also depicted in Fig. ITHTc). Note that the charge correlation always 
exhibits a peak at q=(7r,7r) (not shown here), indicating the checkerboard-type charge 
ordering. 

In experimental results, a peak at (7r/2, 7r/2) in S{q) has been reported, suggesting 
an AF pair of S=l spins across the singly-occupied sites with holes. Moreover, the 
checkerboard-type charge ordering has been experimentally observed jl96| I197| I198| 
I199J . Thus, the spin-charge patterns of CE(II)- and E-type are consistent with the 
experimental results. Our phase diagram has a robust region with a peak at (7r/2, 7r/2), 
both for CE- and E-type phases, although the CE-phase exhibits an extra peak at (vr, 0). 
Whether the & or CE-phases are present in nickelates can be studied experimentally in 
the future by searching for this (vr, 0) peak. Note that if diffuse scattering experiments 
detect the AF correlation along the hole stripe, as has been found at x=l/3 |2()4j . 
the CE(II)-type may be the only possibility. Summarizing, the spin-charge structure 
in x=l/2 experiments can be understood within the Hamiltonian H by assuming a 



Orbital ordering phenomena in d- and f -electron systems 



59 




S=1/2 




^(b) x=1/2 U'=8 A-0.5 
J^6 X.=1.7 



0.0 0.1 0.2 0.3 0.4 0.5 
d (units of g/K) 




Figure 19. (a) Numerically obtained cooperative distortion pattern for an 8-site 
lattice at x=l/2 including H^[ ^^. Black and open circles indicate Ni and O ions, 
respectively. Open symbols indicate eg orbitals in the optimized state, (b) Total 
ground-state energy vs. d for x=l/2. (c) Orbital densities (r^i) and {txi) for sites 1-4. 
See (d) for the site labels. Optimized orbitals at d—Q.3 for sites 1 and 3 are also shown. 



relatively large J. 

Consider now the effect of in-plane oxygen motion. Note that apical oxygen motions 
have already been included as an Cg-level splitting. The extra electron-phonon coupling 
term is Hlf_^-^, Eq. (jl9|) . Since all oxygens are shared by adjacent NiOe octahedra, 
the distortions are not independent. To consider such cooperative effect, in principle, 
the 0-ion displacements should be optimized. However, in practice it is not feasible to 
perform both the Lanczos diagonalization and the optimization of all oxygen positions 
for 6- and 8-site clusters. In the actual calculations, Qu, Q2i, and Qsi are expressed by 
a single parameter d, for the shift of the 0-ion coordinate. Note that the unit of d is 
g/k, typically 0.1~0.3A. Then, the total energy is evaluated as a function of d to find 
the minimum energy state. Repeating these calculations for several distortion patterns, 
it is possible to deduce the optimal state. 

After several trials, the optimal distortion at x=l/2 is shown in Fig. ITW a). The 
diagonalization has been performed at several values of d on the 8-site distorted lattice 
and the minimum in the total energy is found at d=0.3, as shown in Fig. ITW b). As 
mentioned above, even without HJ^^, the checkerboard- type charge ordering has been 
obtained, but the peak at q={7i,7i) significantly grows due to the effect of lattice 
distortions. Note that the distortion pattern in Fig. ITW a) is essentially the same as that 
for half-doped manganites. This is quite natural, since JT active and inactive ions exist 
bipartitely also for half-doped nickelates. Then, due to this JT-type distortion, orbital 
ordering for half-doped nickelates is predicted, as schematically shown in Fig. fTT^ a). 
The shapes of orbitals are determined from the orbital densities, (r^i) and (r^.i), as 
shown in Fig. I19r c). The well-known alternate pattern of Sx^— r^ and 3y'^—r^ orbitals 
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Figure 20. (a) Zigzag 6-sites cluster covering the 2D lattice. Black circles denote Ni 
ions, and dashed lines indicate hole positions, (b) Phase diagram at x=l/3 without 
-^of-ph- Each phase is characterized by the momentum that shows a peak in S{q). (c) 
S{q) and (d) C{q) vs. J for f/'=6 and A=0.5. 



in half-doped manganites is denoted by dashed hnes. Increasing rf, the shape of orbitals 
deviates from 3a;^— r^ and 3j/^— r^, but it is still characterized by the orbitals elongating 
along the x- and y-directions. See insets of Fig. imf c). It would be very interesting to 
search for orbital ordering in half-doped nickelates, using the resonant X-ray scattering 
technique. 

Now let us move to the case x=l/3. If the actual expected stripe structure at x=l/3 
is faithfully considered |196l I197| I198t IT^ . it is necessary to analyze, at least, a 6x6 
cluster. However, such a large-size cluster with orbital degeneracy cannot be treated 
exactly due to the exponential growth of the Hilbert space with cluster size. Then, a 
covering of the two-dimensional (2D) lattice using zigzag 6-sites clusters as shown in 
Fig. EUT a) is considered by assuming a periodic structure along the diagonal direction. 

First we consider the case without H^^ ■^. The phase diagram obtained by analyzing 
the zigzag 6-site cluster for H is in Fig. I^HT b). Typical spin and charge correlations are 



shown in Figs. EDfc) and|2ni;d), where C(g) = (l/A^)^i . e^^^'-J) {{u; - {n))- 



n; 



in))), 



with ni=Y,^ni-y. 

Since the momentum q is defined along the zigzag direction in the unit of V^/a, 
where a is the lattice constant, the phase labeled by g=27r/3 in Fig. I^UT b) denotes an 
incommensurate AF phase with the proper spin stripe structure. The phase labeled by 
g=7r/3 indicates a spin spiral state, which will eventually turn to the FM phase in the 
thermodynamic limit. Thus, the spin stripe phase appears between the commensurate 
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Figure 21. (a) Cooperative distortion pattern for the zigzag 6-sites cluster at x=:l/3. 
(b) Total ground-state energy and (g) C{q) vs. d for x=l/3. 



AF and FM-like phases, similar to the case of x=l/2. However, as seen in Fig. I^UT d). 
C{q) in the spin stripe phase does not show the striped charge structure (g=27r/3). 
Rather, bipartite charge ordering characterized by a peak at q=TT still remains. Namely, 
the Hamiltonian without H^f^ , can explain the spin stripe, but does not reproduce the 
striped charge ordering at x=l/3, indicating the importance of H^p^. 

Consider now the effect of H^^ -^ for x=l/3. After evaluating total ground-state 
energies for several kinds of distortions, the pattern in Fig. l?1T a) has been found to 
provide the optimal state at x=l/3. This type of distortion induces a spatial modulation 
of the level splitting as — (5i/2=(52=53=— (^4/2=55 =^6, where 5i is the level splitting 
caused by the in-plane oxygen motions. Note that this breathing-mode modulation is 
consistent with experimental results j205j . The site numbers are found in Fig. l2]T a). 
The minimum energy is found at (i=0.1, as shown in Fig. 171T b). The modulation of 
level splitting stabilizes the striped charge ordering characterized by a g=27r/3 peak in 
C(g), as clearly shown in Fig. ITlT c). 

Note that (3x^— r^/3y^— r^)-type orbital ordering does not occur in Fig. l?1T a). 
Phenomenologically, such orbital ordering tends to appear in a hole pair separated by one 
site, the unit of the "bi-stripe" of manganites |174H175J . However, such a bi-stripe-type 
ordering contradicts the x=l/3 striped charge-ordering, and the bi-stripe-type solution 
was found to be unstable in these calculations. One may consider other distortion 
patterns which satisfy both (3a;^— r^/3|/^— r^)-type orbital and striped charge-ordering, 
but in such distortions no energy minimum was obtained for d>Q. After several trials. 
Fig- EHa) has provided the most optimal state. 

In summary, possible spin, charge, and orbital structures of layered nickelates have 
been discussed based on the Cg-orbital degenerate Hubbard model coupled with lattice 
distortions. To understand the nickelate stripes, both Hund's rule interaction and 
electron-lattice coupling appear essentially important. At x=l/2, (3x^— r^/3|/^— r^)- 
type orbital ordering similar to that in half-doped manganites is predicted. Even FM 
phases could be stabilized by chemically altering the carrier's bandwidth. For x=l/3, a 
spatial modulation in level splitting plays an important role for stripe formation. 
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4-2. t2g electron systems 

Let us now consider orbital ordering in a system with active t2g-orbital degree of 
freedom. As is well known, t2g-electron systems such as titanates and vanadates have 
been studied for a long time. After the discovery of superconductivity in layered 
cobalt oxyhydrate Nao.35CoO2-l.3H2O |2()(ij . magnetic properties of cobaltites have been 
discussed intensively both from experimental and theoretical sides. In relation with 
cobaltites, one-dimensional t2g electron model has been studied theoretically for the 
understanding of spin and orbital state of t2g electrons |2U7| I2U8J . When we turn our 
attention to Ad electron systems, Ru-oxides have been also focused, after the discovery 
of triplet superconductivity in the layered ruthenate Sr2Ru04 209J. In addition, the 
isostructural material Ca2Ru04 has been studied as a typical stage of spin and orbital 
ordering of t2g electrons. In this subsection, we review the orbital ordering phenomenon 
in Ca2Ru04. 

As mentioned in Sec. 2, Ru^+ ion include four electrons in the low-spin state, 
since the crystalline electric field is effectively larger than the Hund's rule interaction. 
Thus, four electrons occupy t2g orbitals, leading to S=l spin. The G-type AF phase in 
Ca2Ru04 is characterized as a standard Neel state with spin S=l |21()| 12111 I212J . The 
Neel temperature T-^ is 125K. To understand the Neel state observed in experiments, one 
may consider the effect of the tetragonal crystal field, leading to the splitting between 
xy and {yz,zx} orbitals, where the xy-orbital state is lower in energy than the other 
levels. When the xy-orbital is fully occupied, a simple superexchange interaction at 
strong Hund's rule coupling can stabilize the AF state. However, X-ray absorption 
spectroscopy studies have shown that 0.5 holes per site exist in the xy-orbital, while 
1.5 holes are contained in the zx- and yz-orbitals j2THj . suggesting that the above naive 
picture based on crystal field effects seems to be incomplete. This fact suggests that 
the orbital degree of freedom may play a more crucial role in the magnetic ordering in 
ruthenates than previously anticipated. 

First let us briefly review the result by Hotta and Dagotto |214j . The Hamiltonian 
is the t2g Hubbard model coupled with Jahn- Teller distortions, already given by Ht^ , 
Eq. ()63|) in Sec. 2. This model is believed to provide a starting point to study the 
electronic properties of ruthenates, but it is difficult to solve even approximately. To 
gain insight into this complex system, an unbiased technique should be employed first. 
Thus, Hotta and Dagotto have analyzed a small 2x2 plaquette cluster in detail by 
using the Lanczos algorithm for the exact diagonalization, and the relaxation technique 
to determine the oxygen positions. In actual calculations, at each step for the relaxation, 
the electronic portion of the Hamiltonian is exactly diagonalized for a fixed distortion. 
Iterations are repeated until the system converges to the global ground state. 

The ground state phase diagram obtained by Hotta and Dagotto is shown in Fig. 1221 
There are six phases in total, which are categorized into two groups. One group is 
composed of phases stemming from the U'=0 or Ejt=0 limits. The origin of these 
phases will be addressed later, but first their main characteristics are briefiy discussed. 
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Figure 22. Ground state phase diagram for tlie i2g Hubbard model coupled with 
Jahn- Teller phonons for J—iU' /A. The notation is explained in the maintext. 



For £'jT=0, a C-type AF orbital disordered (OD) phase appears in the region of small 
and intermediate U' . This state is characterized by n^y\nyz+n.^^=l/2\?)/2, where n^ is the 
hole number per site at the 7-orbital. Hereafter, a shorthand notation such as "1/2:3/2" 
is used to denote the hole configuration. For large f/', and still ii^jT=0, a FM/OD phase 
characterized by 3/4:5/4 is stable, which may correspond to Sr2Ru04. On the other 
hand, for [/'=0 and small £^jt, a "metallic" (M) phase with small lattice distortion 
is observed, while for large £^jt, a charge- density- wave (CDW) state characterized by 
1:1 was found. In short, the G-type AF phase observed experimentally |212j does not 
appear, neither for ii^jT=0 nor for f/'=0. 

Another group includes two phases which are not connected to either i?jT=0 or 
f/'=0. It is only in this group, with both lattice and Coulomb effects being relevant, 
that for intermediate U' the G-type AF and orbital ordered (00) phase with 1/2:3/2 
found in experiments |213j is stabilized. At larger [/', a FM/00 phase occurs with the 
same hole arrangement. In the FM phase, since an S=l spin with Sz=+1 is formed at 
each site, the up-spin number is unity at each orbital, while the down-spin distribution 
depends on the orbital. In the AF state, the configuration of double-occupied orbitals is 
the same as in the FM phase, but the single-occupied orbital contains 0.5 up- and 0.5- 
down spins on average, since the S=l spin direction fluctuates due to the AF coupling 
between neighboring S=l spins. However, the spin correlations peak at (vr,??), indicating 
the G-AF structure. Except for the spin direction, the charge and orbital configuration 
in the FM/00 phase is the same as in the G-AF/00 state. An antiferro-orbital ordering 
pattern including xy, yz, and zx orbitals has been suggested for these FM and AF phases. 

On the other hand, a ferro "0:2" xy-orbital ordered state has been suggested by 
Anisimov et al. |215j . Fang et al. also predicted the ferro- type orbital ordering j216j . It 
seems to be different from experiments on the hole distribution in Ref. |213j . but due to 
the combination of optical conductivity measurement and LDA+U calculations |217j , it 
has been found that xy-orbital ferro ordering occurs and the change of hole population 
can be explained due to the temperature dependence of electronic structure. 
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Figure 23. (a) Azimuthal angle dependence of the interference term for a main edge 
peak at 305 and 6 K at Q={0,2,6). The thick and thin curves denote the analysis 
results at 305 and 6 K, respectively, (b) Temperature dependence of the interference 
term at Q={0,2,6). As for details, readers refer Ref. |218j 



Recently, Kubota et al. have performed the experiment to determine the orbital 
ordering in Ca2Ru04 by using the resonant X-ray scattering interference technique at 
the K edge of Ru |218j . In this new and skillful technique, it is remarkable that the 
dxy orbital ordering is observed even at room temperature, in which the Jahn- Teller 
distortion is negligible. Note here that the Jahn- Teller distortion is defined as the 
ratio of the apical Ru-0 bond length to the equatorial Ru-0 bond length in the RuOe 
octahedron. 

The resonant x-ray scattering (RXS) measurement has been very powerful method 
to detect the orbital ordering, but the conventional RXS measurement is not useful for 
the observation of a ferro-type orbital state, since it is difficult to extract the signal for 
the ferro-orbital ordered state at F point in a momentum space, which is accompanied 
with a large amplitude of a fundamental reflection by Thomson scattering. However, the 
RXS interference technique can observe the ferro-type orbital ordered state, in which 
the signal is magnified by the interference with a fundamental signal. 

In Figs. I2S1 we show typical results at the K edge of Ru, obtained by Kubota et al. 
FigureEB^a) denotes the azimuthal angle dependence of the interference term for a main 
edge peak. We note that the signal exhibits the characteristic oscillation with the period 
of 360 degrees. Moreover, the significant signal can be found even at 305 K, in which 
the Jahn- Teller distortion is negligible. This fact suggests that the interference term is 
directly related to the orbital ordering. In Fig. l23r b). the temperature dependence of the 
RXS signal is shown. Below 200K, we can observe that the RXS signal for the ferro-type 
xy-orbital ordering is almost saturated due to the occurrence of the G-AF Neel state. 
Above 200 K, the magnitude of the signal is gradually decreased and becomes zero at the 
metal-insulator transition temperature (~357 K). Since the apical bond length of RuOg 
becomes almost equal to the averaged equatorial bond length around at 300K |212j . it is 
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difficult to consider that the Jahn- Teller distortion is the origin of the orbital ordering. 
Thus, the coupling of t2g electron with Jahn- Teller distortion is not the primary term 
of the electron-phonon coupling part. Rather, the tilting and/or buckling modes should 
be included seriously in the model Hamiltonian. 

Finally, let us briefly mention another result of a resonant X-ray diffraction study 
on Ca2Ru04 at the Ru Lu and Lm edges |219j . Zegkinoglou et al. have observed a 
significant enhancement of the magnetic scattering intensity at the wave vector which 
characterizes the AF ordering. Then, they have found a phase transition between two 
paramagnetic phases around 260 K, in addition to the well-known AF transition at 
Tn=110K. Due to the analysis of polarization and azimuthal angle dependence of the 
diffraction signal, Zegkinoglou et al. have concluded that the transition at 260K is 
attributed to the orbital ordering of Ru t2g electrons. This orbital order is characterized 
by the same propagation vector as the low-temperature AF phase. Note, however, that 
the ferro-orbital component of the ordering pattern cannot be ruled out, as mentioned 
by Zegkinoglou et al. 

4-3. Geometrically frustrated systems 

As an important ingredient to understand novel magnetism of actual strongly 
correlated electron materials, thus far we have emphasized a potential role of 
orbital degree of freedom, when electrons partially ffil degenerate orbitals. However, 
on the lattice with geometrical frustration, a subtle balance among competing 
interactions easily leads to a variety of interesting phenomena such as unconventional 
superconductivity and exotic magnetism. The recent discovery of superconductivity 
in layered cobalt oxyhydrate Nao.35CoO2-l.3H2O j2U6j has certainly triggered intensive 
investigations of superconductivity on the triangular lattice. Concerning the magnetism, 
antiferromagnetism on the triangle-based structure has a long history of investigation 
|22Uj . In the low- dimensional system, the combined effect of geometrical frustration and 
strong quantum fluctuation is a source of peculiar behavior in low-energy physics, as 
typically found in the Heisenberg zigzag chain with spin 5*= 1/2. As the strength of 
frustration is increased, the ground state is known to be changed from a critical spin- 
liquid to a gapped dimer phase |22H 1222^ \'2'2'A\ I224J . In the dimer phase, neighboring 
spins form a valence bond to gain the local magnetic energy, while the correlation among 
the valence bonds is weakened to suppress the effect of spin frustration. 

Here we have a naive question: What happens in a system with both active orbital 
degree of freedom and geometrical frustration? It is considered to be an intriguing 
issue to clarify the influence of orbital ordering on magnetic properties in geometrically 
frustrated systems. For instance, significant role of t2g-orbital degree of freedom has 
been remarked to understand the mechanism of two phase transitions in spinel vanadium 
oxides AV2O4 (A=Zn, Mg, and Cd) [223 1223 1223 • I* has been proposed that orbital 
ordering brings a spatial modulation in the spin exchange and spin frustration is 
consequently relaxed. Similarly, for MgTi204, the formation of a valence-bond crystal 
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Figure 24. Lattice location and site numbering of A^-site ladder and zigzag chain. 
The length is defined as L=N/2. 



due to orbital ordering has been also suggested |228, 229J. 

Since (i- and /-electron orbitals are spatially anisotropic, there always exist easy 
and hard directions for electron motion. Thus, it is reasonable to expect that the 
effect of geometrical frustration would be reduced due to orbital ordering, depending 
on the lattice structure and the type of orbital, in order to arrive at the spin structure 
which minimizes the influence of frustration. However, the spin structure on such an 
orbital-ordered background may be fragile, since the effect of geometrical frustration 
never vanishes, unless the lattice distortion is explicitly taken into account. It is a 
highly non-trivial problem, whether such an orbital arrangement actually describes the 
low-energy physics of geometrically frustrated systems. In particular, it is important 
to clarify how the orbital-arranged background is intrinsically stabilized through the 
spin-orbital correlation even without the electron-lattice coupling. In this subsection, 
we review the recent result by Onishi and Hotta concerning the role of orbital ordering 
in the geometrically frustrated lattice [SHHl EHD EHS EHS] • 

Onishi and Hotta have considered an Cg-orbital model on the A^-site ladder or zigzag 
chain, including one electron per site with two orbitals, i.e., quarter filling. The lattice 
they have used is shown in Fig. |^ Note that the zigzag chain is composed of equilateral 
triangles. The Cg-orbital degenerate Hubbard model is already given by Eq. (jUHj) . but 
the electron-phonon term is not considered. Namely, the model is written as 



He, = ^kfn + HZ 



-el- 



(101) 



Here the (i-electron hopping amplitude t^ , for the oblique u direction is defined by 
t^^=ti/4, t^^=t^^=y/3ti/8, tS,=3ti/16. Note the relation of t'^~,''=t'^y. Concerning 
hopping amplitudes along x- and ^/-directions, see Eqs. fjS^ and (|3B| . In this subsection, 
ti is taken as the energy unit. 

In order to analyze the complex model including both orbital degree of freedom 
and geometrical frustration, Onishi and Hotta have employed the finite-system density 
matrix renormalization group (DMRG) method, which is appropriate for the analysis of 
quasi-one-dimensional systems with the open boundary condition |234U235ll236j . Since 
one site includes two Cg orbitals and the number of bases is 16 per site, the size of the 
superblock Hilbert space becomes very large as m'^ x 16^, where m is the number of states 
kept for each block. To accelerate the calculation and to save memory resources, Onishi 
has skillfully reduced the size of the superblock Hilbert space to m^x4^, by treating 
each orbital as an effective site. In the actual calculations, m states up to m=200 were 
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Figure 25. The spin-correlation function measured from the center of the lower 
chain for the PM ground state in (a) the ladder and (b) the zigzag chain. The local 
magnetization for the first spin-excited state in (c) the ladder and (d) the zigzag chain. 



kept in the renormalization process and the truncation error was estimated to be 10~^ 
at most. 

Now we introduce the results on the spin structure of the paramagnetic (PM) 
ground state at J=0, since the zigzag chain is relevant to a geometrically frustrated 
antiferromagnet in the spin-singlet PM phase. Here we refer only the result for J=0, 
but readers should consult with Ref. |23()j about the results for J^O. In Figs. OMa) and 
l2Sfb), we show the DMRG results for A^=40 of the spin-correlation function Cspin(i, j)= 
{S?S^) with 5*? =^ (pi^f— pi^i)/2. Note that a large value of U'=20 was used to consider 
the strong-coupling region, but the results did not change qualitatively for smaller values 
of U'. As shown in Fig. l25f a). we observe a simple Neel structure in the ladder. On the 
other hand, in the zigzag chain, there exists AF correlation between intra-chain sites in 
each of lower and upper chains, while the spin correlation between inter-chain sites is 
much weak. Namely, the zigzag chain is considered to be decoupled to a double chain 
in terms of the spin structure. 

In order to clarify the characteristics of the spin structure in the excited state, 
Onishi and Hotta have investigated the local magnetization Mi=(S'f) for the lowest- 
energy state with 5'^^^^=!, i.e., the first spin-excited state, where S^^^ is the z component 
of the total spin. In the ladder, the total moment of S^^^=l is distributed to the whole 
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Figure 26. Optimal orbital arrangement in (a) the ladder and (b) the zigzag chain. 

system and there is no significant structure, as shown in Fig. EH c). On the other hand, 
the situation is drastically changed in the zigzag chain. As shown in Fig.|23fd), the total 
moment of S^^^=l is confined in the lower chain and it forms a sinusoidal shape with 
a node, while nothing is found in the upper chain. Note that the sinusoidal shape of 
the local magnetization is characteristic of the 5*= 1/2 AF Heisenberg chain with edges 
at low temperatures |237| 1238] . Thus, the double-chain nature in the spin structure 
remains robust even for the spin-excited state. 

Onishi and Hotta have also discussed the orbital arrangement to understand the 
mechanism of the appearance of the spin structures. For the determination of the orbital 
arrangement, orbital correlations are usually measured, but due care should be paid to 
the definition. By analogy with Eqs. (f77j) and (ff^ which have treated the phase of the 
JT distortions, phase-dressed operators are introduced as 

diaa = e''>-^/^[cos{ei/2)diaa + sm{ei/2)dib^], 

diba = e'''l\- sm{e-j2)d,a„ + cos{e-j2)d,,,]. ^ ' 

Then, the optimal set of {^i} is determined so as to maximize the orbital-correlation 
function, which is defined as 

T(q) = (l/iV2) ^{f?f^)e''^-^'-'^\ (103) 

With f.: = Y.M\aAa.-diJ,l,.)l2. 

As shown in Fig. I^UT a). in the case of the ladder, Onishi and Hotta have found 
that a ferro-orbital (FO) ordering, characterized by ^i/7r~1.18, appears in the ground 
state. In the first spin-excited state, the FO structure also appears, but the angle 
characterizing the orbital shape is slightly changed as ^i/7r~1.20 to further extend to 
the leg direction. On the other hand, in the zigzag chain, it is observed that hath in the 
ground and first spin-excited states, T(q) becomes maximum at q=0 with 6'i/7r~1.32, 
indicating a 3x^— r^ orbital at each site, as shown in Fig. |2El^b). Note that the orbital 
arrangement is unchanged even in the spin-excited state. Namely, the orbital degree of 
freedom spontaneously becomes "dead" in low-energy states to suppress the effect of 
spin frustration. 

It is interesting to remark that the spin-exchange interactions become anisotropic 
due to the orbital-arranged background. First let us consider the zigzag chain, which is 
effectively described by the Hubbard model composed of 3a;^— r^ orbital. It is intuitively 
understood that the AF exchange interaction along the u direction Ji should be much 
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weaker than that along the x direction J2, since the orbital shape extends along the 
double chain, not along the zigzag path, as shown in Fig. I26r b). In order to estimate 
the ratio of J1/J2, it is enough to consider the hopping amplitudes between adjacent 
optimal "a"-orbitals, 3x^— r^ in this case, which are given by t^^=l and t"^=l/64. Then, 
taking account of the second-order process in terms of electron hopping between only 
3^2_^2 orbitals, we obtain J1/J2 =[2{i^^y /U]/[2{i^J'^/U] =1/642. This small value of 
J1/J2 clearly indicates that the spin correlation on the zigzag path is reduced due to 
the spatial anisotropy of Sx^— r^ orbital. Thus, the zigzag chain is effectively reduced 
to a double-chain system of the 5"= 1/2 AF Heisenberg chain, suggesting that the spin 
gap should be extremely suppressed, since the spin gap decreases exponentially with the 
increase of J2/J1 in the gapped dimer phase in the zigzag spin chain |224j . 

On the other hand, in the ladder with the ferro-orbital structure as shown in 
Fig. I^M a). the orbital shape extends to the rung direction as well as to the leg 
direction. When we define Jieg and Jrung as the AF exchange interactions along the 
leg and rung directions, respectively, we obtain Jrung/</ieg=0.26, which is much larger 
than Ji/J2=l/642 in the case of the zigzag chain. The spin correlation on the rung 
is considered to remain finite, leading to the simple Neel structure. Thus, the spin 
excitation in the ladder is expected to be gapful similar to the spin ladder |239tl24UJ . 

We have reviewed both ground- and excited-state properties of the Cg-orbital 
degenerate Hubbard model on the ladder and the zigzag chain. It has been found that 
the zigzag chain is reduced to a decoupled double-chain spin system due to the selection 
of a specific orbital. It is considered as a general feature of geometrically frustrated 
multi-orbital systems that the orbital selection spontaneously occurs so as to suppress 
the effect of spin frustration. 

Finally, let us briefly comment on the effect of level splitting between x^—y'^ 
and 3^2— r^ orbitals, which has not been considered in the present Hamiltonian. In 
particular, when "iz^—r"^ orbital is the lower level which is well separated from x^—y^ 
orbital, the hopping amplitude does not depend on the direction and the effect of spin 
frustration revives for the system with isotropic AF interactions. In such a region with 
strong spin frustration, a finite energy gap between ground and first-excited states can 
be clearly observed. Naively thinking, it may be called a spin gap, but we should note 
that the orbital arrangement is significantly influenced by the spin excitation. In general, 
the energy gap between ground and first-excited states in multi-orbital systems should 
be called a spin-orbital gap. As for details, readers consult with Ref. |241j 

5. Model Hamiltonian for /-electron systems 

Thus far, we have reviewed the theoretical results on orbital ordering phenomena of 
(i-electron systems. As typical examples, we have picked up manganites, nickelates, 
and ruthenates. However, there exists another spin-charge-orbital complex system such 
as /-electron compounds. In the latter half of this article, we review orbital ordering 
phenomena of /-electron systems. Before proceeding to the description of the theoretical 
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results, again it is necessary to set the model Hamiltonian for /-electron systems. In 
order to construct such a microscopic Hamiltonian, we must include simultaneously the 
itinerant nature of / electrons as well as the effects of strong electron correlation, CEF, 
and spin-orbit interaction. Among them, the existence of strong spin-orbit interaction is 
essential difference from d-electron systems. The inclusion of the spin-orbit interaction 
is a key issue, when we construct the model Hamiltonian for / electron materials. Since 
this is a complicated problem, it is instructive to start the discussion with a more basic 
level. Namely, we first review in detail the single ion problem focusing on the properties 
of local /-electron states in comparison with those obtained in the LS and j-j coupling 
schemes. Then, we move on to the explanation of the microscopic /-electron model on 
the basis of the j-j coupling scheme. 

5.1. LS vs. j-j coupling schemes 

In the standard textbook, it is frequently stated that for rare-earth ion systems, the 
LS coupling scheme works well, while for actinides, in particular, heavy actinides, the 
j-j coupling scheme becomes better. However, do we simply accept such a statement? 
Depending on the level of the problem in the condensed matter physics, the validity 
of the approximation should be changed, but such a point has not been explained in 
the textbook. It is important to clarify which picture is appropriate for the purpose to 
consider the many-body phenomena in /-electron systems. 

Let us generally consider the /" configuration, where n is the number of / electrons 
included on a localized ion. In the LS coupling scheme, first the spin S and angular 
momentum L are formed due to Hund's rules as S=J2^^i *« ^'^^ L=Y^'^=if-i, where s, 
and £i are spin and angular momentum for i-th f electron. Note that the Hund's rules 
are based on the Pauli principle and Coulomb interactions among / electrons. After 
forming S and L, we include the effect of spin-orbit interaction, given by XL-S, where A 
is the spin-orbit coupling. We note that A>0 for n<7, while A<0 for n>7. Note also that 
a good quantum number to label such a state is the total angular momentum J, given 
by J=L+S. Following from simple algebra, the ground-state level is characterized by 
J=\L—S\ for n<7, while J=L-\-S for n>7. 

On the other hand, when the spin-orbit interaction becomes larger than the 
Coulomb interactions, it is useful to consider the problem in the j-j coupling scheme. 
First, we include the spin-orbit coupling so as to define the state labeled by the 
total angular momentum jfj for the i-th electron, given by ji=Si+£i. For /-orbitals 
with i=3, we immediately obtain an octet with j=7/2(=3+l/2) and a sextet with 
j=5/2(=3— 1/2), which are well separated by the spin-orbit interaction. Note here that 
the level for the octet is higher than that of the sextet. Then, we take into account the 
effect of Coulomb interactions to accommodate n electrons among the sextet and/or 
octet, leading to the ground-state level in the j-j coupling scheme. 

As is easily understood from the above discussion, the LS coupling scheme works 
well under the assumption that the Hund's rule coupling is much larger than the spin- 
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orbit interaction, since S and L are formed by the Hund's rule coupling prior to 
the inclusion of spin-orbit interaction. It is considered that this assumption is valid 
for insulating compounds with localized / electrons. However, when the spin-orbit 
interaction is not small compared with the Hund's rule coupling, the above assumption 
is not always satisfied. In addition, if the / electrons become itinerant owing to 
hybridization with the conduction electrons, the effect of Coulomb interactions would 
thereby be effectively reduced. In rough estimation, the effective size of the Coulomb 
interaction may be as large as the bandwidth of / electrons, leading to a violation of 
the assumption required for the LS coupling scheme. 

Furthermore, even in the insulating state, we often encounter some difficulties to 
understand the complex magnetic phases of /-electron systems with active multipole 
degrees of freedom from a microscopic viewpoint. In a phenomenological level, it is 
possible to analyze a model for relevant multipoles obtained from the LS coupling 
scheme, in order to explain the phenomena of multipole ordering. However, it is difficult 
to understand the origin of the interaction between multipoles in the LS coupling 
scheme. 

From these viewpoints, it seems to be rather useful to exploit the j-j coupling 
scheme for the purpose to understand magnetism and superconductivity of /-electron 
materials. Since individual /-electron states is clearly defined, it is convenient for 
including many-body effects using the standard quantum-field theoretical techniques. 
However, it is not the reason to validate to use the j-j coupling scheme for the model 
construction. In order to clarify how the j-j coupling scheme works, it is necessary to 
step back to the understanding of the local /-electron state. In the next subsection, let 
us consider this issue in detail. 

5.2. Local f- electron state 

In general, the local /-electron term is composed of three parts as 

Hf = /7el-cl + ^so + i^CEF, (104) 

where Hq is the Coulomb interaction term, written as 

-nel-el — / ^ / ^ / ^ -'mi,m2,m3,m4/imio-i/im20-2-'i™30'2/im4o-i- l-'-UOJ 

i mi~m4 f7i,cr2 

Here fima is the annihilation operator for /-electron with spin a and angular momentum 
m(=— 3,- • -,3) at a site i. Similar to the li-electron case, the Coulomb integral Ilii,m2,m3,m4 
is given by 

6 

lLuni2,ni3,mi = S„,,+„,2^m,+^^'^F}c'-''\mi,m4)c^''\m2,m3), (106) 

fc=0 

where the sum on k includes only even values {k=0, 2, 4, and 6), F^ is the Slater-Condon 
parameter for / electrons including the complex integral of the radial function, and c^ 
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is the Gaunt coefficient. It is convenient to express the Slater-Condon parameters as 

F° = A + 15C + 9D/7, 

F| = 225(5 -6C/7 + D/42), 

F/ = 1089(5C/7 + 1^/77), ^ ' 

Ff = (429/5)2 ■(D/462), 

where A, B, C, and D are the Racah parameters for / electrons \V2'A\ . 
The spin-orbit coupling term, Hgo, is given by 

-^so = _^ ^^ KoCm,cr,m' ,cr' fimafim'a' , (108) 



m,a,m' ,(T' 



where Ago is the spin-orbit interaction and the matrix elements are explicitly given by 

Cm+u,m^ = v /l2-m(m + l) /2, (109) 

Cm-i,T,m,i = a/12 -m(m- l)/2, 
and zero for other cases. 

The CEF term Hcef is given by 

-f^CEF = 2^ Am,Tn'fymafim'a, (HO) 

i,m,m' ,(7 

where Am,m' can be evaluated in the same manner as has done in Sec. 2 for d electrons 
with i=2. However, there is no new information, if we repeat here lengthy calculations 
for / electrons with i=3. As already mentioned in Sec. 2, it is rather useful and 
convenient to consult with the table of Hutchings for angular momentum J=3 |12nj . 
For cubic symmetry, Am,m' is expressed by using a couple of CEF parameters, B^ and 
Bl as 

^3,3 = ^-3,-3 = 1805° + 1805°, 

A2,2 = ^-2,-2 = -4205° - 10805°, 

^1,1 = ^-1,-1 = 605° + 27005°, 

Ao,o = 3605° - 36005°, ^ ^ 

^3,-1 = ^-3,1 = 60^15(5° - 215°), 

^2,-2 = ^-2,2 = 3005° + 75605°. 

Following the traditional notation, we define 

5° = W^x/F(4), 

5° = Wil - |a;|)/F(6), ^ ' 

where x specifies the CEF scheme for 0^ point group, while W determines an energy 
scale for the CEF potential. Although F{4) and F{Q) have not been determined 
uniquely, we simply follow the traditional definitions as -F(4)=15 and F(6)=180 for 

j=3 [nni. 

Here we note that the CEF potential is originally given by the sum of electrostatic 
energy from the ligand ions at the position of /-electron ion, leading to the one-electron 
potential acting on the charge distribution of /-orbitals, as expressed by Eq. ()11()|1 . 
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Figure 27. Energies of / electrons as functions of x for (a) the LS coupling and (b) 
the j-j coupling schemes for n=2. The magnitude of the CEF potential energy is fixed 
as \W\/U=0.001. 



Thus, in principle, it is not necessary to change the CEF potential, depending on the 
/-electron number. As we will see later, the CEF schemes for ?t,=1~13 are automatically 
reproduced by diagonalizing the local /-electron term H\oc, once we fix the CEF 
parameters in the form of one-electron potential Eq. ()11()|1 . 

Now we compare the electronic states of Hi^c with those of LS and j-j coupling 
schemes. We believe that it is quite instructive to understand the meanings of the CEF 
potential in /-electron systems. We introduce "?7" as an energy scale for the Racah 
parameters. A, B, C, and D. In this subsection, U is the energy unit, which is typically 
considered to be 1 eV. In /-electron compounds, the magnitude of the CEF potential 
is much smaller than both spin-orbit coupling and Coulomb interactions. Thus, it is 
reasonable to consider that W is always much smaller than Ago and U. However, there 
occur two situations, depending on the order for taking the limits of Xgo/W^oo and 
U/W^oo. When the limit of U/W^oo is first taken and then, we include the effect 
of the spin-orbit coupling Ago, we arrive at the LS coupling scheme. On the other 
hand, it is also possible to take first the infinite limit of Xso/W. After that, we include 
the effect of Coulomb interaction, leading to the j-j coupling scheme. In the present 
local /-electron term Hioc, it is easy to consider two typical situations for /-electron 
problems, \W\<t^\so<U and \W\'^U <\so, corresponding to the LS and j-j coupling 
schemes, respectively. 

Let us consider the case of n=2 as a typical example of the comparison between 
the two schemes. In the LS coupling scheme for the /^-electron system, we obtain 
the ground-state level as ^H with 5=1 and L=5 from the Hund's rules, where S and 
L denote sums of /-electron spin and angular momentum, respectively. Upon further 
including the spin-orbit interaction, the ground state is specified by J=4 expressed as 
^H^ in the traditional notation. Note that the total angular momentum J is given by 
J=\L — S\ and J=L + S for n<7 and n>7, respectively. In order to consider further 
the CEF effect, we consult with the table of Hutchings for the case of J=4. In the LS 
coupling scheme, W is taken as W/U=0.001 and we set F(4)=60 and F(6)=1260 for 
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J=4 by following the traditional definitions |12(Jj . Then, we easily obtain the nine eigen 
values, including Fi singlet, T^ doublet, and two kinds of triplets, r4 and T^, as shown 
in Fig. n7\ a). Note that for odd n, the eigenstate has odd parity, specified by "u" in 
MuUiken's notation and "— " in Bethe's notation, while the even n configuration has 
even parity, labeled by "g" and "+" |242j . When we use Bethe's notation to specify the 
/-electron eigenstate, the "+" or "— " superscript is suppressed for convenience. 

In the j-j coupling scheme, on the other hand, first we take the infinite limit of Ago- 
Thus, we consider only the j=5/2 sextet, where j denotes the total angular momentum 
of one / electron. In the /^-electron system, two electrons are accommodated in the 
sextet, leading to fifteen eigen states including J=4 nontet, J=2 quintet, and J=0 
singlet. Due to the effect of Hund's rule coupling, J=4 nontet should be the ground 
state. When we further include the CEF potential, it is necessary to reconsider the 
accommodations of two electrons in the /^-electron potential with F7 doublet and Fg 
quartet. Thus, in the j-j coupling schemes, except for the energy scale W, only relevant 
CEF parameter is x, leading to the level splitting between F7 doublet and Fg quartet. 
For the j-j coupling scheme, we set -F(4)=60 and W/U=— 0.001. Note that the minus 
sign in W is added for the purpose of the comparison with the LS coupling scheme. As 
shown in Fig. OTTb). the J=4 nontet is split into Fi singlet, F3 doublet, F4 triplet, and 
F5 triplet. The ground state for x>0 is F5 triplet composed of a couple of Fg electrons, 
while for x<0, it is Fi singlet which is mainly composed of two F7 electrons. Note that 
for x>0, the first excited state is F4 triplet, composed of F7 and Fg electrons. 

At the first glance, the energy levels in the j-j coupling scheme seems to be different 
from those of the LS coupling scheme. How do we connect these different results? In 
order to answer to this question, let us directly diagonalize H{ by changing U and A. 
Here it is convenient to introduce a new parameter to connect the LS and j-j coupling 
schemes as 

/, ^-/\w\ ni3) 

u/\w\ + x,j\wy ^''^^ 

where we explicitly show |iy| in this formula, since both U and Ago should be always 
very large compared with \W\ in actual /-electron compounds. Note that k=0 and 
1 are corresponding to the limits of Xso/U=0 and Aso/^=oo, respectively. Then, we 
can control the change of two schemes by one parameter k, by keeping ?7/|iy|^l and 

Aso/|W^|»l. 

In Figs. l2Hta)-(d), we show the energy levels of Hf for several values of k with both 
Aso and U larger than \W\. Racah parameters are set as A/U=10, B/U=0.3, C/U=0.1, 
and D/U=0.05 in the units of U. As described above, the CEF potential is always small 
and here we set W/U=— 0.001. In Fig. EHT a). results for k=0.1 are shown. In this case, 
Aso/f^=0.11, while the condition Aso/|W^|^l is still satisfied. Without the spin-orbit 
interaction, the ground-state level is expressed as ^H with S=l and L=5 due to the 
Hund's rules. When we increase Ago, the multiplets labeled by J are well separated and 
the ground-state level is specified by J=4, as expected from the LS coupling scheme. 
Then, the energy levels in Fig. OHT a) are quite similar to those of Fig. ITfT a). since we 
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Figure 28. Eigen energies of H{ as functions of x for (a) fc=0.1, (b) fc=0.8, and (c) 
fc=0.99. Racah parameters are set as A/U=10, B/U=0.3, C/U=0.1, and D/U=Om. 
The energy scale for CEF potentials are given by W/U=— 0.001. (d) Overlap integral 
between the eigenstate of H{ and that in the LS coupling scheme for the case of the 
Fi ground state. Solid squares at fc=:0 and 1 are obtained separately from the LS and 
j-j coupling schemes, respectively. 



are now in the region where the LS couphng scheme is appropriate. 

Even when Ago is further increased and k is equal to 0.5, the structure of the 
energy levels is almost the same as that of the LS coupling scheme (not shown here). 
However, when k becomes 0.8, as shown in Fig. EHl^b), the energy level structure is 
found to be deviated from that of the LS coupling scheme. Rather, it becomes similar 
to the energy level structure of the j-j coupling scheme. To see the agreement with 
the j-j coupling scheme more clearly, we consider very large A which gives k=0.99. As 
shown in Fig. |2Hfc), we can observe the energy level structure similar to Fig. ITfT b). In 
particular, the region of the Fs ground state becomes very narrow, as discussed later. 
Thus, it is concluded that Hf actually reproduces the energy levels both for the LS and 
j-j coupling schemes. We also stress that H^ provides correct results in any value of 
/-electron number. 

A crucial point is that the structure of energy levels is continuously changed, as 
long as Ago and U are large compared with the CEF potential. Namely, the states 
both in the LS and j-j coupling schemes are continuously connected in the parameter 
space. Thus, depending on the situation to consider the problem, we are allowed to use 
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either LS or j-j coupling scheme. In order to clarify this point, we evaluate the overlap 
(\E'|\E'ls), where |\E') and I^E'ls) are the eigenstate of iff and that in the LS coupling 
scheme, respectively. In Fig. OHT d). we show the overlap for the case of Fi ground state 
for x=-0.1 and W/U=-Omi. For k=0, (^|^ls)=1 due to the definition. The overlap 
is gradually decreased with the increase of k, but it smoothly converges to the value 
at fc=l, i.e., the j-j coupling scheme. Note that the overlap between the eigenstates 
of the LS and j-j coupling schemes is as large as 0.865, which seems to be larger than 
readers may naively anticipated from the clear difference between Figs.EZfa) and l27r b). 
It is not surprising, if we are based on the principle of adiabatic continuation, since the 
eigenstates of the LS and j-j coupling schemes are continuously connected. 

Remark that we can observe the common structure around at the value of x, in 
which singlet and triplet ground states are interchanged. Namely, essential point of 
the singlet-triplet crossing can be captured both in the two schemes. However, the F3 
non-Kramers doublet cannot be the ground state in the j-j coupling scheme, since the 
doublet in the J=4 nontet is composed of degenerate two singlets formed by F7 and Fg 
electrons. As easily understood, such singlets are energetically penalized by the Hund's 
rule interaction and the energy for F4 triplet composed of F7 and Fg electrons is always 
lower than that of the singlets. Thus, in the j-j coupling scheme, F3 non-Kramers 
doublet does not appear as the ground state except for x=0. 

Of course, if j=7/2 octet is explicitly included and A is kept finite, it is possible to 
reproduce F3 doublet. Namely, taking account of the effect of j=7/2 octet is equivalent 
to consider the local /-electron term H{, as we have done in this subsection. If we simply 
expand the Hilbert space so as to include both j=5/2 sextet and j=7/2 octet, we lose 
the advantage of the j-j coupling scheme considering only j=5/2 sextet. However, for 
an actual purpose, it is enough to consider perturbatively such effect in the order of 
1/Aso- In fact, quite recently, Hotta and Harima have shown that the result of the 
LS coupling scheme can be reproduced quantitatively even in the j-j coupling scheme, 
when the effect of j=7/2 octet is included as effective one- and two-body potentials up 

to the order of 1/Aso |21S]. 

One may claim that it is possible to reproduce the result of the LS coupling scheme 
even within the j-j coupling scheme, just by assuming that the CEF potential for J=4 
in the LS coupling scheme also works on the J=4 /^-states composed of a couple of 
/ electrons among j=5/2 sextet. However, such a procedure is not allowed due to the 
following two reasons. First it should be noted that the CEF potential is not determined 
only by the value of J. For instance, the results of the energy levels for n=7 and 13 are 
apparently different, even though both of the ground-state multiplets are characterized 
by J=7/2, since the CEF potential depends also on the values of L and S. Note that 
for n=7, S=7/2 and L=0, while for n=13, L=3 and S'=l/2. For the case of n=2, even 
if the /^-state is characterized by J=4 in the j-j coupling scheme, we cannot simply 
validate the application of the CEF potential in the LS coupling scheme to the J=4 
/^-state in the j-j coupling scheme. 

Second we should note again that the CEF effect appears only as a one-electron 
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potential. The CEF potential working on the two-electron state should be given by 
the superimposition of the one-electron potential. Thus, when we use the basis which 
diagonalizes the spin-orbit interaction, it is necessary to consider that the CEF potential 
should work on the state labeled by the z-component of j. This is the only way to define 
the CEF potential in the j-j coupling scheme, even though the Fa non-Kramers doublet 
is not reproduced. As mentioned in the above paragraph, in order to reproduce the 
results of the LS coupling scheme including the non-Kramers doublet, it is necessary 
to consider appropriately the effect of i=T /'2, octet, leading to the effective potential 
among j=5/2 states. 

5.3. Local f- electron term in the j-j coupling scheme 

In the previous subsection, we have shown the relation between LS and j-j coupling 
schemes on the basis of the local term including correctly the Coulomb interaction, 
spin-orbit coupling, and CEF potential terms. In order to make further steps to the 
construction of a microscopic Hamiltonian, let us first discuss the local /-electron state 
on the basis of the j-j coupling scheme. For the purpose, it is necessary to define the 
one /-electron state, labelled by /x, but in the j-j coupling scheme, the meaning of [i is 
clear. In the case of n<7, jj, should be the label to specify the state in the i=5/2 sextet, 
namely, the z-component of the total angular momentum j=5/2 and takes the values 
of /i=— 5/2, —3/2, ■ ■ -, 5/2. Note that for 3<n<7, j=7/2 octet is not occupied, since 
we presume that the effect of spin-orbit interaction is larger than that of the Hund's 
rule coupling in the j-j coupling scheme. On the other hand, for the case of n>7, /x 
should be considered to specify the state in the j=7/2 octet, since j=5/2 sextet is fully 
occupied. Note again that spin-orbit interaction is larger than that of the Hund's rule 
coupling. In this paper, we concentrate only on the case of n<7. Thus, in the following, 
IX indicates the z-component of the total angular momentum which specifies the state 
in the j=5/2 sextet. 

In the j-j coupling scheme, the local /-electron term should be composed of two 
parts as 

H\oc = Hqef + Hci-ch (114) 

where -ffei-ei and Hqef are Coulomb interactions among / electrons and the CEF term, 
respectively. Note that the spin-orbit interaction has been already included, when we 
define the one /-electron state in the j-j coupling scheme. In order to express -ffei-ei and 
HcEF, it is useful to define the annihilation operator of / electron in the j-j coupling 
scheme, ai^, which is related to fima with real-spin a and orbital m (=— 3,- ■ -,3) as 

a 

where Cfj_cr is the Clebsch-Gordan coefficient, given by 



c„„ = -aJ'-^^^. (m 
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with 0"=+! (—1) for up (down) real spin. 

As mentioned in the the previous subsection, in the j-j couphng scheme, we should 
take into account the CEF effect as the one-electron potential. Multi-electron state is 
obtained by simply accommodating electrons due to the balance between the Coulomb 
interaction and the one-electron potential, as has been done in d-electron systems. Then, 
the CEF term is given by 

Hci^Y = ^B^^a\^aiy, (fl7) 

where B^y is expressed by the CEF parameters for J=5/2. For the case of cubic 
structure, we can easily obtain 

5±5/2,±5/2 = 605°, 

-B±3/2,±3/2 = -18052, 

-B±l/2,±l/2 = 1205°, 

-B±5/2,t3/2 = -B;p3/2,±5/2 = 60v5i?4, 

and zero for other /x and v. For the case of tetragonal structure, we obtain 

5±5/2,±5/2 = 1050 + 605°, 



:ii8) 



-B±3/2,±3/2 = -25° - 1805°, 
5±i/2,±i/2 = -85° + 1205°, 



?° + 1205°, ^^^^^ 

-B±5/2,=f3/2 = -St3/2,±5/2 = 12a/55|, 

and zero for other /i and v. 

Note that the coefficients 5™ are, in actuality, determined by the fitting of 
experimental results for physical quantities such as magnetic susceptibility and specific 
heat. Note also that the above formulae have been obtained from the case of J=5/2. 
In general, the CEF term is expressed in matrix form, depending on the value of J; for 
J larger than 5/2, higher terms in 5™ should occur. However, as already mentioned 
above, since in this paper the effect of the CEF is considered as a one-electron potential 
based on the j-j coupling scheme, it is enough to use the CEF term for J=5/2. 

Next we consider -ffei-ei in the j-j coupling scheme. It is easy to understand that 
the Coulomb interaction term is given in the form of 

He\-e\ = - ^ /(/i,z/;z/',/i')aJ^aL«i!^'«iM'' (120) 

where / is the Coulomb interactions. The point here is the calculation of /, which is 
the sum of two contributions, written as 

/(/i, v] u', /i') = Kf^y^f,/ - Kf^y^^i^,,, (121) 

with the Coulomb integral K. The former indicates the Coulomb term, while the latter 
denotes the exchange one. It should be noted that / vanishes unless ^+v=^'+v' due 
to the conservation of z-component of total angular momentum. The matrix element 
^MiM2,/^3M4 is explicitly given by 

-'^MiM2,/i3M4 — _^^Mio"^M20-'L^^30-'t^^40--'^^_o-/2,^2-o-72,^3-o-'/2,^4-cr/2' (122) 

cr.cr' 
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where /■^ is the Coulomb matrix element among / electrons, already defined in the 
previous subsection. 

When two electrons are accommodated in the j=5/2 sextet, the allowed values 
for total angular momentum J are 0, 2, and 4 due to the Pauli principle. Thus, the 
Coulomb interaction term should be written in a 15x15 matrix form. Note that "15" is 
the sum of the basis numbers for singlet (J=0), quintet {J =2), and nontet (J=4). As 
is easily understood, this 15x15 matrix can be decomposed into a block-diagonalized 
form labeled by Jz, including one 3x3 matrix for Jz=0, four 2x2 matrices for Jz=±2 
and ±1, and four 1x1 for J2=±4 and ±3. We skip the details of tedious calculations 
for the matrix elements and here only summarize the results in the following by using 
the Racah parameters E^ {k=0,l,2) in the j-j coupling scheme |244j . which are related 
to the Slater-Condon parameters F^ as 
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For the sectors of J,=4 and 3, we obtain 



/(5/2,3/2;3/2,5/2) = Eo-5E2 



and 



1(5/2, 1/2; 1/2, 5/2) = Eo- 5^2, 
respectively. For Jz=2 and 1, we obtain 



7(3/2,1/ 
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and 

lif^, V- v\ fx') = li-u, -/x; -/i', -z/'). (132) 

By using these two relations and Eqs. ()12filll3())l . we can obtain all the Coulomb matrix 
elements |245j . 

It is instructive to understand how the /^ configuration is determined by the 
Coulomb interaction in the j-j coupling scheme. We will discuss later the local /- 
electron state determined by Hqef+Hci-cI- In the j-j coupling scheme, two electrons 
are accommodated in the j=5/2 sextet. When we diagonalize the 15x15 matrix for 
Coulomb interaction terms, we can easily obtain the eigen energies as Eq—5E2 for the 
J=4 nontet, Eq+9E2 for the J=2 quintet, and Eo-\-3Ei for the J=0 singlet. Since 
the Racah parameters are all positive, the ground state is specified by J=4 in the j-j 
coupling scheme. In the LS coupling scheme, on the other hand, we obtain the ground- 
state level as ^H with S=l and L=5 from the Hund's rules. On further inclusion of the 
spin-orbit interaction, the ground state becomes characterized by J=4, expressed as ^H4 
in the traditional notation. Note that we are now considering a two-electron problem. 
Thus, if we correctly include the effects of Coulomb interactions, it is concluded that 
the same quantum number as that in the LS coupling scheme is obtained in the j-j 
coupling scheme for the ground-state multiplet. 

In order to understand further the physical meaning of Racah parameters, it is 
useful to consider a simplified Coulomb interaction term. In the above discussion, the 
expressions using Racah parameters are not convenient, since they depend on the orbitals 
in a very complicated manner, even though they keep the correct symmetry required 
by group theory. To clarify their meanings, let us step back to the following simplified 
interaction form among i = 3 orbitals: 

-nint = U y ^ P\m\P\mi + U / ^ PimuPim'a' 

im i,(T,cr',m>m' 

+ J 2_^ JimaJim'a'Jima'Jim'a, (133) 

i,cr,cr' ,m>m' 

where Pimo-=/imcr/imo-- ^^ this equation, we include only three interactions; intra- 
orbital Coulomb interaction U, inter-orbital Coulomb interaction U', and the exchange 
interaction J. We ignore the pair-hopping J' for simplicity. Since we set J'=0 in the 
relation of U=U'+J+J', the relation U=U'+J holds among Coulomb interactions to 
ensure rotational invariance in the orbital space. 

By using Clebsch-Gordan coefficients, fima with real-spin a can be related to Oi^ as 



hma = -(y\l = O'im+u/2- (134) 

Note here that we consider only the j = 5/2 sextet. The Coulomb interaction term for 
j=5/2 is given by 

Hint = UeS Yl ^i/^^i/^' " '^'^f^'^i + (35 Jeff/4) A^;, (135) 

i^i>fj,' 
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where ni^=al Oi^, Ni=^ Ui^, Ues=U'—J/2, Jeflf=<//49, and Ji is the operator for total 
angular momentum with j=5/2. Explicitly, J? is written as 

+ (0;t0;'aS;,+iai/.av-i«iM' + </>/! '/'Jaf/.-i«iMaL'+i«iM')/2], (136) 

with 0J = ^j{j + l)-fi{fi±l) =v/35/4-/i(/i±l). 

For two electrons in the j=5/2 sextet, based upon the simplified Coulomb 
interaction term, we can easily obtain the energy levels as Ucs—^Jcsf^ for the J=4 
nontet, Ues+^SJesf^ for the J=2 quintet, and ?7efr+35 Jefr/2 for the J=0 singlet. When 
we compare these energy levels with the results obtained using Racah parameters, we 
understand the correspondence such as Egr^Ues and i?2~<^ff- Namely, Eq is the effective 
inter-orbital Coulomb interaction, while E2 denotes the effective Hund's rule coupling. 
Note that Ei does not appear, since it is related to the pair-hopping interaction which 
is not included here. 

We also note the smallness of Jqs, given as Jcff=J/49. The origin of the large 
reduction factor 1/49 is, in one word, due to the neglect of j=7/2 octet. In the Coulomb 
interaction term Eq. pHHj) . the Hund's rule term is simply written as —JS^. Note the 
relation S={gj — 1) J with gj the Lande's ^f-factor. For j=5/2, we easily obtain gj=6/7, 
indicating S=—{1/7)J. Thus, the original Hund's rule term is simply rewritten as 
-(J/49)J2. 

5.4- Level scheme in the j-j coupling scheme 

Before proceeding to the exhibition of the model Hamiltonian obtained by further 
considering the kinetic term of / electrons, it is instructive to show how the j-j coupling 
scheme works to reproduce the local level scheme of actual /-electron materials. It is an 
important point that we can resort to the analogy with the d-electron-like configuration, 
as discussed in Sec. 2. As a typical example, here we consider the /-electron state for 
the case of cubic CEF potential. 

After some algebraic calculations, we obtain two degenerate levels under the cubic 
CEF. One is T-j doublet with Kramers degeneracy and another is Fg quartet including 
two Kramers doublets. It is quite useful to define new operators with "orbital" degrees 
of freedom to distinguish two Kramers doublets included in Fg as 



/iaT = ^/W^ai-5/2 + VT76ai3/2, .^2^. 

/iai = \/5/6ai5/2 + A/l76ai_3/2, 
for "a" -orbital electrons and 

/ibt = '^i-1/2, fihi = ail/2, (138) 

for "b" -orbital electrons, respectively. The F7 state, defined as "c" orbital, is 
characterized by 

/icT = A/T76ai_5/2 - A/576ai3/2, 



/ici = A/T76ai5/2 - A/576ai„3/2. 



(139) 
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Figure 29. Views for (a) r|, (b) r|, and (c) Fy orbitals. T^ indicates the irreducible 
representation of point group in Bethe's notation |242j . 



For the standard time reversal operator }C=—iayK, where K denotes an operator to 
take the complex conjugate, we can easily show the relation 

ICfira = (Tf\r-a- (140) 

Note that this has the same definition for real spin. 

In Fig. |2ni we show the shape of three orbitals. As intuitively understood from the 
shape of Fy orbital, this keeps the cubic symmetry, indicating A- representation. In fact, 
in the group theory, it is characterized by a^. Note that the subscript "u" indicates 
ungerade, since we consider / electron with £=3. On the other hand, degenerate Fg 
and Fg orbitals seems to be similar to x^ — y"^ and ?>z^ — r^ orbitals of 3(i electrons, 
respectively, indicating E- representation. In the group theoretical argument, these are 
classified into Cu- Concerning the similarity between Fg and Cg orbitals, it is quite natural 
from a mathematical viewpoint, since we recall the fact that Fg is isomorphic to F3 x Fg, 
where F3 indicates E representation for the orbital part and Tq denotes the spin part. 
This point is quite impressive when we consider the orbital physics for d- and /-electron 
systems. Namely, by exploiting this mathematical similarity, it is possible to understand 
the complex /-electron phenomena with the use of the microscopic Hamiltonian in 
common with that of the d-electron multiorbital model. We will see later this point in 
the construction of the model Hamiltonian. 

Now we discuss the /-electron configuration in the F7 and Fg levels in the manner in 
which we have considered the d-electron configuration. First, we pick up the AuCus-type 
cubic crystal structure. A typical AuCua-type material with one / electron per site is 
is Celna, in which F7 and Fg are the ground and first excited states, respectively |24(ij . 
If we accommodate one more electron to consider the p configuration, immediately 
there appear two possibilities, "low" and "high" spin states, as we have discussed in the 
d-electron configuration. When the CEF splitting energy between F7 and Fg levels is 
smaller than the Hund's rule coupling, the second electron should be accommodated in 
the Fg levels. In the situation in which one is in the F7 and the other in the Fg, a F4 
triplet appears for the p state in the j-j coupling scheme. As mentioned in the previous 
subsection, F3 non-Kramers doublet does not appear in the j-j coupling scheme. On 
the other hand, if the CEF splitting is larger than the Hund's rule interaction, then the 
P ground state is formed from two F7 electrons, leading to a Fi singlet state. When we 
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Figure 30. Electron configurations in tlie j-j coupling scheme for (a) Celna, (b) 
Prlus, and (c) Ndlna. T^ denotes the irreducible representation of point group in 
Bathe's notation 



compare this Fi state with that in the LS couphng scheme, we notice that it is given 
by a mixture of J=0 and J=4 states, but the J =4 component is found to be dominant. 
Note also that Fi is the antisymmetric representation of F7 x F7. 

Since we do not know the exact value of the Hund's rule interaction in /-electron 
compounds, it is difficult to determine the /^ state by purely theoretical arguments. 
In this case, we have to refer to the data on actual materials. Fortunately, we have 
the example of Prlus, a typical /^ material with AuCus-type crystal structure. From 
several experimental results, Fi has been confirmed to be the ground level in Prina 
|247| I248J . Thus, the low-spin state should be taken for the AuCus-type structure in 
the j-j coupling scheme. 

Here the reader may pose a naive question: Is the Hund's rule interaction really so 
small in /-electron systems? However, we have already discussed this point in the 
previous subsection. Namely, the effective Hund's rule interaction Jefr is given by 
Jeff=J/49 in the j-j coupling scheme, where J is the original Hund's rule interaction 
among / electrons. Note again that the magnitude of the Hund's rule interaction is 
effectively reduced by the factor 1/49 in the j-j coupling scheme. Even if J=leV, Jgfr is 
reduced to be about 200K, which is comparable with the CEF splitting energy. Thus, 
it is possible to have the low-spin state in the j-j coupling scheme. 

Next, we take a further step to the /^ state by adding one more / electron. Since 
F7 is fully occupied to form Fi, the next electron should be placed in the Fg state, 
as shown in Fig. IHUT c). clearly indicating that there exists an active orbital degree of 
freedom. The /^ state composed of two F7 and one Fg electron is expressed as Fg in the 
terminology of group theory. When we again consider actual materials, Ndlna is found 
to be a typical /^ material with the AuCua-type crystal structure. In experiments, it 
has been confirmed that Fg is the ground level |249| I25()[l25lj . as we have found with 
the present j-j coupling scheme. 

Let us turn our attention to another crystal structure, in which Fg is lower than 
F7 in the /^ configuration. Typical materials are the rare-earth hexaborides RBg with 
R=Ce, Pr, and Nd. As is well known, the ground level of CeBg is Fg, indicating 
that the quadrupolar degree of freedom plays an active role in this material |252j . In 
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Figure 31. Electron configurations in the j-j coupling scheme for rare-earth 
hexaborides, (a) CeBg, (b) PrBg, and (c) NdBg. F^ is the irreducible representation 
of point group in Bcthc's notation |242j . 



fact, anomalous behavior related to quadrupolar ordering has been suggested by several 
experimental results. 

First, we note that the level splitting between Fg and F7 is assumed to be larger 
than the Hund's rule interaction. When we accommodate two electrons in Fg orbitals, 
the triplet (F5), doublet (F3), and singlet (Fi) states are allowed. Among these, owing 
to the effect of the Hund's rule interaction, even if it is small, the F5 triplet should be 
the ground state. This has actually been observed in PrBg |253|I254] . Further, in order 
to consider NdBg, another electron is put into the Fg orbital, making a total of three. 
Alternatively, we may say that there is one hole in the Fg orbital. Such a state is found, 
again, to be characterized by Fg . Experimental results on NdBg have actually been 

('9^ a^-^p-n p^i^p^p- p^i^--n i^-^p^sa 

reported which lead to the ground state of Fg |253| I255[ I256| 125 7j . Thus, when Fg is 
the ground state for the one /-electron case, we obtain F5 for the /^ and Fg for the /^ 
configurations. 

We have shown that the ground states deduced from the j-j coupling scheme are 
consistent with experimental results. However, in order to explain the experimental 
results quantitatively, it is unavoidable to analyze the CEF levels using the LS coupling 
scheme. As mentioned above, it is possible to reproduce the result of the LS coupling 
scheme by considering the effective potentials from the j=7/2 octet, but what we would 
like to stress here is that even in a localized system, the symmetry of the ground level can 
be understood via the simple j-j coupling scheme. We need to recognize the limitations 
of the simple j-j coupling scheme when we treat a local electronic state. For instance, to 
consider the /^ state, we simply put three electrons into the CEF level scheme which is 
determined with the /^ configuration. Thus, the wavefunction of the /^ state is uniquely 
determined. However, in an actual situation, the dectet labeled by J=9/2 (L=6 and 
S=3/2) is split into two Fg and one Fg orbital. The ground-state wavefunctions will then 
depend on the two CEF parameters B^ and Bq |258j . In order to explain experimental 
results on localized /-electron materials, one should analyze the Hamiltonian which also 
includes the complex effective potentials from j=7/2 octet. In this paper, however, 
the electronic states are considered with an itinerant picture based on the simple j-j 
coupling scheme. Thus, it is important to check that the local electronic state formed 
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by / electrons in this way is consistent with the symmetry of the state obtained with 
the LS couphng scheme. 

In summary, it has been shown that the ground states of the /^ and /^ 
configurations can be quahtatively reproduced by accommodating / electrons in the 
CEF levels of a corresponding /^ material, provided that the CEF level splitting is 
larger than the Hund's rule interaction. Thus, the j-j coupling scheme works even in 
the localized case. Accordingly, we believe that a microscopic theory can be developed 
in which we discuss the magnetism and superconductivity of /-electron compounds in 
terms of the j-j coupling scheme. 

5.5. Model Hamiltonian 

In previous subsections, we have explained in detail that the j-j coupling scheme works 
for the local /-electron state. Now let us include the effect of kinetic motion of / 
electrons [Z21- In this article, we are based on a itinerant picture for / electrons. From 
our experience, this picture seems to be valid, when we consider actinide compounds, 
in particular, heavy actinides. On the other hand, for rare-earth materials, it has been 
firmly believed to take the localized picture for / electrons. In particular, in order to 
consider the heavy fermion behavior, it is indispensable to consider the system including 
both the conduction electron with wide band and the almost localized / electron, which 
are hybridized with each other. 

It is believed that the hybridization of / electrons with conduction electron band 
is important to understand the magnetism of /-electron systems. In fact, in the 
traditional prescription, first we derive the Coqblin-Schrieffer model from the periodic 
Anderson model by evaluating the c-f exchange interaction J^i within the second-order 
perturbation in terms of the hybridization between /- and conduction electrons. Then, 
we derive the RKKY interactions again using the second-order perturbation theory 
with respect to J^. In general, the RKKY interactions are orbital dependent and 
interpreted as multipole interactions. Such orbital dependence originates from that 
of the hybridization. Note that the hybridization should occur only between /- and 
conduction band with the same symmetry. Here we emphasize that the symmetry of 
/-electron state is correctly included in our calculations. Thus, the structure in the 
multipole interactions will not be changed so much, even if we consider the effect of 
hybridization with conduction band, as long as we consider correctly the symmetry of 
/ electron states. 

In this paper, we consider the tight-binding approximation for /-electron hopping 
motion. The kinetic term can be written as 

^kin = Y^ t;^,at^ai+a., (141) 

i,a,/i,i^ 

where tj^^, is the overlap integral between the /i- and //-states connected by the vector a. 

Let us consider the hopping motion of / electrons based on the j-j coupling scheme. 

In order to evaluate tj^^, which is hopping of / electrons between the /i-state at i site 
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and the i/-state at i + a site, again it is convenient to step back to /-electron operators 
in the i=3 multiplet, defined as fima- Since the real spin should be conserved in the 
hopping process, tj^^, is given as 

tfiiy = ^^Ct^crC,yaT^^f_i-a/2-3,iy-a/2^ (142) 

a 

where T^^.^, ^, is the hopping amplitude of electrons between {i, m)- and (£', m')-states 
along the a-direction. 

Now the problem is reduced to the evaluation of T^^.^,^,. Although we can 
simply consult the paper of Slater and Koster |1H5[ I259J , a convenient formula has been 
obtained by Sharma for the overlap integral between two orbitals, (£, m) and {£', m') 
[26Qj. connected by unit vector a. It is expressed as 



Att / Anr 

Tt,m;l',m' = (^^V) ^ ^^^ ^ ^^^^F/^l^, ^)F,w(^, ^), (143) 

where {ii'a) denotes Slater's two-center integral through the a bond, for instance, it is 
{ffcr) for i=i'=3 and {fpc) for £=3 and i'=l. 9 and ip are polar and azimuth angles, 
respectively, to specify the vector a as 

a= (sin 6' cos (y9, sin 6' sin y), cos 6'). (144) 

Here we consider the hopping between / orbitals in nearest neighbor sites by putting 
£=£'=3. After some algebraic calculations, we obtain the hopping amplitudes as follows. 
For diagonal elements, we obtain 

^±5/2,±5/2 = StoSm e, 

^±3/2,±3/2 = tosm^e{l + lhcos'^e), (145) 

^±i/2,±i/2 = 2to(l-2cos2^ + 5cos^^), 
where the energy unit to is given by 

to = (3/56)(//a). (146) 

Here {f fu) is the Slater-Koster two-center integral between adjacent / orbitals. Note 
that t^_ =0. For off-diagonal elements, we obtain 

^±5/2,±i/2 = -toyiOeT^''^sin2^(l-3cos2 0), 

i^5/2!T3/2 = tov^e^^^^sin^^, (147) 

^±i/2,T3/2 = -to\/2e'^^'^siYi^e{l + hcos^e), 



and 



^5/2,-1/2 = -^1/2,-5/2 = toVlQe-^''^ sin^ e sin 26, 
t^/9-,/9 = -t't./o./o = -2toV5e-'^sm'^esm2e, (148) 



''5/2,3/2 = -^-3/2,-5/2 = -^^OVoe ^smf7SmzC7, 

"1/2,3/2 ''-3/2,-1/2 



^1/2 3/2 = -i'^3/2-1/2 = ^ov^e^'^sin2^(l-5cos2| 



Note that tl^=tll. 
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5.6. Fg model 

In the previous subsection, we have completed the construction of a model Hamiltonian, 
which is expected to be a basic model to investigate the microscopic aspects of 
magnetism and superconductivity of /-electron systems. Since it includes six states 
per site, i.e., three Kramers doublets, the analysis may be difficult. Of course, even 
if the calculations seem to be tedious, it is necessary to carry out analytical and/or 
numerical research on the basis of such a three orbital model. However, for practical 
purposes, it is convenient to simplify the model, if possible. In this subsection, as an 
effective model for actinide compounds, we introduce a Fg model, by discarding F7 
orbital [12]. 

A simple explanation to validate the ignorance of F7 is to assume large CEF splitting 
energy between F7 and Fg levels. This simplification is motivated by the fact that 
the possibility of exotic octupole ordering has been actively discussed in Ce^Lai-^jBg 
and Np02 with Fg ground state. Here readers may be doubtful of the reality of our 
assumption, since the Coulomb interaction among / electrons is naively thought to be 
larger than the CEF level splitting in any case. However, it should be noted again that 
we are now considering the /-electron state in the j-j coupling scheme, not in the original 
/-electron state with angular momentum i=3. As already mentioned, the Hund's rule 
interaction in the j-j coupling scheme is effectively reduced to be 1/49 of the original 
Hund's rule coupling. Even when the original Hund's rule coupling among / electrons 
is 1 eV, it is reduced to 200 K in the j-j coupling scheme. For instance, the CEF level 
splitting in actinide dioxides is considered to be larger than 1000 K. We also recall that 
the CEF level splitting in CeBg is as large as 500 K. Thus, we safely conclude that our 
present assumption is correctly related to some actual materials. Of course, in order 
to achieve quantitative agreement with experimental results, it is necessary to include 
also F7 level, since the magnitude of the CEF splitting is always finite, even if it is large 
compared with the effective Hund's rule interaction. However, we strongly believe that 
it is possible to grasp microscopic origin of unconventional superconductivity as well as 
spin and orbital, i.e., multipole, ordering in /-electron systems on the basis of the Fg 
model, since this model is considered to be connected adiabatically from the realistic 
situation. It is one of future tasks to develop more general theory to include all the 
j=5/2 sextet in future. 

Concerning the /-electron number, typically we treat the case with one / electron 
in the Fg multiplet per site. However, this restriction does not simply indicate that 
we consider only the Ce-based compound. In the j-j coupling scheme, in order to 
consider /"-electron systems, where n indicates local / electron number per site, we 
accommodate / electrons in the one-electron CEF levels due to the balance between 
Coulomb interactions and CEF level splitting energy, just as in the case of d-electron 
systems. Thus, as shown in Fig. IH^ a). the Fg model is applicable to the cases for 
n = 1 ~ 4 in the Fg-F7 system, where Tr^-Ty symbolically denotes the situation with F^. 
ground and Ty excited states. Furthermore, we should note that due to the electron-hole 



Orbital ordering phenomena in d- and f -electron systems 



88 



(a) 



r^ 



r.^^ 



W) 



U: 



^7 



rs^\i 



r; n 



^7 






i*f 



i 



% 



^7 



f 






t 



Figure 32. Electron configurations in the j-j coupling scheme to which the Fg model 
is applicable, (a) /^ ~ f^ configuration when Fg is lower, (b) /^ ~ /^ configuration 
when Fy is lower, 



symmetry in the Fg subspace, the Fg model is also applicable to the cases for ra = 3 ~ 6 
in the Fy-Fg system, as shown in Fig. |22fb). 

Now let us define the Fg model. First we consider the hopping part. For simplicity, 
here we set the cubic-based lattice. Later we will discuss the hopping amplitude of 
other lattice structures. We include the hopping in the xy plane and along the z- 
axis for a=x= [1,0,0], y= [0,1,0], and z= [0,0,1], respectively. To evaluate the hopping 
amplitude, we simply set {9, ip) to be (7r/2,0), (7r/2,7r/2), and (0,0) for x, y, and z 
directions. Then, by using the general results in the previous section, we easily obtain 
t^^ between neighboring / orbitals in the xy plane and along the z axis. Further we 
transform the basis by the above definitions for Fg operators with orbital degrees of 
freedom. The results are given as 



t-., = t 



3/4 -v^/4 \ 
-v^/4 1/4 )' 



for the x-direction, 



t^ 



3/4 v^/4\ 
" v^/4 1/4 ]' 



for the y direction, and 



t 




(149) 



(150) 



(151) 



for the z direction. Note that t=8tQ={3/7){ffa). Here we stress that the hopping 
amplitudes among Fg orbitals are just the same as those for the Cg orbitals oiSd electrons. 
Intuitively, readers can understand this point due to the shapes of Fg orbitals as shown 
in Fig. 123 which are similar to Cg orbitals. 
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After transforming the basis of the Coulomb interaction terms, the Hamiltonian for 
Fg orbitals is given by 

^rs = XI ^rr'/ia/i+arV + U^ P\r]Piri + U' ^ piaPib 
i,a,(T,T,r'' i,T i 

+ JY1 fUL'i\a.'h. + J'Y. fllfllf-^r'lfir'^, (152) 

where pira= firafira and piT-=^^pio.^. In the Coulomb interaction terms, U, U', J, and 
J' denote intra-orbital, inter-orbital, exchange, and pair-hopping interactions among Fg 
electrons, respectively, expressed by using the Racah parameters Ek as 

U = Eo + E^ + 2E2, 

U' = Eo + {2/3)E,, 

J = 5E,, ^^^^' 

J' = Ei-{11/3)E2. 

Note that the relation U=U'+J+J' holds, ensuring rotational invariance in pseudo- 
orbital space for the interaction part. For rf-electron systems, one also has another 
relation J=J\ as mentioned in Sec. 2. When the electronic wavefunction is real, this 
relation is easily demonstrated from the definition of the Coulomb integral. However, 
in the j-j coupling scheme, the wavefunction is complex, and J is not equal to J', in 
general. For simplicity, we shall assume here that J=J', noting that essential results 
are not affected. Since double occupancy of the same orbital is suppressed owing to the 
large value of U, pair-hopping processes are irrelevant in the present case. 

We believe that this Fg Hamiltonian provides a simple, but non-trivial model to 
consider superconductivity and magnetism in /-electron systems. Note again that it is 
essentially the same as the model for Cg electron systems such as manganites, although 
the coupling with Jahn- Teller distortion is not included in the present model. Due to the 
complex interplay and competition among charge, spin, and orbital degrees of freedom, 
a rich phase diagram has been obtained for manganites. Thus, it is definitely expected 
that a similar richness will also be unveiled for /-electron systems based on the Fg model 
Hamiltonian. 

6. Orbital physics in /-electron systems 

We have constructed a microscopic model Hamiltonian for /-electron systems in the 
previous subsection. In particular, we could obtain the Fg orbital degenerate model 
as an effective Hamiltonian for actinide compounds. In this section, we review the 
theoretical results of the spin and orbital structure, i.e., multipole order, based on the 
Fg model. 

6.1. Spin and orbital structure of actinide compounds 

In this subsection, we review the theoretical effort to understand the magnetic structure 
of uranium and neptunium compounds with HoCoGas-type tetragonal crystal structure 
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Figure 33. (a) Crystal structure of AnTGas. Schematic views of magnetic structures 
at low temperatures, composed of magnetic moments of U and Np ions for (b) UNiGas, 
(c) UPtGas, (d) NpFeGas, (e) NpCoGas, and (f) NpNiGag. For NpFeGag, magnetic 
moments at Fe sites are also depicted. 



The crystal structure is shown in Fig. lH^ a). The varieties of magnetic structure 
of U-115 and Np-115 explained in Sec. 1 are summarized in Figs. l33f b)- l33f f ) In order to 
set up the microscopic model for actinide 115 materials, it is useful to consider UGas and 
NpGa3, which are the mother compounds of U-115 and Np-115. Among them, it has 
been reported that UGas exhibits a G-type AF metallic phase in the low-temperature 
region |261j . but a "hidden" ordering different from the magnetic one has been suggested 
by resonant X-ray scattering measurements |262j . Unfortunately, orbital ordering in 
UGas is not yet confirmed experimentally, but it may be an interesting possibility to 
understand the result of resonant X-ray scattering experiment on UGas based on the 
orbital-ordering scenario. 

Although it is difficult to determine the valence of actinide ions in the solid state, 
for the time being, we assume that the valence is U^"*" or Np^+, including three or four 
/ electrons per ion. By considering the CEF potential and Coulomb interactions, we 
then assign three or four electrons to the states in the j=5/2 sextet. In order to proceed 
with the discussion, it is necessary to know which is lower, Ti or Fg, in the one /- 
electron picture. For some crystal structures it is possible to determine the level scheme 
from intuitive discussions of /-electron wavefunctions and the positions of ligand ions. 
However, this is not the case for the AuCua-type crystal structure. For this case, we 
again invoke experimental results on Celns, a typical AuCus-type Ce-based compound, 
where Fy and Fg have been reported as ground and excited states, respectively, with 
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Figure 34. Level schemes for (a) Celna, (b) UGas, and (c) NpGaa based on the 
j-j couphng scheme. Here we assume trivalent actinide ions as U'^"'" (5/^) and Np'^+ 
(5/'*). It should be noted that up and down arrows denote pseudospins to distinguish 
the states in the Kramers doublet. Note also that for NpGaa, a couple of electrons in 
Fg orbitals form a local triplet, leading to T^. 



an energy difference of 12nieV |246j . Tlius, we take Fy to be lower for tlie present 
considerations, as sliown in Fig. lH^ a). 

In tlie j-j coupling scheme for UGaa and NpGaa, we accommodate three or four 
electrons in the one-electron energy states Fy and Fg. We immediately notice that 
there are two possibilities, i.e., low- and high-spin states, depending on the Hund's rule 
interaction and the splitting between the Fy and Fg levels. As discussed in the previous 
subsection, the effective Hund's rule interaction can be small in the j-j coupling scheme 
and thus, the low-spin state should be realized, as shown in Figs. IH^ b) and (c). We 
emphasize that this low-spin state is consistent with the LS coupling scheme. In fact, 
for NpGas, the observed magnetic moment at Np ion has been found to be consistent 
with F5 triplet |263j . 

In the electron configuration shown in Figs. IMT b) and (c), the Fy level is fully 
occupied to form a singlet. If this Fy level is located well below the Fg, the occupying 
electrons will not contribute to the magnetic properties. Thus, we can ignore the Fy 
electrons for our present purposes. In order to validate this simplification, it is useful to 
introduce the results of band-structure calculations for Celus |264j and UGas j265j . Note 
that both results have been obtained assuming the system is in the paramagnetic state. 
In order to focus on the / electron components of the energy band, we concentrate 
on the bands around the F point near the Fermi level. For Celus, the energy band 
dominated by Fy character is found to be lower than the Fg-dominated band, consistent 
with the local level scheme in Fig. IMT a). An important point is that the Fermi level 
intersects the Fy-dominant band, indicating that the Fermi surface is mainly composed 
of Fy electrons hybridized with Ga-ion p electrons. On the other hand, for UGas, the 
Fy band is also lower than the Fg band, but the Fermi level crosses the Fg band. Thus, 
the Fy band appears to be fully occupied, consistent with the j-j coupling level scheme, 
as shown in Fig. IMT b). Since the main contribution to the Fermi surface comes from 
Fg electrons, it is natural to dwell on the Fg bands and ignore the occupied Fy bands in 
giving further consideration to many-body effects. 
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So far, we have considered the model in the cubic system, but as mentioned before, 
115 materials exhibit tetragonal crystal structure. To include the effect of tetragonality, 
here we introduce two ingredients into the model Hamiltonian. One is non-zero A, 
which is the level splitting between two orbitals. Under the tetragonal CEF, the local 
electronic levels are given by two T^ and one Tq states. Among them, Tq is just equal to 
Fg in the cubic system. Two T^ states are given by the linear combinations of ai_|_3/2|0) 

and Oizr5/2|0)' which can be expressed also by the mixture of Tj and Fg. Here for 
simplicity, we introduce A, splitting energy between Fg orbitals, by ignoring the change 
of wavefunctions from cubic to tetragonal systems. Another is the change in the hopping 
amplitude along the 2;-axis. In AnTGas (An=U and Np), AnGas layer is sandwiched 
by TGa2 blocks, as shown in Fig. l33f a). indicating that the hopping of /-electron along 
the 2;-axis should be reduced from that in AnGas. However, it is difficult to estimate 
the reduction quantitatively, since it is necessary to include correctly the hybridization 
with d-electrons in transition metal ions and p-electrons in Ga ions. Thus, we change 
the hopping t as t^ in the definition of t^^,. 

Then, the Hamiltonian is the sum of iifrg and the level splitting term, given by 

H = Hr,-A J^iPia - Pib)/2. (154) 

i 

Concerning the hopping amplitudes in the xy plane, they are given by Eqs. p49|) and 
()150|) . but along the z-axis, it is necessary to include the change from cubic to tetragonal 
case. Namely, the effective hopping along the z axis is expressed as 



trr'=tA . , , (155) 




where t^ is the reduced hopping amplitude along the z-axis. The ratio tz/t is less 
than unity. We note that in actuality, A should be related to the value of t^, since 
both quantities depend on the lattice constant along the z axis. However, the relation 
between tz and A is out of the scope at present and thus, here we simply treat them as 
independent parameters. 

Among several methods to analyze the microscopic model, we resort to an exact 
diagonalization technique on a 2x2x2 lattice. Although there is a demerit that it is 
difficult to enlarge the system size, we take a clear advantage that it is possible to deduce 
the magnetic structure by including the effect of electron correlation. In order to discuss 
the ground-state properties, it is useful to measure the spin and orbital correlations, 
which are, respectively, defined by 






S{q) = {l/N)}_^{a?a;)e''^-^^-^\ (156) 

with cr?=X;^(^irT-^iri)/2, and 

nq) = il/N)J2{rlr,ne"^-^'-^\ (157) 

with r[=J2ai^iaa-—niba)/'2- Here N is the number of sites. 
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Figure 35. (a) Phase diagram for UGaa obtained by the exact diagonaUzation. The 
region of J>U' is ignored, since it is unphysical. See Fig. |S{a) for the definitions of 
abbreviations. Here "PM-G" indicates the PM phase with enhanced (tt, tt, tt) spin 
correlation, (b) Phase diagram of the magnetic structure in the (A,tz) plane for J—0 
and U' = 3.5. (c) Ferro orbital pattern in the A-type AF phase, (d) Antifcrro orbital 
pattern in the G-type AF phase. 



6.1.1. U-115 Let us review the results for n=l [H2], in which one electron is 
accommodated in Fg orbital, corresponding to uranium compounds. First we consider 
the cubic case (A=0 and tz/t=l). After we evaluate spin and orbital correlations for 
several parameter sets, the ground-state phase diagram is completed on the (f/', J) 
plane, as shown in Fig. IHKT a). In the small- J region, the paramagnetic (PM) phase 
exists for large parameter space and in the boundary region between PM and G-type 
AF states, we can see the PM phase with dominant (vr, vr, vr) spin correlation. Note that 
such a PM-G region is not specific to the case of J=0, since it appears even when we 
increase the Hund's rule interaction. 

Here we briefly discuss the phases in the large- J region. We observe an interesting 
similarity with the phase diagram for undoped manganites RMnOa, in which mobile 
Cg-electrons are tightly coupled with the Jahn- Teller distortions and the background t2g 
spins. Note that the present Hamiltonian is just equal to the Cg electron part of the 
model for manganites [72]. In the so-called double-exchange system with large Hund's 
rule coupling between Cg and t2g electrons, the Jahn- Teller distortion suppresses the 
probability of double occupancy and it plays a similar role as the interorbital Coulomb 
interaction U'. The AF coupling among t2g spins, Jaf, controls the FM tendency in 
the Cg-electron phases. Roughly speaking, large (small) Jaf denotes small (large) J. 
Then, we see an interesting similarity between Fig. IHsT a) and the phase diagram for 
manganites, except for the PM region. See Figs. IHl^b) andEfa). In particular, a chain of 
the transition, FM— i^A-AF— >C-AF^G-AF, occurs with decreasing J (increasing Jaf)- 
Again we stress that the present Fg model for /-electron systems is essentially the same 
as the Cg orbital model in the d-electron systems. It is interesting to observe common 
phenomena concerning orbital degree of freedom in /-electron systems. 

Now we include the effect of tetragonality. Since the mother compound UGas is AF 
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metal, in our phase diagram, it is reasonable to set the parameter region corresponding 
to "PM-G". Then, we choose f/'=3.5 and J=0. Again we evaluate spin and orbital 
correlations by changing A and tz, and obtain the phase diagram in the {A,tz) plane, 
as shown in Fig. ESfb) for J=0 and U'=3.5. Note that the ground state for A=0 and 
tz=l is magnetic metallic, as seen in Fig. IHHT a). It is found that an A-type AF phase 
appears in the negative A region for tz>0.68. Note that the appearance of the A-AF 
phase is not sensitive to t^ as long as tz>0.68. Rather, A seems to play a key role 
in controlling the change of the magnetic phase. Here we recall the experimental fact 
that UNiGas exhibits a G-type AF phase, while UPtGas shows an A-type. Thus, it is 
necessary to relate the effect of A to the difference in magnetic structure found between 
UNiGas and UPtGas. Although t^ may differ among U-115 compounds, we focus here 
on the effect of A. 

From the orbital correlation, we obtain the ferro-orbital and anti-ferro orbital 
patterns for A-AF and G-AF phases, respectively, as shown in Figs. IH^T c) and lHM d). 
Let us now discuss the reasons for the appearance of an A-AF phase. For negative 
values of A, we easily obtain ferro-orbital (FO) pattern composed of Fg orbitals, as 
illustrated in Fig. ItlKT c). For electrons to gain kinetic energy of motion along the z-axis, 
it is necessary to place the AF spin arrangement along this same axis. In the FM spin 
configuration, electrons cannot move along the 2;-axis due to the Pauli principle, since 
hopping occurs only between Fg orbitals along the z-axis. On the other hand, in the xy 
plane 6-orbital electrons can hop to neighboring a-orbitals with a significant amplitude, 
which is larger than that between neighboring 6-orbitals. Thus, in order to gain kinetic 
energy, electrons tend to occupy a-orbitals even in the FO state composed of 6-orbitals, 
as long as |A| is not so large. When we explicitly include the effects of the Hund's 
rule interaction J, electron spins should have FM alignment between neighboring sites 
in order to gain energy in hopping processes from b- to a-orbitals. Consequently, a 
FM spin configuration is favored in the xy plane. In the case with antiferro orbital 
correlations, spin correlation tends in general to be FM, as has been widely recognized 
in orbitally degenerate systems. 

Here we mention a relation of A to the magnetic anisotropy in U-115 materials. 
For UPtGas with the A-AF phase, Xa is larger than Xc, whereas this anisotropy is not 
pronounced in UNiGa5 with the G-AF phase jSHl- An analysis for the high-temperature 
region based on LS coupling yields the Jz=±l/2 Kramers doublet as the ground state 
among the dectet of J=9/2 (L=6 and 5=3/2). The states with J^=±l/2 in the LS 
coupling scheme have significant overlap with /ifet/ict/icjO) and /ibj/ict/icjO) i^ the j-j 
coupling scheme. Accordingly, by the present definition A should be negative to place 
Fg below Fg. If the absolute value of A(<0) becomes large, Fg is well separated from 
Fg and the magnetic anisotropy will consequently become large. Thus, a change from 
G- to A-type AF phase is consistent with the trends of magnetic anisotropy in UNiGas 
and UPtGas. 

Finally, we make a brief comment about the effect of t^. Following the above 
discussion, the A-AF phase should appear even for small t^- However, in the present 
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calculation it disappears for t2<0.68, a critical value which seems to be rather large. 
Such a quantitative point depends on the system size, and we note that it is necessary 
to perform the calculation in the thermodynamic limit. 

While such investigations are just beginning, we already see a number of 
opportunities for future work along this path. Concerning issues directly related to 
the present context, it is highly recommended that calculations be carried out in the 
thermodynamic limit, in order to confirm the present exact diagonalization results. 
For instance, the magnetic susceptibility should be evaluated in the random phase 
approximation or fluctuation-exchange method. With such an approach, the magnetic 
structure can be discussed by detecting the divergence in the magnetic susceptibility. 
This is one of our future tasks. Another problem is how to establish the effective 
reduction of tz in considering the case of UTGas. In such systems, TGa2 blocks are 
interspersed between UGas layers, but the main process may occur through the Ga 
ions. To analyze this, it is necessary to treat a three-dimensional f-p model with explicit 
consideration of U and Ga ions. This is another problem for future investigation. 

6.1.2. Np-115 Now we review the theoretical results for magnetic structure of Np-115 
[5^ . First we consider the case of n=2, which is corresponding to the trivalent Np 
ion. At tz=1 and A=0 (cubic case), local triplet composed of a couple of / electrons 
is formed at each site and the G-type AF structure is stabilized due to the so-called 
superexchange interaction. Even when A is introduced as the tetragonal CEF effect, 
the G-AF structure remains robust for |A| < 1. When |A| is larger than unity, two 
electrons simultaneously occupy the lower orbital, leading to the non-magnetic state 
composed of local Fi, irrelevant to the present study to consider the magnetic phase. 
When we change t^ for A=0, again the G-type AF structure is stabilized, but we find 
that the spin correlation of g=(7r,7r, 0) comes to be equivalent to that of q=(7r,7r,7r) 
with the decrease of t^, since the AF structure is stabilized in each xy plane due to 
superexchange interaction and the planes are decoupled for small tz. 

At the first glance, it seems difficult to understand the variety of magnetic phases 
observed in NpTGas even in a qualitative level, when we consider only the trivalent Np 
ion. However, there is no a priori reason to fix the valence as Np'^"*". In NpTGas, d- 
electron band originating from transition metal ions may significantly affect the valence 
of Np ion ISO! • In addition, we also stress that the actual compounds exhibit AF metallic 
behavior. In the band-structure calculation, the average number of / electrons at Np ion 
is easily decreased from four. Thus, we treat the local /-electron number as a parameter. 

We may consider another reason to decrease effectively the number of / electron 
from n=2 in NpGas. In the present two-orbital model, the G-AF structure is robust, 
which is natural from the theoretical viewpoint within the model. However, in the 
experimental result on NpGas, the low-temperature ground state is ferromagnetic, 
although the AF phase has been observed around at T~60K. In order to understand 
the occurrence of the FM phase in the two-orbital model, it is necessary to inject some 
amount of "hole" in the AF phase, since the double-exchange mechanism works to 
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Figure 36. (a) Ground-state phase diagram in the ([/', J) plane for n—\.h, t^^l, and 
A=0. Ground-state phase diagrams of the magnetic structure in the (A, t^) plane for 
n^l.5, C/'=5, (b) J=0, (c) J=0.5, and (d) J=4. 



maximize the kinetic motion of electrons. It is difficult to determine the amount of 
doped holes to obtain the FM phase, but at least qualitatively, the effective decrease of 
n seems to be physically meaningful in NpGaa as well as NpTGas. 

Then, we consider the case of n=1.5. In Fig. lHUT a). we show the ground-state phase 
diagram in the {If, J) plane at n=1.5 for the cubic case with tz=l and A=0. At J=0, 
a G-type AF structure is stabilized due to superexchange interaction in the same way 
as the case of n=2. However, the G-AF structure is immediately changed to a C-AF 
structure only by a small value of the Hund's rule interaction. With increasing J, the 
magnetic phase changes in the order of G-AF, C-AF, A-AF, and FM phases, except for 
the C-AF phase in the large J region. Concerning the spin structures, see Fig. |Hfa). 
This result is quite natural, since we are now considering the magnetic structure based 
on the two-orbital model, in which the FM tendency is due to the optimization of kinetic 
motion of electrons. 

After calculations of the spin and orbital correlations for several parameter sets, 
we obtain the ground-state phase diagram in the (A,^^) plane, as shown in Figs. IHHT b). 
for f/'=5 and J=0. In the region of large positive A, we find that G-AF(I) phase with 
ferro-type arrangement of Fg orbital extends in a wide range of the phase diagram. It is 
found that the C-AF(I) phase appears in the region for small positive A and 0.5<tz<l, 
The C-AF(I) phase exhibits the dominant component of (vr, vr, 0) in the spin correlation. 
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When we further decrease A, we find G-AF(II) phase, which may be considered as 
orbital disordered, since there is no dominant component in the orbital correlation. For 
small tz and small negative A, we find another C-AF phase, which we call C-AF(II), in 
which the spin correlation of (vr, 0, vr) and (0, vr, vr) are dominant. For small t^ and large 
negative A, there appears yet another G-AF phase, called G-AF(III) with ferro-type 
arrangement of Fg orbital. In any case, for J=0, we can observe several kinds of C- and 
G-AF phases, but A-AF phase does not occur. 

Although we increase the value of J as J=0.5, no new phases appear in the phase 
diagram for n=1.5, as shown in Fig. IHUT c). There are three phases, but they are two 
C-AF and one G-AF states. As labeled explicitly in the phase diagrams, C-AF (I), C- 
AF(II), and G-AF(I) are the same as those in the phase diagram of Fig. EH^b). Due to 
the effect of J, G-AF (II) and G-AF (III) disappear, since the number of FM bond should 
be increased to gain the kinetic energy. As shown in Fig. lHHT d). when we further increase 
the value of J as J=4, the G-AF phase completely disappears and instead, we observe 
the A-AF phase sandwiched by two C-AF phases. By analogy with the various phases 
of manganites, the A-AF phase is considered to appear due to the double-exchange 
mechanism in the two-orbital model, when J is increased. 

In the experiments for NpTGas, C-, A-, and G-AF magnetic phases have been 
found in NpFeGas, NpCoGas, and NpNiGas. Here we have a naive question: What 
is a key parameter to understand the change of the magnetic structure? In the case 
of UTGas, it has been claimed that the level splitting A is important to explain the 
difference in magnetic structure as well as the magnetic anisotropy for a fixed value of 
n=l. Roughly speaking, A is positive for T=Fe, small positive for T=Co, and negative 
for T=Ni. Among UTGas with T=Ni, Pd, and Ft, when we assume that the absolute 
value of A is increased in the order of Ni, Pd, and Pt, it is possible to understand 
qualitatively the change in the magnetic anisotropy, in addition to the change in the 
magnetic structure of G-AF for T=Ni and A-AF for T=Pd and Pt. It has been found 
that the value of t^ is not so crucial to explain qualitatively the magnetic properties of 
U-115 based on the two-orbital model for n=l. 

For n=2, we always obtain the G-AF phase. However, for n=1.5, we have observed 
three kinds of AF magnetic structure in the phase diagrams. Let us summarize the 
change in the magnetic structure for a fixed value of ^2=0.8. Note that this value is 
larger than ^^=0.1, which we have considered to reproduce two kinds of cylindrical 
Fermi-surface sheets of Np-115. However, in the small-sized cluster calculations, it is 
difficult to compare directly with the values in the thermodynamic limit. Thus, we do 
not discuss further the quantitative point on the values of t^. As shown in Fig. IHM b). 
for J=0 and tz=0.8, we see the change in the magnetic structure as G-AF (A<0), C- 
AF(0<A<0.4), and G-AF (A>0.4). For J=0.5 and tz=0.8, as shown in Fig.EH^c), the 
C-AF phases are always observed, but they have different orbital structures. Finally, for 
J=4 and ^^=0.8, we observe C-AF (A < -0.15), A-AF(-0.15 < A < 0.3), and C-AF 
(A>0.3), as shown in Fig. IHHTd). 

In order to understand the appearance of three types of the AF phases, we may 
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consider an explanation due to the combination of the changes in A and n. For instance, 
by assuming that J=4 for NpTGas and the change in A for NpTGas is just the same 
as that for UTGa5, we consider that n^2 with A<0 for T=Ni, n~1.5 with AfaO for 
T=Co, and n~1.5 with A> for T=Fe. Then, it seems to be possible to relate our 
theoretical AF phases with the experimental observations in NpTGas. However, it is 
difficult to claim that the above parameter assignment for three Np-115 materials is the 
best explanation for the magnetic structure of Np-115, since in actual compounds, there 
are other important ingredients which have not been included in the present model. For 
instance, we have never discussed the direction of the magnetic moment of Np ion. In 
particular, the canted AF structure cannot be considered at all for the G-AF phase of 
NpNiGas. Thus, we need to recognize some distance between the actual magnetic states 
and the theoretically obtained phases. Our theory should be improved by taking into 
account other realistic ingredients of 115 structure. 

Finally, let us remark a possible future direction of the research. In Fig. IHUT b). we 
have explained the G-AF (II) as the orbital disordered phase, which is considered to be 
related to the metallic phase. Of course, we cannot conclude the matallicity only from 
the present calculations for small-size cluster. However, it seems to be an interesting 
concept that the /-electron state is controlled by orbital degeneracy. Quite recently, 
Onishi and Hotta has pointed out the orbital incommensurate state appearing between 
two kinds of localized AF phases, by using the density matrix renormalization group 
method to the model of the j-j coupling scheme |266j . We may throw new lights on the 
long-standing issue concerning the competition between localized and itinerant nature 
of /-electrons, when such competition is controlled by orbital degree of freedom. 

6.2. Multipole ordering 

In the previous subsection, we have discussed the spin and orbital structure of 115 
materials on the basis of the Fg model. Again we note that such "spin" and "orbital" 
are not real spin and orbital, but pseudo spin and orbital. In principle, in /-electron 
systems, real spin and orbital are not independent degrees of freedom, since they are 
tightly coupled with each other due to the strong spin-orbit interaction. Then, in order 
to describe such a complicated spin-orbital coupled system, it is rather appropriate to 
represent the /-electron state in terms of "multipole" degree of freedom, rather than 
using spin and orbital degrees of freedom as in (i-electron systems. We show a table 
to summarize the multipole operators up to rank 3. The relation of the irreducible 
representation with pseudo spin and orbital are also shown. When we use the term of 
multipole degree of freedom, pseudo spin a and pseudo orbital r are r4u dipole and 
Fag quadrupole, respectively, from Table 3. Thus, in the previous section, we have 
considered the ordering tendencies of dipole and quadrupole degree of freedom. Note 
that it is important to show explicitly the parity of the multipole operator to complete 
its symmetry. In this section, for convenience, we use the subscripts "u" and "g" for 
odd and even parity, respectively. 
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Table 3. Multipole operators in the Fg subspace, shown in Ref. |2f)7| . The 
first, second, and third lines denote 7 of the irreducible representation T-y, multipole 
operator X^"! , and pseudospin representation, respectively. The multipole operators 
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In this article, we pick up octupole ordering in Np02. First we briefly discuss the 
level scheme for actinide dioxides with CaF2 cubic crystal structure. Due to the CEF 
effect, the sextet is split into Fg quartet and r7 doublet. In this case, due to the intuitive 
discussion on the direction of the extension of orbital and the position of oxygen ions, the 
Fy state should be higher than the Fg level. Here we define the splitting energy as A. As 
discussed in Sec. 5, in order to make the low-spin state, we accommodate two, three, and 
four electrons in the Fg level. Then, the ground states are F5, Fg , and Fi, respectively 



[12], consistent with the CEF ground states of UO2 EHU M, NpOz ^M I2ZI1, and 
Pu02 |272|I27I-?] . respectively. Then, A is estimated from the CEF excitation energy in 
Pu02, experimentally found to be 123 meV |272[ 1273] . On the other hand, as mentioned 
repeatedly in this article, the Hund's rule coupling Jh between Fg and F7 levels is 1/49 
of the original Hund's rule interaction among / orbitals. Namely, Jh is as large as a few 
hundred Kelvins. Thus, we validate our assumption that the F7 state is simply ignored. 
From the qualitative viewpoint, unfortunately, this simplification is not appropriate 
to reproduce experimental results, since the ground-state wave-function is not exactly 
reproduced in the Fg model. However, we believe that this approximation provides a 
qualitatively correct approach, in order to understand the complex multipole state from 
the microscopic viewpoint. 

Let us here review the recent theoretical result on octupole order of Np02 
|113| 1116] I117t ITTH] I119J . We set the Hamiltonian for actinide dioxides as the Fg model 
Eq. fll52j) on an fee lattice. The form of the Hamiltonian is already given, but the 
hopping should be estimated on the fee lattice. For instance, the hopping amplitudes 
between /-orbitals at (0, 0, 0) and (a/2, a/2, 0) (a is the lattice constant) are given by 

4 2VSi 
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(158) 
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Table 4. Coupling constants in the effective model. The energy unit is (l/16)t^/([/' — 
J). 
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where t = (//o")/28 and t~^.^,^, =t^^.^,^,. Note that the hopping integrals depend on 
/x and they are intrinsically complex numbers in the fee lattice. It is in sharp contrast 
to the previous model on a simple cubic lattice, in which the hopping integrals are real 
and the same form as for Cg orbitals of d electrons. In other words, there is no difference 
between the Fg model for / electrons and the Cg orbital model for d electrons on the 
simple cubic lattice in the absence of a magnetic field. For Cg orbitals, we can always set 
the hopping integrals to be real irrespective of the lattice type, by selecting appropriate 
basis wave-functions, while Fg orbitals on the fee lattice appear to be complex in nature, 
specific to /-electron systems with strong spin-orbit coupling. 

In order to discuss multipole ordering, we derive an effective multipole model in 
the strong-coupling limit using standard second-order perturbation theory with respect 
to t, which is exactly the same as that used in the derivation of the Heisenberg model 
from the Hubbard model in the research field of transition metal oxides. It is one of 
advantages of the /-electron model on the basis of the j-j coupling scheme that we can 
exploit the technique developed in the d-electron research. We consider the case of one 
electron per / ion in the Fg orbitals, but the effective model is the same for the one-hole 
case, which corresponds to Np02, due to an electron- hole transformation. Among the 
intermediate /^-states in the perturbation theory, we consider only the lowest-energy 
F5 triplet states, in which the two electrons occupy different orbitals, assuming that 
other excited states are well separated from the p ground states. In fact, the excitation 
energy from the F5 ground state of p in UO2 is 152 meV |268| I269J . Note that the 
CEF excitation energy is considered to be larger than the triplet excitation one, since 
the Hund's rule interaction is effectively reduced in the j-j coupling scheme. Thus, it is 
reasonable to take only the F5 states as the intermediate states. 

After straightforward, but tedious calculations |I113|I116] . we arrive at an effective 
model in the form of 

-f^eff = /_^{Hlq + H2q + IH-Aulq + -f^4«2q), (160) 

q 

where q is the wave vector. Hiq denotes the interactions between quadrupole moments, 
given by 

Hiq = a\{0\^qOlqC^Cy -h c.p.) 
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Figure 37. (a) The 8-site fee eluster, shown by sohd spheres, taken in the ealeulation 
of exaet diagonaHzation. (b) Correlation functions for the 8-site cluster for q — (0, 0, 0) 
(triangles), q — (0,0, 1) (squares), and q = (1/2, 1/2, 1/2) (diamonds) in units of 27r/a. 



+ asiO^^^qO^y^qS^Sy + C.p.) 

+ ClAiOxy-qOxy^qCxCy + C.p.), (161) 

where c.p. denotes cyclic permutations, Ci, = cos{qua/2), and s^ = sm{q^a/2) {u = x, y, 
or z). The definitions of the multipole operators and values of the coupling constants Oj 
are given in Tables El and El respectively. Note that 0°^ transforms to ("\/30|q — 02g)/2 



and {—y/^Ol — 02„)/2 under c.p. {x,y, 



{y,z,x) and {x,y,z) 



{z,x,y), 



respectively. T-C2q and H^^nq {n = 1 or 2) are the interactions between dipole and 
octupole moments, given by 

H2q = 68[Tf_,Tf^(c,c, + c,c,) + c.p.] 
+ 69[T,%T^^>.s, + c.p.] 

+ bioTxyz,^qTxyz,q{CxCy + c.p.), (162) 



and 



K 



iunq 



+ 4" 

in 



+ K 

where values of the coup 



[J'-Vtrc^cy + c.p.] 



4un jAun , 



[jT-qJtTs^Sy + c.p.] 
['''Xyz-q\'J zq -^x-^y + C.p.jJ 



(163) 



\TtqJtTc.{c^-Cy) + c.v.)] 

ling constants hi and h^ are shown in Table El The above 

Eqs. ()160p - ()163|) are consistent with the general form of multipole interactions on the fee 

lattice derived by Sakai et al. |274j . We follow the notation in Ref. |274j for convenience. 

When a mean-field theory is applied to the effective model, due care should be 

taken, since in an fee lattice with geometrical frustration, the effect of fluctuations may 
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Figure 38. The triple-q r5„ octupole state in the fee lattice. 



be strong enough to destroy the state obtained within the mean-field theory. Thus, we 
first evaluate the correlation function in the ground state using an unbiased method such 
as exact diagonalization on the A^-site lattice. Here we set A^ = 8, as shown Fig. ITTT a). 
The correlation function of the multipole operators is defined by 



A.C 



(l/iV)^e^''-(^-^')(X^^X;;7) 



(164) 



where (■ ■ ■) denotes the expectation value by the ground-state wave- function. 
Figure ITTT b) exhibits the results for correlation functions. As shown in the Table 
4, the absolute value of the interaction between T2u moments (6io) is the largest 
among multipole interactions, but the correlation function of the r2u moment is not so 
enhanced, suggesting that the frustration effect is significant for an Ising-like moment 
such as T2u- Rather, large values of correlation functions are found for J^"^; '^z^^ ^^^ O^y 
moments at q = (0, 0, 1) in units of 211 /a. Note that the effective model does not include 
the term which stabilizes O^y quadrupole order at q = (0,0, 1). The enhancement of 
this correlation function is due to an induced quadrupole moment in Tiu2 or T^u moment 
ordered states. Thus, the relevant interactions of the system should be 62 and b^, which 

stabilize the J|"^ and T|" order, respectively, aX q= (0, 0, 1). In the following, then, we 

(2) 
consider a simplified multipole model which includes only b^ and b%. 

Now we apply mean-field theory to the simplified model to specify the ordered state. 

As easily understood, the coupling constant fog is slightly larger than 62 , indicating that 

Fsu order has lower energy than F4„2 order. The interaction fog stabilizes longitudinal 

ordering of the T^u moments, but their directions are not entirely determined by the 

form of the interaction. Here we point out that in the Fg subspace, the T^^ moment 

has an easy axis along [111] |100| 11011 IllOj . Thus, taking the moment at each site 

along [111] or other equivalent directions, we find that a triple-q state is favored, since 

it gains interaction energy in all the directions. In fact, as shown in Fig. EHl the 

ground state has longitudinal triple-g Fs^ octupole order with four sublattices, i.e., 

((T|;f), (Ty^^), (T|^)) oc (exp[i27ra;/a],exp[i27r|//a],exp[i27rz/a]). Note that this triple-g 
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structure does not have frustration even in the fee lattice. The ground state energy 
is —468 per site, and the transition temperature Tq is given by fce^o = 468- Another 
important message of Fig. IHHl is that both up- and down-spin densities are anisotropic 
with different distribution. It may be natural, if we consider octupole as the combined 
degree of freedom of spin and orbital, but such an intuitive explanation of octupole 
becomes possible from the microscopic viewpoint on the basis of the j-j coupling scheme. 

Let us briefly introduce new experimental evidence of octupole ordering in Np02. 
Quite recently, in order to clarify the nature of the ordered phase of Np02, Tokunaga 
et al. have performed ^^0-NMR measurement of Np02 below Tq |lllj . Tokunaga et 
al. have observed the occurrence of two inequivalent oxygen sites below Tq from the 
^^0-NMR spectrum. They have also found that the hyperfine interaction at the oxygen 
sites are explained by invoking a hyperfine interaction with field-induced AF moments 
due to the longitudinal triple-g octupole order. Thus, the NMR results strongly support 
the occurrence of the longitudinal triple-q multipole structure in Np02. 

Here we briefly mention the physical properties of the phase in the high-temperature 
region. The magnetic susceptibility of UO2 follows the standard Curie- Weiss law |275j . 
while that of Np02 is significantly deviated from the Curie- Weiss behavior well above 
the transition temperature Tq |276|I277] . As pointed out by Kubo and Hotta |278j . the 
difference in the temperature dependence of magnetic susceptibilities in these materials 
is naturally explained due to the fact that dipole and octupole moments coexist in Np02 
while only dipole exists in UO2. When the magnetic moment consists of two independent 
moments, it is intuitively understood that in such a situation, the magnetic susceptibility 
is given by the sum of two different Curie- Weiss relations, leading to non-Curie- Weiss 
behavior. It is one of characteristic issues of the system with active orbital degree of 
freedom. 

We have not included the effect of oxygen in the model, but it has been already 
shown that the conclusion does not change |117j . even if we analyze the so-called f-p 
model which include both 5/ electrons of Np ion and 2p electrons of O ions. It is also 
possible to perform similar analysis for another lattice structure. In fact, we have found 
a r^g antiferro-quadrupole transition for the simple cubic lattice and a r2„ antiferro- 
octupole transition for the bcc lattice |116|I118] . 

7. Conclusions 

We have reviewed orbital ordering phenomena in d- and /-electron compounds starting 
from a basic level for the construction of the model. Since this subject includes so many 
kinds of materials such as transition metal oxides, rare-earth compounds, and actinide 
materials, it is almost impossible to cover all the results concerning the orbital-related 
phenomena in these compounds. Thus, we have picked up some typical materials and 
attempted to explain how and why the orbital ordering occurs. 

In order to summarize this article, we would like to emphasize three points. One 
is the understanding of the orbital ordering from a band-insulating picture, which 
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beautifully explains the appearance of E- and CE-type spin structure. The heart 
of the explanation is the interference effect of electron phase originating from the 
anisotropic orbital. By including further the realistic effect of Jahn- Teller distortions 
and/or Coulomb interactions, it is possible to obtain spin, charge, and orbital ordering 
of transition metal oxides. 

Second is the similarity between d- and /-electron orbital, when we consider the /- 
electron model on the basis of the j-j coupling scheme. In particular, we have remarked 
significant correspondence between Cg orbital degenerate model for d electrons and Eg 
model for / electrons. The large variety of manganite-like magnetic structure observed 
in U-115 and Np-115 has been understood qualitatively on the basis of the Eg model. 

The third point is the understanding of multipole order in a microscopic level on the 
basis of the j-j coupling model. As a typical example, we have reviewed the octupole 
ordering in Np02. By applying the (i-electron-like analysis on the Eg model on the fee 
lattice, it is possible to understand naturally the stability of octupole ordering. It is 
also remarkable that octupole is clearly understood by the anisotropic spin-dependent 
charge distribution. 

Of course, there still remain more kinds of orbital ordering in d- and /-electron 
compounds, which cannot be explained in this review article. We could not cite lots of 
important papers of other authors on the issue of orbital ordering phenomena. However, 
it is not the main purpose of this article to introduce orbital order in d- and /-electron 
systems with complete references. We would like to convey the unique viewpoint that 
orbital ordering (including multipole ordering) is the common phenomenon in d- and /- 
electron systems. By developing further microscopic theory on orbital ordering using the 
orbital degenerate model, we hope that it is possible to understand complex magnetism 
among transition metal oxides, rare-earth compounds, and actinide materials, from a 
unified viewpoint. 
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